A Final Report 
Grant No. NAG-3-909 

June 1, 1988 - April 17, 1992 


>r 

UT\ 



& 

\ O 




\ vn 


o 

r- 

3 r- 


p- 

rvi 

O C fM 

t V 

o 

i 

X 1 


rsi 

m 

K *> 

o 

sO 

0 s 

\ O 

c 

r-t 

z 

\ z 

o 

o 




O' 




fO 








ro 




o 


f* c 


MAGNETIC ACTUATORS AND SUSPENSION FOR 
SPACE VIBRATION CONTROL 

Submitted to: 

National Aeronautics and Space Administration 
Lewis Research Center 
21000 Brookpark Road 
Cleveland, OH 44135 

Attention: 

Dr. David P. Fleming, M/S 23-3 
Structural Dynamics Branch 

Submitted by: 


Carl R. Knospe 
Assistant Professor 


r-t Q. 

I *1 




Paul E. Allaire 

Mac Wade Professor and Chairman 

David W. Lewis 
Professor 

SEAS Report No. UVA/528292/MANE93/101 
May 1993 


> 


DEPARTMENT OF MECHANICAL, AEROSPACE 
AND NUCLEAR ENGINEERING 



The Journal of the Astronautical Sciences, Vol. 40, No. 2, April-June 1992, pp. 241-259 


/.yv i 
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Introduction 1 


R.D. Hampton , 2 C. M. Grodsinsky , 3 P.E. Allaire , 2 D.W. Lewis , 2 

and C. R. Knospe 2 


Abstract 

Certain experiments contemplated for space platforms must be isolated from the a«el- 
’rations of the platform. In this paper an optimal active control is developed for micro- 
gravity vibration isolation, using constant state feedback gains (identical t0 ' hose obtamed 
from the Linear Quadratic Regulator [LQR] approach) along with constant feedforward 

cosl function fo, .his control algorithm effectively weights accelerations 
of the platform due to external disturbances by a factor proportional to (1/m) ■ Low fre- 
quency accelerations (less than 50 Hz) are attenuated by greater 

tude The control relies on the absolute position and velocity feedback of the experiment 
and the absolute position and velocity feed-forward of the platform, and generally derives 
the stability robustness characteristics guaranteed by the LQR approach to optimality. 

The method as derived is extendable to the case in which only the relative positions 
and the absolute accelerations of the experiment and space platform are available. 

Introduction 

A space platform experiences local, low frequency accelerations (0.01-30 Hz) 
due to equipment motions and vibrations, and to crew activity [1, 2], as indicated 
in Fig. 1 [3]. Certain experiments, such as the growth of isotropic crystals, re- 
quire an environment in which the accelerations amount to only a few micro-gs 
f41 fSee Fie 2 T21.) Microgravity requirements to accommodate such experi- 
ments have been specified for Space Station Freedom (SSF). (Fig. 1.) [2] Such an 
environment is not presently available on manned space platforms [2J. 

'Presented at the Workshop on Aerospace Applications of Magnetic Suspension Technology, 
September 25-27, 1990, under the title "Microgravity Vibration Isolation: An Optimal Control 

J Depanment oTMMhan'ical and Aerospace Engineering, University of Virginia, Charlottesville, 
VA 22901. 

'NASA Lewis Research Center, Cleveland. OH 44135. 
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FIG. 1. SSF Microgravity Requirements and Anticipated Acceleration Environment. 
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Since the experiment and space platform centers of gravity do not coincide, a 
means is needed to prevent the experiment from drifting into its own orbital 
motion and into the space platform wall. Additionally, some experiments require 
umbilicals to provide power, experiment control, coolant flow, communications 
linkage, or other services. Unfortunately, such measures also mean that unwanted 
platform accelerations will be transmitted to the experiments. This necessitates 
experiment isolation. Passive isolators, however, cannot compensate for umbilical 
stiffness, nor deal adequately with direct disturbances, nor can they achieve low 
enough corner frequencies even if umbilicals are absent [5]. Active isolation is 

therefore essential. . . . 

The problem, then, is to design an active isolation system to minimize these 

undesired acceleration transmissions for a tethered payload, while achieving ade- 
quate stability margins and system robustness. Spatial and control energy limita- 
tions must also be accommodated. Although microgravity isolation systems have 
been developed and tested [6], no controller offered to date takes into account 
the effect of umbilicals in isolator control design [7]. 

Mathematical Model 

The general problem has three translational and three rotational degrees of 
freedom For simplicity, however, this analysis will consider only the one-dimen- 
sional problem. The general problem could be treated in an analogous manner. 
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FIG. 2. Ex 
a Variety of 


Deriment Isolation Requirements (Tolerable g-level as a Function of Frequency for 
Materials Science Experiments as Predicted by Order-oEMagnitude t Analysts for a 
Single-Frequency Disturbance [after Demel, 1986, from Nelson, 199 ]). 


Let the experiment be modeled as a mass m, with position x(t). Assume that the 
space station has position </«. and that umbilieals with stiffness * and damping 
c connect the experiment and space station. Suppose further that a magnetic ac- 
tuator applies a control force proportional to the applied current i(f), with pro- 
portionality constant a. Such a model is shown in Fig. 3. 

The system equation of motion is 
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Division by m and rearrangement yields 

x = — -{x — d) — ~ d) 

m m 

In state space notation this becomes 

x = Ax + bti + f 


a . 
m 


( 2 ) 


(3) 


where 


x = 


A = 


H, 

W W hi 




A77 


b = 


a 


m 


u = i, 


f = 


0 


k . c • 
— d + — d 
m m 


The objective is to minimize the acceleration i:(r). 

Optimal Control Problem 

The optimal control problem is that of determining the control current u(f) = 
i(t) which minimizes a suitable performance index 

J = J(x,u,t) (4) 

for the system described by equation (3) subject to the state variable conditions 

x(0) = Xo (5a) 

lim x(/) = 0 

Another reasonable condition is that f(t) is bounded and it will be found mathe- 
matically advantageous to assume that f(r) is also a dwindling func i 

lim f(r) = 0 < 5c) 

In actuality, f(t) is not dwindling^and so neither is x(f). However for bounded 
f(r) (and a controllable system [8] such as this), the optimal contro ler developed 
in this paper has only a vanishing dependence on f(r) for times in the dista 
lure and't depend/on x(0 only in a causa, fashion. Further, with only minor 
changes in theperformance index [9] (and more complicated mathematics) the 
dwindling assumptions can be removed without affecting t e resu t,ng ’ 
using the § basic method of the present paper. [Tomizukas dynamic programm g 
approach to the command-following problem leads to corresponding results [1 ].] 
Hence, the above simplifying assumptions are justifie . 

A quadratic performance index 


J = — [ [x r IT,x + w>u z ]dt 
2 Jo 


( 6 ) 
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, , thit ipnrlc itself well to the variational approach to opti- 

has been chosen, as one that lends itseii weii to i ^ defi . 

mat controls since an analytical solution is desired. The upper lirni of the de 
nite integral has been selected so as to yield a time-invariant controller. Here W x 
is a square 2 x 2 constant weighting matrix while w, is a weighting const ^ n ^ 
Although W\ could be a full 2 x 2 matrix, for this problem a diagonal form 

has been employed for the sake of simplicity. 

\w Xa 0 
0 


w t = 


Wu, 


(7) 


The performance index consequently reduces to 


/ = — f [w u x f + w'u.jc; + W}U 2 ]dt, 

2 J{| 


( 8 ) 


so that each state is weighted independently. 

If sinusoidal motion of the experiment is considered, so tha 

x(f) = B , sin wt 

and x(f) = -w 2 x(t), the cost function can be expressed in terms of the accelera- 
tion and control as 

(9) 


.i r[(^ + ^W 

2 J tl \oj 


x 2 + WjU 2 


dt 


It is apparent that this performance index conveniently weights accelerations at 
low frequencies much more than at higher frequencies. 

Solution 

Finding the optimal control to minimize equation (4) is a variational problem 
of Lagrange for which the initial steps of the solution are well-known (e g 
Elbert [11]) The variational approach is outlined below, following w ic 
S eals added b, .he nonhomogeneous term «.) »■« * M***-0£ 
rent optimal control texts either assume that f« 0 * • ' that 

auire that it have a restricted range space (e.g., [12], p. 238). 
follows provides an analytical optimal without imposing such restriction . 
TSS v the “St function 7 from equation (6) is augmented b, the 
Lagrange multiplier A times the system equation of motion (3) where 

(10) 


_ A i 

l A 2, 


The result J can be expressed as 


where the Hamiltonian H is 


-f 


Hdt 


H = Tx r M / ,x + w,if) + A r (x - Ax bu f) 

It is desired to obtain an optimal solution u = u* which minimizes J. 


HD 


( 12 ) 
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A. , 

The first variation of J {\ , u, x) is 

1 \dH 


8J 


= r + 

J 0 L 5x du dx J 

which is set equal to zero to minimize J. However, integrating by parts, 
J |^r6xjdf = — | (\ T 8x)dt 

A 

so that the above expression for 8J becomes 


8J 


-r[g-') 


aw f 

8x + “ 5 m 

du 


dt = 0 


Both 8x and 8u are arbitrary variations, so 5/ - 0 only if 

dx 

dU 

The conditions given by equation (5) still apply. 

Solving equations (14a) and (14b) yields 

A = W\\ - AX 
u* = — A r b = b r A 

w 3 VV 3 

Temporarily eliminating u* produces the result 

{a} - 4*1 + W 

where 


A = 


(13) 

(14a) 

(14b) 

(15a) 

(15b) 

(16) 


/4 

-bb T 

Wj 

W x 

-A 1 


If equation (16) is now solved for A in terms of x and f, equation (15 ) w 

furnish an expression for the optimal control u*. 

As noted before, optimal control texts generally treat the homogeneous prob 
lem (where f(f) s 0), but they do not provide an analytical solution to the non- 
system described b, equations (5) and (16). Salukvadze has treated 
the nonhomogeneous problem [13,14], but his difficult fmatment seems Uigely lo 
have remained either uncomprehended or under-appreciated, l “ 

pecially well-suited to low-frequency disturbance rejection, and has been app 
below to the present problem. 
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(17) 


The four et 


igenvalues of A may be found to be, in ascending order of real parts. 


= 


M 2 = 


- fl. + (ff? - 4/3 i ) i,2 \ 1/2 


Zjx - (gj - 

2 


r - 4/3;)' y 7 


JU.3 = 

p,4 = “M2 

where 0 i and j 3 2 are defined as follows: 

2k _ _ Wit 

~ 7n m 2 tmv 3 


(18a) 

(18b) 

(18c) 

(18d) 

(19a) 


and 


The eigenvectors 
chosen to be 


, Ja-wu k\ (19b) 

* = rrv) 

of A corresponding to the respective eigenvalues M* may be 


P k = 


74 


y? Yi(72 + M-*) 

flic 73^* ?3 

Yi + (72 + 

73 


( 20 a) 


Using equations (18) 


and 74 are defined below: 


k 

y\ = TT 

/71 

( 20 b) 

C 

yi= m 

( 20 c) 

a 2 

( 20 d) 

m : w 3 


y 4 = vv la 

( 20 e) 

(18) through ( 20 ) with equation 

(17) the solution to the homoge- 
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neous system is 



Cie MI 'pi, + c 2 e M; 'p 2 , + c 3 e Ml 'p 3 , + c 4 e M; 'p4, 
c 1 e' 11 'pi, 4- c 2 e M; 'p 2 , + c 3 e"' l ' , p 3 , + c 4 e'“ ; 'p 4 . 


( 21 ) 


with p* = {p* - , p* 7 ,} 7 ", k = 1 4 and where are arbitrary constants. 

Application of the variation of parameters method with the terminal con- 
ditions (5b, c) leads to the general solution of the non-homogeneous system, with 
two constants of integration yet undetermined. 

If the two constants of integration are eliminated by solving for A in terms of 
x and f, the general solutions for A, and A 2 become: 

A, = £,jc, + f 2 Jr 2 + f 3 e _, “' + f.e'" 2 ' (22a) 

A 2 = fjJfi + £>*2 + + ^e^-’ 1 (22b) 

in which the £,’s are functions of the eigenvalues and eigenvectors of A, and of 
the disturbance f(t). 


The Solution Form 


Using the fact that 


«*(,) = — b r A [cf. (15b)] 

W 3 


(23) 


the optimal control is found to be 

u*(t) = tj.x, + T 72* 2 + 7 + ij.e - '* 2 ' | e M; '/ 2 (r) dt (24a) 


where 


-ml k \ 

V\ = Mi M: 

a \rn J 

- m / c \ 

Vi ~ — + Mi + M 2 1 

a \m } 

a \fi x - in) \ m mf 

a \fi i - M2/ \ m m / 


Vi 


V 4 


(24b) 

(24c) 

(24d) 

(24e) 


(It should be noted that the feedback gains 77 , and tj 2 are those which would re- 
sult from applying standard LQR theory to the homogeneous system equation 
x = Ax + bul) In equation (24) and m are the eigenvalues of A with nega- 
tive real parts, [see equations (18a, b)] and 


MD = —d + —d. 

mm 


(24f) 


For = nz = n [the critically damped case] equations (24a-e) reduce to 
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U*(t ) = -("tM 2 ~ k )x i + -C 2m ^ C) 

a a 


*2 


- e J < 




no* 


By repeated application of the method of integration by parts, the control may 
be re-expressed in terms of an infinite sum. 


/ * (-l) r mO \ f y (-l) r CT\ 

= T7 1 + r\iXi + 1)3^2 ^r¥\ I ^ 4 \ r=0 ix r i x / 

Rewriting fr in terms of d and d, the control function becomes 


(25) 


n-l 

+ 2 

< = l 


(-D 




(-i r 


m 


m \fi\ n) m V#*» ^ I 

\m" Pi / 


d { "\t) + higher order terms 


(26) 


(27) 


This may be written in a more appealing form as 

„•(,) = c„x(t) + c,x(t) + C d0 d(t) + c M m + higher order terms 

in which the constant coefficients c p , c„, c M , and c d i may be defined from equa 
tions (24) and (26). Clearly, if the infinite sums converge rapidly enough, P 
timal control can be approximated by 


u*(t) = c p x{t) + c„i(f) + c J0 d(t) + c dl d(t) 


(28) 


For very low frequency disturbances the higher order terms in equation (26) are 

Control Evaluation 

Physical Realizability of the Control 

tie ^ ,* -*« - *■*«<»> 
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TABLE 1. Optimal F/B and F/F Gains for Selected State Variable and Control Weightings 
System Parameters: m = 100 Ibm, c = 0.000622 Ibf-sec/ft = 0.1%), k = 0.3 Ibf/ft, 
a * 10 Ibf/amp 


Weights 

F/B Gains 

F/F Gains 

w la 

w lb 

w 3 

c p 


C d 0 






2 

1 

1 

1.3845 

1.3637 

0.0294 

- 0.0006 

- 0.0070 

- 0.0067 

- 0.0049 

- 0.0032 

10 

2 

1 

3. 1324 

1.9863 

0.0297 

- 0.0001 

- 0.0030 

- 0.0019 

- 0.0009 

- 0.0004 

23 

3 

1 

4 . 7659 

2.4413 

0.0298 

- 0.0000 

- 0.0020 

- 0.0010 

- 0.0004 

- 0.0001 

41 

4 

1 

6.3732 


0.0299 

0.0000 

- 0.0015 

- 0.0007 

- 0.0002 

- 0.0001 

64 

5 

1 

7.9701 

3.1544 

0.0299 

0.0000 

- 0.0012 

- 0.0005 

- 0.0001 

- 0.0000 

92 

6 

1 

9 . 5617 

3.4552 

0.0299 

0.0000 

- 0.0010 

- 0.0004 

- 0.0001 

- 0.0000 

126 

< 

1 

11.1950 

3.7354 

0.0299 

0.0000 

- 0 . 0008 

- 0.0003 

- 0.0001 

- 0.0000 

165 

8 

1 

12.8153 

3.9949 

0.0299 

0.0000 

- 0.0007 

- 0.0002 

- 0.0001 

- 0.0000 

209 

9 

1 

14.4269 

4.2380 

0.0299 

0 . 0000 

- 0.0006 

- 0.0002 

- 0.0000 

- 0.0000 

258 

10 

l 

16.0324 

4.4674 

0.0299 

0.0000 

- 0.0006 

- 0.0002 

- 0.0000 

- 0.0000 

581 

15 

l 

24.0740 

5.4729 

0.0300 

0.0001 

- 0.0004 

- 0.0001 

- 0.0000 

- 0.0000 

1034 

20 

1 

32.1259 

6.3209 

0.0300 

0.0001 

- 0.0003 

- 0.0001 

- 0.0000 

- 0.0000 

1617 

25 

1 

40.1819 

7.0680 

0.0300 

0.0001 

- 0.0002 

- 0.0000 

- 0.0000 

- 0.0000 

2329 

30 

1 

48.2297 

7.7431 

0.0300 

0.0001 

- 0.0002 

- 0.0000 

- 0.0000 

- 0.0000 

3171 

35 

1 

56.2816 

8.3640 

0.0300 

0.0001 

- 0.0002 

- 0.0000 

- 0.0000 

- 0.0000 

4143 

40 i 

1 

64.3361 

8.9420 

0.0300 

0.0001 

- 0.0001 

- 0.0000 

- 0.0000 

- 0.0000 

9325 

60 

1 

96.5360 

10.9526 

0.0300 

0.0001 

- 0.0001 

- 0.0000 

- 0 . 0000 

- 0.0000 


which the closed loop system is closest to being critically damped without being 
underdamped. Corresponding controller feedback and feed-forward gains (for 
the first five derivatives) are also included. 

The states x{t) and i(/) and the derivatives d m {t), d a \t ), and d n \t) are clearly 
available for an Earth-based system. However, in space, the only absolute mea- 
surements which can be directly available ar cx(t) and from which Jc(r), d(t) 

and x{t), d(t) are obtainable only by successive integration(s). Rearrangement of 
equation (28) into 

u*(t) = (Cp + c d0 )x(t) + (c v + c di )x(‘) ~ c M [x(t) - d(t)] - Crfi[i(r) - d{t)} (29) 

or 

u*(t) = (c p + c M )d(t) + (c v + C d] )d(t) + Cp[x(t ) - d{t)] + c v [i(f) - d(0] (30) 

obviates the need for one accelerometer, but one accelerometer plus two inte- 
grations remain necessary for either the platform or the experiment. Since 
[jc(t) - d{t)} (or one of its integrals) has not been weighted in the performance 
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index J, experiment drift will be a problem that must be corrected either by an- 
other control loop or by a change of system states^ The latter could be accom- 
plished by incorporating an accelerometer attached to the experiment .into 
state equation. Alternatively, one could append an integrator to the : p ant in- 
clude the current i(t) as a third state, and optimize the control di/dt. But for the 
sake of simplicity (i.e., fewer states) the former has been assumed (without devel- 

^The'higher'o^der 6 terms of the control [equations (25) and J 26) ] 
elected for low frequencies, if the eigenvalues and are of sufficient mod 
Ulus These eigenvalues, in turn, are under the control of the designer 
determined by his choice of weights w la , w,>, and w 3; It is apparent from 
equation (25) that u*(t) essentially reduces to two alternating power series. For a 
shi^soVdal disturbance of frequency „ the series form of the control converges for 
Tl <10 = 1,2). It can be shown that each alternating power sene con- 
verges like 2* 0 (-l)Wu) 2 '. with ‘ t|ow ” frequency disturbances (i.e., small rela 
live 8 to s'ysVem closed loop eigenvalues) a control formed by senes truncation very 

'^For^TampklTuppos^ JhatThe normalized frequencies M for a sinusoidal 
disturbance are less than 1/5, and that only the feedforward control tenns^J 
and cnd(t) are included with the feedback terms. Even so, the feedforward por 

tion of the truncated control, at any time f, will be a current that he' normal” 

4% [i e , (1/5) 2 ] of the feedforward portion of the actual optimal. If the nor 
ized 'frequencies are below 1/10, this approximation error will be _ less . than U 
Table 1 shows that the gains c* of higher order derivatives d (t) [see equa 
tion (26) for algebraic representations] are, in fact, quite small. 

In some circumstances there may be design constraints which prevent the 
siener from selecting weights that will lead to sufficiently rapid convergence. 
Howevlr^ ^convergence occurs rapidly even for eigenvalues of relatively 
small modulus (|o»/ M ,| < 1/3), in a great many cases the designer wtl [ have muc 
latitude in his choice of weights. For “low’ frequency disturbances in these 
cases, a control which includes only one or two feedforward terms will be close 
to the optimal. These frequencies will be well-attenuated. 

Higher frequency disturbances will also be well-attenuated, provided the 
input-to-output transfer function(s) are proper in the Laplace Transform van 
5 This will not be the case for the present problem if more than three feedfor- 
ward gains (c.c.c) are included in the control. Practically ^ ^ ^ 
only proportional and first-derivative feedforward [equation (25) with r - 0,1 or 
eauauon P (26) with n = 2] should be added to the feedback control terms. As will 
how.v«! adding even .he proportiona, feedforward 
dramatically improve the disturbance rejection over that afforded by LQR feed 
back alone. 

Transfer Function and Block Diagram 

Neglecting the higher order terms, the transfer function between input and 
output accelerations or displacements is 
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s^X(s) 

s 


X(s) 


(I - * (« : 


CdO 


(31) 


= DW - <*» = + (i + 4 ♦ (i + c) 

and a block diagram of the controlled system can be drawn as in Fig. 4. 

Control Stability, Stability Robustness, and General Robustness 

Since the control feedback gains are the same as those obtained by solution of 
the standard Linear Quadratic Regulator (LQR) problem, the closed loop system 
is stable and enjoys the stability robustness characteristics guaranteed by the 
(LQR) approach to optimality, viz., a minimum of 60° phase margin, in mite 
positive gain margin, and 6 dB negative gain margin [12], Additionally, numeri- 
cal checks indicate that it enjoys substantial insensitivity, or genera ro ustness 
to uncertainties in k, c, and m, as indicated by Table 2 and Figs. 5 through 12. By 
comparing the Bode plots of Figs. 5, 7, 9, and 11 (corresponding to controls using 
both LQR F/B and proportional F/F) with those of Figs. 6, 8, 10, and 12, respec- 
tively (corresponding to controls using LQR F/B only), one can see that adding 
feed-forward substantially improves disturbance rejection at low frequencies, or 
example a comparison of Fig. 5 with Fig. 6 indicates that the optimal control 
method described above can lead to acceleration reductions of greater than four 
orders of magnitude for all frequencies. This reduction is more than two orders 
of magnitude below that afforded by LQR feedback alone at the lower frequen- 
cies, i.e., those most heavily weighted in the performance index. 

The order of the reduction is eventually limited by control cost, of course, 
probably in terms either of actuator-related limitations (such as heat-remova or 
force-generation requirements) or of power limitations (especially in a space- 





FIG. 4. System Block Diagram. 
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FIG. 8. Bode Plots of Transfer Function G4 (Optimal Control with F/B Only. System Stiffness 

Estimate Poor). 
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Concluding Remarks 

This paper has applied an existing method for obtaining an optimal control to 
the microgravity platform isolation problem, for which the disturbances to be 
rejected are low-frequency accelerations. The system was assumed to be repre- 
sentable in the form x = Ax + bu + f, with quadratic cost function J = 
i/ ( 7(x r ^,x + w^u 2 )dt and diagonal weighting matrix W x . The resultant control 
law was found to be simple, stable, robust, and physically realizable. Further it 
was shown to have excellent acceleration- and displacement-attenuation charac- 
teristics, and to be frequency-weighted toward the low end of the acceleration 
spectrum. 

By making an appropriate choice of states, along with the use of frequency 
weighting, the method can be extended to the case for which only relative posi- 
tions and absolute accelerations are available. With this extension one can then 
weight relative displacements in the performance index. [9] Additionally, since 
absolute positions and velocities will then not appear as states, accelerometer 
drift (which affects integrations) will no longer be problematic. Any one of a 
number of noncontacting relative displacement sensors would be suitable for this 
application (e.g., eddy current probes, photoelectric sensors, capacitance probes, 
Hall-effect probes). 

The approach as presented is algebraically intensive, but symbolic manipula- 
tors can be used to ease the algebraic labors. Further, since the method produces 
feedback gains identical to those obtained by the LQR approach to optimality, 
numerical computation to those gains is easily accomplished, even for large sys- 
tems. The feed-forward gains can be found numerically with comparable ease. 
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SUMMARY 


Ho synthesis techniques are developed for a general multiple-input-multiple— output 
(MIMO) system subject to both stochastic and deterministic disturbances. The H 2 
synthesis is extended by incorporation of anticipated disturbance power-spectral— density 
information into the controller— design process, as well as by frequency weightings of 
generalized coordinates and control inputs. The methodology is applied to a simple 
single— input— multiple-output (SIMO) problem, analogous to the type of vibration 
isolation problem anticipated in microgravity research experiments. 


INTRODUCTION 


The vibration environment onboard current and planned manned orbiters requires 
isolation for microgravity science experiments. The disturbance frequencies are sufficiently 
low, and the attenuation requirements sufficiently great, so as to preclude a solely passive 
isolation system (ref. 1). 

Since the disturbances to be attenuated are three-dimensional (ref 2, p.2), the 
isolation actuator must be capable of acting over six degrees of freedom. The requisite 
multiple— degree-of-freedom (MDOF) controller is much more difficult to design than a 
single— degree— of— freedom (SDOF) controller, because the isolation system has many 
inputs (actuator forces) and outputs (measured displacements and accelerations). 
Multiple— input— multiple— output (MIMO) designs can be very susceptible to unmodeled 
cross— coupling between channels of input or output (ref. 3), a problem not encountered in 
SDOF design. The control forces used must therefore be properly coordinated if the 
controller's performance is to be sufficiently insensitive to unmodeled dynamics (i.e., 
robust). The design of a robust MIMO control system requires the iterative use of 
synthesis and analysis tools, the former for controller design and the latter for system 
performance and stability evaluation (ref. 4). 

A particular vibration isolation problem may involve different kinds of undesirable 
outputs, such as excessive absolute accelerations and unacceptable relative displacements. 
Some of these undesired outputs may be more important than others, and the degree of 
undesirability may varv with direction or frequency. For example, rattlespace constraints 
may be highly directional. Or a crystal-growth experiment may be particularly sensitive 
to accelerations at certain frequencies (ref. 2, p. 7) or in certain directions. One of the 
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Croats, then, must be to design a controller capable of minimizing selected plant outputs as 
dictated by these considerations. 

Plant outputs, however, cannot be minimized apart from consideration of the 
associated control costs, because any active control both consumes power and releases heat 
Since both of these costs are of concern in a space environment, the control effort used 
should not be excessive. And at higher frequencies control effort should also be minimized 
in order to limit controller bandwidth for the sake of robustness concerns (ref. 5, p. 218). 

This paper describes a design procedure, known as extended Hj synthesis (ref. 5, p. 
207). for developing active isolation system controllers. A single-input-multiple-output 
design problem is then addressed using the presented procedure. 

BASIC PROBLEM AND SOLUTION 


Problem Statement 

We will use Linear Quadratic Gaussian (LQG) theory to design the MDOF 
controller. This theory has been extensively studied and ussed. LQG is chosen as a 
synthesis procedure since the quadratic performance index relates well to 
root-mean-square statistics and power spectral density. 

When linearized, the differential equations of motion of the plant can be 
representable in state-space form by the first order system of equations 


x = Ax + Bu + Edfd + E s ws 
y = Cx + Du 
z = y + Mn 

where x is the state vector, y is the output vector, z is the measurement vector, ujs the 
control vector, fd is a known or measurable disturbance vector, and vys and n are process— 
and sensor noise respectively. We begin by making a series of reasonable mathematical 
assumptions. Assume that not all states are accessible, so that rank C < dim x. Let the 
initial conditions on the state vector be x (0) = xo; let xo t ws, n, and fd be independent and 
bounded; let xo be Gaussian (ref. 6, p. 272); and let n and Ws be zero— mean white 
Gaussian, with cov[wa(t),w s (r)] = Vi<5(t— r) and cov[n(t), n(r)] = V 3 £(t- t ) (ref. 6, p. 27-). 
Assume that {A,B} and {A,E S VP/ 2 } are stabilizable, w’here Vi = Vi 1 / 2 Vp/ 2 (the 
asterisk here means "conjugate transpose"); and that {C,A} is detectable (ref. 5, p. 226). 
Let Vi and V 3 be positive semidefinite (PSD) and positive definite (PD), respectively. 

We choose a performance index of the form 


C-H 

11 

T 

1 i m i f (W u'l 
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rO 

dt • 
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cr> 
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IT 

U j 


( 2 ) 


where W, is PSD and W 3 is PD (ref. 6; pp. 272, 276). 'T is the expected-value operator, 
needed since the system is excited stochastically by w s . The cost rate functional form 
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(with "1 i m i") is used to allow both for the white noise disturbance w s and for the 
T-oo L 

non— dwindling disturbance f d . 

If Z(t) is defined by Z(t)= (z(r), 0 < r < t}; and if u(t) = a[t.Z(t),f d ] defines the set 
of admissible controls (ref. 6, p. 272), where x is a vector operator that is linear in terms of 
its arguments; the basic problem objective is to find an admissible control function u*(t) 
which minimizes J with respect to the set of admissible control functions u(t). (The 
asterisk here indicates optimality, in the sense defined by Eqn. (2).] 


Problem Decomposition 


The basic problem, as stated in Eqns. (1) and (2), can be decomposed into two 
parallel subproblems, one stochastic and the other deterministic. Suppose that x = Xs + 
Xd, where x s is the portion of the system response due to disturbance w*. and where Xd is 
the portion of the response due to fd- Let v s , yd, z s , z d , Z s , Zd, u s , and Ud be 
correspondingly defined. 


, T 

Then J — 1 i m rp J 
T->oo 1 0 


r , 

Wi WT 


^([x 1 + x'l 

|Xs+ x d} \ 

l Us + UdV 

| \ s -d J 

wi w 3 _ 


dt 


(3a) 


can be reduced to J = Js + Jdi where 


J s = 1 i m rp 
T->oo 1 
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and 
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Jd — 1 i m ixr j 
T-*oo 1 0 
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(3b) 


(3c) 


The problem is now separable into a stochastic— and a deterministic subproblem, each of 
which has an analytical solution. The two subproblems are stated, and their solutions 
presented (without development) below. 


Stochastic Subproblcm and Solution 


Statement : 

Given: x s = A x s + B u s + E s w s (4a) 

v s = C x s + D u s (rank C < dim x s ) (4b) 

z s = y s + Mn (4c) 

{A,B} is stabilizable, {C,A} is detectable 
x s (Q) = x s o is Gaussian with zero mean 
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x . 0i w$? and n are independent and bounded 
such that cov[w s (t), w s (r)] = Vi<5(t— t ) 
and cov[n(t), n(r)] = V 3 £(t-r) 


J s = 1 i m m 

T-oo 


where Vi is PSD and V 3 is PD 

T , fWi W 2 ] 

[x£ lii,] {{£}) 

-v?' w Ms 



dt 


where Wt is PSD and W 3 is PD 
Z s (t) = {z s (r), 0 < t < t}, u s (t) = « s [t,Z s (t)] 


defines the set of admissible controls 


(4d) 

(4e) 

(40 

(4g) 


An admissible control function u s *(t) which minimizes J s with respect 
to the admissible control functions u s (t) 


Solution (See ref. 6, pp. 272-277; and ref. 7, ch. 11): 

u s * (t) = -K x s (t) ( 5a ) 

where x s is an estimate of x s using a Luenberger observer (ref. 7, pp. 

288— 289) having observer gain matrix L 

K = W 3 l (B'P + Wo') ( 5b ) 

P is the unique PD solution to 

PA + A'P — (PB + Wo) W 3 1 (PB + W 3 )' 4- Wi = 0 (5c) 

L = QC' (M V 3 M') _1 ( 5d ) 

Q is the unique PD solution to 

AQ + QA' - QC' (M VaM'^CQ + E s E; = 0 (5e) 

P exists if { A,B } is stabilizable and {C,A} is detectable 
or if the system is asymptotically stable 

Q exists if {A, E s VP/ 2 } is stabilizable and {C,A} is detectable 
or if the system is asymptotically stable 


Statement : 

Given: 


Deterministic Subproblem and Solution 


Xd = A Xd + B Ud + Ed fd 

(6a) 

y d = C xd + D ud (rank C < dim x d ) 

(6b) 

Zd = Id 

(6c) 

{ A,B} is stabilizable, {C,A} is detectable 


ld(0) = IdO 
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Xdo and fd are independent and bounded 


Jd 
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1 i m rpr 

T->oo 1 0 


T / 

< i 


*d lid 


Wi Wo 


m w 3 . 


<g>) 


dt 


where Wi is PSD and W3 is PD 
Zd = {zd( r), 0 < t < t}, Ud(t) = *d[t, Z d (t). fd] 

defines the set of admissible controls 


(6d) 

(be) 


Find: An admissible control function Ud*(t) which minimizes J d with respect to the 

set of admissible control functions Ud(t) 


Solution (refs. 8; 9; and 10, pp. 156—157): 

u d * (t) = -K x d - W3 1 B 1 j^ 1 exp[— X'(t— r)] PE d f d (V) dr 

where K = Wj 1 (B'P + WJ) 

P is the unique PD solution to 

PA + A'P - (PB + W 2 ) Wi' (PB + W 2 )' + W, = 0 

P exists if { A,B } is stabilizable and {C,A} is detectable 
or if the system is asymptotically stable 


(7a) 

(7b) 

(7c) 


Combined Solution to Basic Problem 


When rank C < dim xd, an estimate Xd of Xd must be used in the feedback. If one uses 
an asymptotic (i.e., Luenberger) observer, with gains L chosen to give an optimal solution 
to the stochastic subproblem, he can then combine the stochastic and deterministic 
subproblem solutions so as to use the same observer and regulator. This allows the optimal 
solution (feedback portion) to be realized physically. If such a choice is made, 

u*(t) = u*(t) +u*(t) = -K x(t) -Wi 1 B'j^ exp[— A'(t— r)] PE d f d (r)dr (8a) 
d d u 

where x is an estimate of x using a Luenberger observer 
having observer gain matrix L 

K = Wi 1 (B'P + W 2 ') ( 8b ) 

L = QC'(M V 3 M') _1 (Sc) 

P,Q, and X are as defined previously 

If f s and n are correlated by $ [f s (t), n(r)] = V 2 <^(t— r), then the above solution has the 

modification (ref. 7, pp. 414—417) that 

L = (QC 1 + E s V 2 )(M V3MT 1 . ( 8d ) 

where Q is the unique PD solution to 

AQ + QA' - QC'(M V 3 M') -1 CQ + E s ViE; = 0 (8e) 
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PROBLEM EXTENSIONS 


Frequency Weighting 

Suppose now that it is desired to frequency weight the states x and the control u in 
the cost rate functional, so that the weightings vary with frequency (ref. 11). Let x be 
considered to be the input to a filter Fi(s) of which be is the output, and let *i(s) have a 

state-space representation defined by {Ai,Bi,Ci,Di} [i.e., ^l(s) = Ci(sl— Ai) Bj+ Di]. 


Then 

it = At zi + B t x ( 9a ) 

ix = Ci zi + D| x ( 9 ^) 

expresses J x in terms of x, employing pseudostates zi. Similarly, if u is considered to be the 
input to a filter r 3 ( s) of which «u is the output, and if r 3 (s) has a state-space 
representation defined by {A 2 ,B 2 ,C 2 ,D 2 }, 'u can be expressed in terms of u, employing 


pseudostates z 2 : 


Z 2 = A‘2 Zo + Bo U 
J U = Co Zo + D 2 U 2 


(10a) 

(10b) 


Suppose now that these frequency-weighted states (‘x) and controls Oil) are further 
weighted by constant weighting matrices W t and W 2 , respectively. The resulting state 
equations and performance index are as follows: 


x = *A x + l B u + 'Edfd + *E S w.s 


(11a) 


v = 'C x + D u 


z = y + M n 
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Disturbance Accommodation 

Suppose further that the stochastic disturbance is not w s but f s , where f s is a 
stochastically modeled disturbance with power spectral density 

S(M = S( ./2 (j w )S ( 1,2 *(iui). Defining H,(j<u) by S,'' 2 (ju-) V, 1 ' 2 , one can consider fs to be 
the output of a filter Hi(s) excited by zero-mean white Gaussian noise w s (ref. 12) with 
power Vi (i.e., cov[ws(t), Ws(r)] = Vi <5(t— r). 

In state— space form, 



£ = As £ + w s 

(12a) 


f s = Cs(sl - As) -1 

(12b) 

such that 

H f (s) = C s (sl - A s ) -1 

(12c) 


Incorporating these new pseudostates (£) into the state equations and performance index 


371 



yields 


x = 2 A x + 2 B u + 2 Ed fd + 2 Es w s 


2 J 


v = 2 C x + D u 
z = y + M n 
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1 i m rrr 
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1 (lx' u 1 ] 
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_i\v- tw a 



dt 


(13a) 

(13b) 

(13c) 

(13d) 


where 
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(13e) 

(13Q 

(13g) 

(13h) 

(13j) 

(13k) 
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o 

2\Vo = 0 

I C 2 ' W3D2 1 
0 

2\V 3 = [D 2 'W 3 D 2 ] 

The solution to this problem has been given previously. 


(13m) 


(13n) 


SYNTHESIS MODEL 

The model given at the close of the previous section is the model from which the 
controller is synthesized. The synthesis involves the determination of observer gains L and 
regulator feedback gains K. Preview gains Kp p can also be determined, if desired, to 

approximate the Duhamel integral term of the optimal control One approach to 
determining these preview gains has been presented in reference 9. Further study of the 
determination and use of these gains is needed. 

ANALYSIS MODEL 

Once the controller has been selected, it must be connected to the actual plant and 
the resulting "analysis model" used to evaluate closed-loop-system performance and 
stability. For constant gain matrices K, L, and K FF the open loop transfer function from 

to Upn [— — K X] is 


^ Y = 


— FB— 


2 A- 2 BK-L 2 C 


-K 


0 


(14a) 


where the form 


A 

B' 

c 

D 


indicates C(sl - A) -1 B + D. The closed loop transfer functions, respectively, from Fd and 
F s to X, are 


^ L p d (s) - 


and ^ L f (s) = 


- -s 
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(14b) 


(14c) 


The return ratio matrices (ref. 4) at the Y(s) and U(s) nodes, respectively, for D = 



A 2 -BK-L 2 C 

0 

L ‘ 
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" 2 A— 2 BK— L 2 C 

LC 
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(15a) 


(15b) 


The corresponding return difference matrices and inverse return difference matrices (ref. 4) 
are as follows: 



2 A— 2BK— L 2 C 
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’ 2 A— 2 BK— L 2 C 
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! + lT(s) = I + [K(sl - 2 A + 2 B K + L 2 C) 1 Lj '(C(sl-A) 1 B] 1 
I + I,;‘(s) = I + [C(sl - A) 1 B] ‘(KlsI — 2 A + 2 B K + UCfL] -1 


(15c) 


(ISd) 


(15e) 

(150 


The singular values of these matrices can be used to evaluate system noise and disturbance 
attenuation, stability margins, and sensitivity (ref. 4). Iterative application of the 
synthesis- and analysis models can be used to produce the desired controller. 


EXAMPLE PROBLEM 

Suppose one wishes to develop a controller to isolate a space experiment of mass m 

and position x(t), from a unidirectional acceleration disturbance 3(t). Assume that a wall 
having position d(t) acts on m through an umbilical with stiffness k and damping c. (See 
figure 1). Suppose further that rattlespace constraints require the transmissibility to be 

unity below lO -3 IIz, and that it is desired to attenuate the disturbance by at least two 
orders of magnitude between 0.05 and 10 Hz. Let a linear actuator, applying a force that 
varies with control current i, be connected between the wall and the experiment in parallel 
with the umbilical. 

For this problem, it is desirable at low frequencies to penalize the relative 
displacement of the experiment heavily, so that the experiment "tracks" the wall. At 
intermediate frequencies, however, the absolute acceleration of the experiment should be 
heavily penalized to accomplish the desired disturbance rejection. The state space model, 
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then, should have relative position x— d and absolute acceleration x as states, allowing them 
to be frequency-weighted in the performance index. 


The system equation of motion is 

A ^ * * * 1 * * 

x = — k(x— d)— c(x— d) - oci, where k = c = and * = ^ 


In state-space form, the equations can be written as 
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where 


xi(t) = x(t) d(t) 
x 2 (t) = x(t) - d(t) 

X 3 (s) = (^rrr) s 2 X(s), high, 


S + ' 

so that x 3 (t)~ x(t) for ^ccu^ 

Frequency— weighting the states so that 
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(16a) 


(16b) 

(16c) 

(16d) 


(16e) 


(17a) 


(where wi < a- 1 ’) results in a performance index that penalizes xi more highly at low 
frequencies and" x 3 more highly at intermediate frequencies. If the control is 
frequency— weighted so that 

HJ(s) = u(s) ’ . (17b) 

at higher frequencies the control will be more heavily penalized. This is desirable both for 

the sake of robustness and since x 3 approximates x only at frequencies sufficiently below 
a; h . Finally, let the input acceleration be considered to come from zero-mean uaussian 

white noise filtered through g ^ --. 
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The resultant state equations are as indicated on page 8, where 
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(18b) 
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(13d) 


Ao = — U/'4 

B> = 1 
Co — — uft 

4 

Do = 1^4 
A s = -o>f 
B s = 1 

Cg ttr'f 

D s = 0 


Assume that 


cov 

cov 

and cov 


w s (t), w 2 (r 
ni(t), ni(r) 
n 2 (t), n 2 (r). 


| = 1 ?(t- r) 

= 0.001 S(r- t) 
= 0.001 S(t- t). 


q inrp A, has a zero f 1st! column, 2 A will have a corresponding zero [4th] column. To make 
S free ,e^v-webfitei system PC JA} observable, obtain \ (x-k1 dt as a measured state 
(U. Se K pseudostate, St Zu) and modify >C accordingly. Let the measurement no.se 

associated with Zn be n 3 , such that 

cov [n 3 (t), n 3 (r)] = 0.0001 

Gain matrix W, can be varied to "tune" the optimal control to give the most 

satisfactory results. The transmissiblity between d(t) and x(t) is ; given in figure ; 2. The 
control uses feedback (and observer) gams obtained from system paiameters ana 
weightings as indicated on the figure. Note that the low-frequency transmissib.hty is 
unity, as^ desired, and that for intermediate frequencies the transmissibility rolls off with a 

slope of —1. 
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If a different frequency-weighting of X 3 is used, it is to be anticipated that the 
_ transmissibility curve will change as well. 

For 1 X 3 (s) = ( s+Wl ) 2( s+a , 2 )-2 X 3( s ) (19) 

the resultant selected transmissibility curve is given in figure 3. The low-frequency 
transmissibility again, is unity; but now for the intermediate frequencies the 
transmissibility rolls off with a slope of — 2 , as expected. Adding another pole at uj \ and at 
u> to the Xa(s) frequency weighting would further improve the intermediate— frequency 
roll— off. The present controller, however, meets the design specifications. 

If state frequency— weightings of 1 Xi(s) = g + ^ ■ Xi(s) ( 20 a) 

and 'X 3 (s) = X 3 (s) (20b) 

are used, the results (figure 4) are similar to those given previously in figure 3. Note that 
_ with this latter choice of frequency weighting however, (i.e., without any"rigid body poles"), 

the frequency-weighted system { 2 C, 2 A} is observable, without augmenting the actual plant 
output v as was previously necessary. Consequently this is the preferred control. 


DISCUSSION 


Ho synthesis, as the example problem indicates, provides a highly versatile 
loop-shaping tool. It is especially useful in controller development for SIMO and MIMO 
systems, where classical loop-shaping methods are most lacking. Once the designer has 
expressed the system equations in terms of states for which he has an intuitive feel, and of 
measurable outputs, the design process becomes relatively easy. He frequency weights (i.e., 
filters) the states and control inputs according to his engineering experience and intuition, 
to indicate the relative importance of each as a function of frequency. Then he weights 
these frequency— weighted states and controls relati/e to each other. The II 2 synthesis 
methodology automatically provides him with a set of regulator and observer gains that are 
optimal with respect to the chosen weightings, given a quadratic performance index. 
Known aspects of the input disturbances and sensor noise can be incorporated into the 
design as well. Singular value checks provide the ability to evaluate system robustness. 
With a few iterations, the skillful engineer can complete his design. Excellent computer 
software packages already exist to assist in the task. 

The frequency weighting tells the Ho synthesis machinery how much "cost" to place 
on a state or control input at any frequency, relative to its cost at other frequencies. If, for 
example, absolute acceleration is undesirable only in a particular frequency range, that is 
where it should be most heavily weighted. The subsequent weighting of the 
frequency— weighted states and control inputs tell the synthesis machinery how much cost 
to place on each frequency— weighted state or control relative to the others. 
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Subscripts, Superscripts, and Diacritical Marks 


postsubscript 0 
postsubscript 1 


postsubscript 2 


postsubscript 3 

postsubscript d 
postsubscript f 
postsubscript s 
postsuperscript 1/2 
postsuperscript ' 
postsuperscript —1 
postsuperscript * 
underline _ 
over bar — 


over hat 

overtilde ~ 
presuperscript 1 

presuperscript 2 


VMdf A,EfC,D.z: related to state-frequency weighting 

stated-space description 
With L: return ratio matrix at control node 
With V* nrocess noise covariance . v * 

With w,W,F; state (or pseudostate) weightings, applie 
subsequent to any frequency weighting 
With A,B,C,D,z: related to state-frequency weight g 
stated-space description 
With L: return ratio matrix at output node. 

With V,W: cross-weightings 

With V: measurement noise covariance 

With W , V: control weightings 

Related to deterministic disturbance 

Related to filter for stochastic disturbance 

Related to stochastic disturbance 

Square root or spectral factorization 

Transpose 

Inverse 

Optimum or conjugate transpose 

Vector 

With A: closed loop system dynamic matrix 

With x: augmented with frequency— weighting 

AugmenteTwith frequency-weighting- and 
disturbance-accommodation pseudostates 
Estimated or associated with cross-correlation 

$}£ augmented by 

Related 7o%S by frequency-weighting 
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1. ABSTRACT 


It is well-known that the spacecraft environment deviates from a state of 
zero gravity due to various random as well as repetitive sources. Science 
experiments that require a microgravity environment must therefore be isolated 
from these disturbances. Active control of noncontact magnetic actuators enables 
such isolation. A one — degree— of— freedom test rig has been constructed to 
demonstrate the isolation capability achievable using magnetic actuators. A 
cylindrical mass on noncontacting electromagnetic supports simulates a 
micrograLvity experiment on board an orbiter. Disturbances generated b\ an 
electrodynamic shaker are transmitted to the mass via air dashpots representing 
umbilicals. A compact Lorentz actuator has been designed to provide attenuation ot 
this disturbance. 


2. INTRODUCTION 


Space exploration was initiated for the investigation of space itself, ranging 
from the planetary system to the limits of the universe. Resulting benefits of this 
effort include satellite communications and earth observation and imaging systems. 
The scope of space exploration widened in the early eighties with the development 
of the space shuttle - a system capable of transporting a large cargo to a low earth 
orbit, and recovering the payload or frequently servicing it in space. A parallel 
development was the gradual change in the role of man in space, starting with the 
primarily technical function of a pilot and evolving into a more active involvement 
encompassing interactive work and scientific experimentation in space. 

Space— based laboratories like the Skylab and the Spacelab were flown to utilize the 
"vanishingly low" gravitational forces available for extended periods of time, the 
results, however, were mixed at best, and disappointing in certain cases. This can 
be explained in part by the fact that the environment aboard the spacecrafts 
deviates considerably from the ideal of zero gravity due to disturbances produced by 
machinery and people on board, thruster fire, and other factors. 

The incentives for performing science experiments in space include the 
investigation of phenomena that are influenced by gravity on earth, the 
development of novel materials and the improvement of processes like crystal 


This work was supported in part by the NASA Lewis Research Center and the 
Center for Innovative Technology of the Commonwealth of Virginia. 


PREGEDiHG 


PAGE BLAUK * 


JOT FILMED 


437 



growth [1], In theory, a freely orbiting spacecraft offers a state of zero gravity to 
objects inside it, since the gravitational force is balanced by the centrifugal force [2]. 
However, in practice, there are various residual forces that'disturb the environment. 

Attempts to estimate these residual forces have been made in the past few 
years [3-10]. The orbital microgravity environment can be divided into three 
classes, as detailed in Table 1. Quasi-steady accelerations are generated by three 
sources — gravity gradient, aerodynamic drag and rotational acceleration.' Any 
point of an orbiting structure that is at a distance from the structure's center of 
mass experiences a gravitational field that is different from that at the spacecraft 
center of mass. Aerodynamic drag due to the earth's atmosphere represents the 
absolute lower limit of the achievable background microgravity level, if the effect of 
light pressure is neglected. Finally, in order to keep the same vertical orientation on 
the space station with respect to the earth, the station must maintain a constant 
pitch rate about its center of mass. This creates a centripetal force that results in a 
rotational acceleration. 


Orbital thruster fire and the steady— state operation of machinery like fans 
and pumps on board a spacecraft are among the sources of periodic accelerations, 
w'hich occur at known frequencies. Impulsive disturbances like crew push— off and 
the start— up and shut— down of machinery create non-periodic accelerations. The 
irregular, unpredictable nature of these accelerations complicates attempts at 
isolation. 

Theoretical acceleration requirements for various processes and experimental 
conditions have been investigated [3,7.11], The common feature of curves depicting 
the frequency— dependent requirements is that, for a given process, the acceleration 
threshold is lowest from steady— state to about (0.01 — 0.1) Hz. depending upon the 
type of experiment. I he acceptable acceleration then increases linearly with 
increasing frequency, up to (1 — 10) Hz. Subsequently, it increases as the square of 
the frequency. The acceleration tolerance also typically scales inversely with the 
volume that characterizes the process. The slopes and breakpoints result from 
fundamental aspects of a process, and the shape of the curve can be considered to be 
characteristic of a family of experiments. The acceleration level thresholds range 
from an extreme level of (10 1 — 10'*) g 0 for some material science and fluid science 
experiments, to only 10 3 g 0 for the majoritv of biologv and biotechnology 
experiments. v ' 

A comparison of the micrdgravity requirements with the actual environment 
available on the spacecraft indicates the need for vibration isolation. Moreover, the 
frequency range of interest spans several decades, thus requiring the use of multiple 
strategies for isolation. 


For the high frequency range, passive isolators can serve adequately. Since 
these are relatively simple and cheap, they can be placed at each interface' between a 
disturbance source and the space station. It should be noted that the sensitivity of 
various categories of experiments to high frequency disturbances is also 
comparatively low. 

In the quasi— static frequency range, the extremely low stiffness and lar^e 
motion required make attempts at isolation very difficult. Rattlespace constraints 
prohibit the occurrence of such large relative motions between the pavload and the 
spacecraft. This imposes a fundamental limitation upon vibration isolation. 
Consequently, efforts at minimizing the input disturbances, like reducing the surface 
area presented to the atmosphere so as to reduce atmospheric drag and locating 


pav loads as close to the spacecraft's center 

Owens and Jones have also suggested the possibility of cancelin e such disturbances 
by continuous thruster control of the whole spacecraft [1 ]. 

At intermediate frequencies — approximately between 0.01 Hz and 1 Hz — 
nA cci V e isolation scheme can be effective due to the displacements and isolation 
fevelf required “only active vibration isolation at the t m erface 

allows the synthesis of the desired isolator properties and the adjustment o 
properties using a control loop. 

The actuator used to implement an active controi scheme in the 
intermediate frequency range should ideally be noncontact ^ ba been 

acoustic, electrostatic and electromagnetic ( nl miterials in a low gravity 
considered in the contexl ; o fqj es ar e1imited to small objects. Lorentz 

currents also occurs. 

Noncontact magnetic actuators, utilizing electromagnets or pe rmanent 
Sid to the isolation of a variety of science experiments. An active magnetic 

SotmXeTh^^^ 

s s^tssssirtsrer^ w-c-a. 

purpose^Vh*^^^ 

™ ions not exceeding about a hundred miles Since the forc« cated ^ «e oftta 
order of a few pounds at most, such actuators can be quite compact, tne 8 

primarily determined by the stroke required. 

been de^crlbedTml^^he'au^ors^ornp^^'v^rimt^ 

to shaker and accelerometer limitations, the . . simi | ar isolation system 

ssss 

Magnetic Isolation System (FEAMIS) Izzj. 

An interesting dual-mode approach to vibration isolation of large pa>, l °ads 



i 4. ^ motion &ctUcLtor w&s controlled &s & 

systems. In such a tandem system a coarse -t displaccment (or the magnetic 

followup actuator, always a tte ™P l S g t5c actuator functioned as a 'Tine 

s— » *»« 

* lhc fi " e actuator - 

The operation of microgravit^scienc a h'S'sVa^and ^ US<> 
of an umbilical. An examp le is P ’ tQ re ; ect disturbances caused by the 

carrying a coolant. Acce er theoretically investigated [24]. The umbilical 

compliance of the umb heal has damping "The microgravity quality 

was assumed to have stiffness, but not P s ‘ expec ted Acceleration control 
deteriorated with increasing urn 1 1 ,“ when compared to position-only control, but 
improved disturbance re Je c i ? f or " m 0 f a m0 re complicated control system. An 
there was a price to be paid n thelormot a mo q{ the umbiUcal 1S 

active umbilical control strategy . ^ a j s0 analyzed. It was found to 

in performance to the acceleration control 

loon technique. 


3. experimental rig 

An experimental rig to demonstrate 
section, will be used to implement the isolation scheme. 

An electrodynamic shak^ with a ^on^^e^-^^P^^stroke ^ Qn 

represents the space platform. Jhe shaLe u ^ generate sinusoidal, 

concrete block resting on the thus simulating the 

dtsturtancestj^toHyT 1 ^'*^ a'space station that require active isolauon. 

The umbilicals connecting a semnee ex^riment^o^h^pa P (lash t5 w j t h 
expected to be flexible hoses and wires -These will > has been 

adjustable stiffness and damping cw 11^^^ mass-spring-damper 

evaluated at NASA Lewis in a The test indicated the 

system in a fixed-fixed m ®^ in | /dam ® in g mechanism in these air dashpots. The 
possibility of a nonlinear stiffness/ damp g mec kinds 0 f umbilicals may be 

vibration isolation rig has been < J® Ser. This is important, since 

employed, including actual hoses used for « J> d ation t0 microgravity levels 

" lhe experimenl and the 

space platform. 

The long stroke Lorentz actuator Jn tv^^nc^exp^ 

shaker armature to lhe it^^cvUncler 0 weighing 8 75 p6unds, which is a typical 
space. This mass is a solid st .. , er - |g horizontally suspended in space b> 

weight for such an experime • - , e [ ec tromagnetic support system, 

the magnetic forces generated ^by a noncontart elect o „ of lw0 eighl -p 0 le 

Similar to radial magnetic bearings the supp« r ( the cylindrical mass, 

structures, mounted on a concrete base at he two ^en ^ ^ ^ ^ 

™ block on which the shaker is 
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mounted. Eddy current probes sense the radial position of the cvlinripr 
complete the feedback loop supplying current to the electromagnets. When the 
e ectromagnetic support system is turned off, the cylinder rests°on a pair of 
touchdown pedestals made of delrin. p 


™ acce leration of the cylinder will be sensed off a sensory nlate 

resSo^ofn^P^’ USing , a v , er >' low frequency accelerometer with a maximum 
° n 1 , Pr l 0vlsl0n has been made for the use of other types of 

accp|pmmi terS ’ and i the i, S u nS i n ? of other states of the system, like position. The 

current rpniUrpH° na in W1 i de fed t0 a fee^^ack control circuit that determines the 
current required in the electromagnetic actuator to isolate the cylinder from the 
disturbance generated by the shaker. A control strategy fofsuch an isSion 
s\stem with multiple degrees-of-freedom is discussed m [26], 

The background vibration levels on the concrete base on whirh ti ia ^ r a 
.s mounted have been measured over twenty-four-hour S.t lS GF 
oriental and the vertical directions. These vibrations are of the order of milli 
the quietest period occurring from late in the night to earlv in the morninl § 

Figures* ^and^ shn^tT’f y ‘ eld the highest de ? ree of reproducibility in our results. 

horizontal and vertical dir^ns^ThTS^ acceleration ‘ he 

components corresponding to natural frequencies of the mounting plate The 

nwr Z mode. VlbratIOn ^ significant content at 45 Hz - The authors believe this is a 


4. THE LONG-STROKE LORENTZ ACTUATOR 


in tho r^K COr ^ paCt 1 ° n Sr stroke Lorentz Actuator has been designed, built and tested 
m the laboratory. An intermediate version of the design was presented at the 

T an a ip ° P 0 ° ‘^ er °sP ace Applications of Magnetic Suspension Technology at NASA 

0 h^ame S S? r bm h The fi , nal d “ ^ ™ 

accomplished through modification of the flux distribution. 

used i si h o w nO n 1 F i rO \ ^ orentz Actuator, along with the terminology 

core ’ a lOnn. 1 F gUre 5 ; The current car rymg coil moves in and out along the 
tho n r i • ? R ermanen ^ ma g n et in the shell maintains a constant magnetic flux in 

tdtSndeS limifaf ^ ^ irreSpectite of the ^emlfthe c 0 , 

IS "JuhSl curre’nt ps) ent2 """ therefore ' “ n “■»* 

* re Quirements for the laboratory prototype were fixed at a total sfmkp nf 

umbllicaUal? daKp 5? s Capabiiitv t0 > solale * of TMte coVnec^ b v an 
i- • v, na.hpot) to a source generating verv low frequency vibrations Force 

meaxity with position and with current were also'required^ Moreover in view of 
iah, U |S, te r' ° f depio - vment such a dev?ce had ,o Compact T,d 

were also important^ ^ C0nsumptl0n and low heat generation during operation 

r ^ computer program was written to implement a simple design algorithm for 
a Lorentz Actuator. The steps of this algorithm are presented in Figure 6. 

Using a permanent magnet material with a very high maximum ener<*v 
product of 35 MGOe (mega-Gaus s -Oersted) (29) resultedTn Tde^hat juiced 



r o on nntpr diameter. The magnet manufacturer, however, could 
a ring magnet of 3- ■ . j -f |ts outer diameter were less than 2 in.; 

make such a magne^one p.ece^ ^ ^ ^ h&ve reqllire d the costly 

assembly ofmultiple segments, with an escalation of costs. 

The possibility of designing a Lorentz ^ c ‘ u n ^ constrained 

requirements, witf ^^““fore^p ored' The significant parameter in this context 

in order 

by the design program sug Q este thumb of using a gap ratio of 5:1 or 

with the additional constraint, « which the computer 

more were drastically violated. PP : the°failure of such a design. It 

program did not take into accou , neC essarily be the case if the core of the 

was then hypothesized that this w o ti T h e permeability of a saturated 

actuator were saturated during normal . o^rrtion. uiepe^ ^ ,^ kage that 

wouhThav^oc^ a low gap ratio (less than 2d) and an nnsaturated 

core, would be prevented. 

A number of designs ‘TThffinil^ 
finite element analysis package. The finite Umitati £ ns . The mesh 

achieve as much accuracy as possible, f . air elements used earlier, were 

consists predominantly of quad elements. about 

found to cause severe restrictions on mesh finen^ A ^mesh with 0^ ^ ^ 

100 elements could be used. An air thickness o determine d t0 be as accurate 
asymmetric model was specified ms ea£ « “ while a 

as having infinite air elemeM ■ all brees d^tor a^mode cnc0UIUering 

the finest mesh allowed by the configuration of 

our 386-based personal computer was used for the ana } 

Position linearity was improved, relative to the |f ou ter 

t^vare^ W °' The u* of » iSr^StlfiS.Sd^ 

the design of a compact and hgtoe.J ht MtuMO ^ bc f made = in onc piece , lhus 

have required a large diameter magnet t , f caused bv the increase in 

increasing costs. The decrease in flux and the ^“{^^^uction in the 

diameter is only 1.95 in.. 
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The salient features of the final design of the compact Lorentz Actuator are 
described below: 

• Long Stroke — The requirement of two inches of total stroke is satisfied. 

• Position Linearity — Over the whole two inches of stroke, the actuator exhibits a 
high degree of linearity. For a constant coil current, this means that the actuator 
force is the same irrespective of the axial position of the coil, within the stroke 
bounds. Figures S and 9 depict this relationship for positive and negative coil 
currents respectively. Note that flux leakage has been reduced to almost zero over 
the shell— to— core gap to achieve this. The maximum flux density across the 
shell-to— core gap is only about 7% of the maximum flux density across the 

pole— face gap. 

• Current Linearity — This requires that the average flux density in the effective 
air gap remain constant with variations in the coil current between the upper and 
the lower limits. This is indeed the case, resulting in force vs. current linearity, 
Figure 10. 

• Force — A maximum force of 1.25 lbs is produced by this actuator, which is 
sufficient for our needs. This peak force requires a coil current of 2.5 A. 

• Weight — At 2.28 lb., this actuator is only a tenth of a pound heavier than the 
previous design. 

• Current Density — A value of 1000 A/sq. in. in continuous use ensures cool 
operation. For peak loads, a fivefold increase in current density is possible. 

• Materials — The magnet is made of neodymium iron boron, which has a very 
high maximum energy density product of 35 MGOe. Selection of such a material 
helped make the design compact. The high permeability circuit material is a 48% 
nickel— iron alloy that saturates at 15 kG. The B— H curve for this material, 
provided by the manufacturer, was input to MAGGIE as a table of a large number 
of points on the curve. This was necessary because a nonlinear material 
characteristic was being modeled. 

The design specifications of the Lorentz Actuator are detailed in Table 2. 

This actuator was built and tested in our laboratory. Figure 11 compares the 
measured magnetic flux density in the radial direction along the shell— to— core and 
pole— face— to— core gaps with the values predicted by finite element analysis, for no 
current in the coil. The measured peak value is lower, but is spread over a wider 
axial distance. There is good agreement, especially over most of the shell— to— core 
gap, where near— zero values of flux density are crucial to achieve force versus 
position linearity. The actual actuator force is plotted against, position for a number 
of values of coil current in Figure 12. The measured values of force are greater, in 
each case, than the predicted values since most of the small amount of leakage flux 
across the shell— to— core gap was neglected in calculating the predicted forces. 
Moreover, since the coil does see slightly greater total flux as it moves into the 
actuator, because of the small amounts of leakage, the forces measured increase 
somewhat with such motion. However, for low values of current and for coil 
positions that do not place it very near the closed end of the actuator, the actual 
forces deviate by less than 10% from the predicted values. 


5. CONCLUSION 


The riff designed to demonstrate vibration isolation to rnicrogravity levels in 
nn^imension has been built and assembled. Measurements of the background 
acceleration levels have also been made, and the quietest period for °P erat ’° a ^ , 
heen determined A compact, long stroke Lorentz actuator has also been de. e ’ 
S and tested' Its performance has been shown to match that predicted by finite 
element anaksis very well. Microgravity isolation experiments wtll be conducted m 
the very near future. 


444 



6. REFERENCES 


^Siiatsass. sgsig^ga as" s — 

„ • • "i Arrplpration Characterization of Space 

Teledyne Brown Engineering. Oct. 1935. 

Station Environment," Report No. SP 8 o 8 

Chase T L "Repon on Micro-g Measurements for Space Shuttle 
Experiments," NASA lewis Research Center, Dec. 1985. 

m • • „ AKctrartq WorkshoD on Measurement and 

JSrS»“ Sr.;^ Board the Space 

Station. Guntersville, Alabama, Aug. 11-14, 

TT , ur Tiler and V Merbold. "Analysis of Microgravity 
Hamacher, H., R. Jug anci l . • p ■ Svmposium on 

Bordeaux, Dec. 2-5. 

1986. 

Booz, Allen & Hamilton, M °verview^f Space Sution office 

Requirements," report presented to Dr. J. D. Bartoe, v. 
of Space Station. July 14, 1989. 

Dynamics), Dec. 1985. 

, ^ , n A Winrpr "The Effects of r— Jitter and Surface 

“induced Convection on Float Zones.^ 28th Aerospace Sciences 
Meeting, AIAA, Reno, Nevada. Jan. 8 11, 199U- 

\ treant r w "Two— Bodv Dynamic simulation of Space Station 

S hrill 1 Transient Effects," Memo No. 

A 96 —J 749 —STN—M—RW A— 900006. Jan. 10, 1990. 

Owen, R.C., and D.I. 

kh^l of Etecttnic Engineering Science, Bangor. Gwynedd. Sep. 1988. 

V „„„ R i and D D Elleman, "Containerless Processing Technology." 

in space , Ed. B. Fenerbacher. et ah. 

Springer— Verlag, Berlin. 1986. 

Hendricks, C.D., "Levitation. Coating, andjransport of P F ^,;^ e pnlJ2 L 
Materials," Materials Processing in the Redu j j < Scie n^Tublishmg 
Space , Ed. G.E. Rindone. et al„ vol. 9, pp- otf, tisevin o 

Co., Amsterdam. 1982. 

■ ^Development of Eh^tro^aBc^Positiopr fo^^^ 

Rev. Sci. Instrum., vol. 06 , no. 2, p. 30i, 1985. 



27. 


Banerjee. B.B., P.E. Allaire and C.R. Knospe, "Vibration Isolation of 
Science Experiments in Space — Design of a Laboratory Test Setup," 
Workshop on Aerospace Applications of Magnetic Suspension Technology, 
NASA CP-10066 Part 2, Paper 26, March 1991. 

28. Carlson, A.B., D.G. Gisser and F.K. Manasse, "Magnetics and 
Electromechanics," Chapter 17, Electrical Engineering: Concepts and 
Applications , Addison-Weslev Publishing Company, Reading Massachusetts, 
1989. 

29. McCaig, M. and A.G. Clegg, Permanent Magnets in Theory and Practice . 
John Wiley L Sons, New York. 1987. 


44 


15 Frost, R.T. and C.W. Chang, "Theory and Applications of Electromagnetic 
Levitation," Materials Processing in the Reduced Gravity Environment of 
Space, Ed. G.E. Rindone, et ah, vol. 9, pp. 71, Elsevier Science Publishing 
Co., Amsterdam, 1983. 

16. Grodsinskv, C.M. and G.V. Brown, "Nonintrusive Inertial Vibration 
Isolation Technology for Microgravitv Space Experiments," submitted for 
publication to Journal of Spacecraft and Rockets, AIAA, Jan. 1990. 

17. Allaire. P.E., M.A. Scott and B.B. Banerjee, "Magnetic Actuators for 
Microgravity Space Isolation," Workshop on Vibration Isolation Technology 
for Microgravity Science Experiments," NASA Lewis, Cleveland, Ohip, Sep. 
28-29, 19S8. 

IS. Banerjee. B.B., "Analysis and Design of Magnetic Thrust Bearings." M.S. 
Thesis, Univ. of Virginia, Charlottesville, May 1988. 

19. Allaire, P.E., A. Mikula. B.B. Banerjee, D.W. Lewis and J. Imlach, "Design 
and Test of a Magnetic Thrust Bearing," Journal of the Franklin Institute, 
vol. 326, no. 6. pp. S3 1-847, 19S9. 

20. Allaire, P.E., J. Imlach. J.P. McDonald, R.R. Humphris, D.W. Lewis, B.B. 
Banerjee, B.J. Blair, J. Claydon and R.D. Flack, "Design, Construction and 
Test of Magnetic Bearings in an Industrial Canned Motor Pump," Proc. of 
the Sixth International Pump Users Symposium, Texas A&M Univ., College 
Station, Nov. 1988. 

21. Havenhill, D.D. and K.D. Krai, "Payload Isolation Using Magnetic 
Suspension," A AS 85-014, Annual A AS Guidance and Control Conference, 
Keystone, Colorado. Feb. 2—6, 1985. 

22. Allen. T.S., D.D. Havenhill and K.D. Krai, "FEAMIS: A Magnetically 
Suspended Isolation System for Space— Based Materials Processing," AAS 
86-017, Annual AAS Guidance and Control Conference, Keystone, Colorado, 
Feb. 1—5, 1986. 

23. Hamilton, B.J., J.H. Andrus and D.R. Carter, "Pointing Mount with Active 
Vibration Isolation for Large Payloads," AAS 87—033, Annual AAS 
Guidance and Control Conference, Keystone, Colorado, Jan. 31 — Feb. 4, 
1987. 

24. Jones, D.I., A.R. Owens, R.G. Owen and G. Roberts, "Microgravity 
Isolation Mount: Design Report," Technical Note No. BTN-009, University 
College of North Wales, School of Electronic Engineering Science, Bangor, 
Gwynedd, Sep. 19S9. 

25. Sutliff, T., "Vibration Isolation: Airpot (Dashpot) Damping 
Evaluation," PIR No. 89-9. NASA, Lewis, Cleveland, Ohio, 19S9. 

26. Hampton, R.D. and C.R. Knospe, "Extended H.-, Synthesis for 

Multiple— Degree— of-Freedom Controllers," International Symposium on 
Magnetic Suspension Technology, NASA CP— 3152, 1992. 


446 



1 e- 7 g 
1 e -8 g 
1e-7 g 

1 e -2 g 

le-3 g 
le-4 g 

le-4 g 
le-4 g 


(0 ,0 1.-3) H, Aerodynamic Drag 

(0 to 1e-3) Hz 

(0 to 1e-3) Hz 

pcpinniC ACC ELERATIONS 
9 Hz 


Light Pressure 
Gravity Gradient 


Thruster Fire 
(Orbital) 


(5 to 20) Hz 
17 Hz 


Crew Motion 
Ku Band Antenna 


^-pfriODIC _ A££ELEBAIiQNS 
1 Hz 


Thruster Fire 
(Attitudinal) 


1 Hz 


Crew Push-Off 


Table 2: Design Speci fications for m^Lorentz^ctuator 


I QRENTZ_ AdJiMQB_ 

Total length 

Magnet outer diameter 

Magnet inner diameter 

Magnet length 

Shell outer diameter 

Pole-piece thickness 

Core diameter 

Air gap 

Shell-to-core gap 
Gap ratio 

Coil length 
Coil wire diameter 
Number of turns 
Number of layers 
Maximum coil current 

Air gap flux density 
Max. force generated 
Actuator wt. (no coil) 


final DESIGN 
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1.95 
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600 turns 
4 layers 
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Fig. 3: Background Vibration in the Laboratory -- Horizontal 



Fig. 4: Background Vibration in the Laboratory -- Vertical 
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Lorentz Actuator 



Axis of Symmetry 


Closed 

End 


Fig. 5: Schematic of a Lorentz Actuator 


LORENTZ ACTUATOR . DESIGN EQUATIONS 

1. Assume permanent magnet operating point for maximum 

energy product : (-HI, B1). 

2. Compute magnet flux, fm ■ B1 * Am. 

3. Compute circuit flux, fc ■ HI * Lm / R, where R is 

the circuit reluctance. 

4. Compare fm and fc. 

5. Adjust operating point until fm ■ fc ■ f, the actual 

operating point. (When saturated, f ■ saturation 
flux in saturated segment of circuit.) 

6. Calculate air gap flux density, Bg - f / Ag. 

7. Compute force capability, F ■ i * I * Bg, where i is 

the actuator current and I is the total length of 
coil wire in the air gap. 

8. Change actuator geometry or circuit / magnet 

material until desired force level is achieved. 


Fig. 6: A Simple Algorithm for Designing a Lorentz Actuator 



Fig. 7: The Compact, Long-Stroke Lorentz Actuator 
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Fig. 8: Force vs. Position (for Positive Coil Currents) 
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Fig. 9: Force vs. Position (for Negative Coil Currents) 
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Fig. 10: Force vs. Current (for Three Different Coil 
Positions) 
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Abstract — Active vibration isolation systems contemplated for microgravity space experiments may be 
designed to reach given performance requirements in a variety of ways. An analogy to passive isolation 
systems proves to be illustrative but lacks the flexibility as a design tool of a control systems approach 
and may lead to poor designs. For example, it is shown that a focus on equivalent stiffness in isolation 
system design leads to a controller that sacrifices robustness for performance. Control theory as applied 
to vibration isolation is reviewed and passive analogies are discussed. The loop shaping trade-off is 
introduced and used to design a single degree of freedom feedback controller. An algebraic control design 
methodology is contrasted to loop shaping and critiqued. Multi-axis vibration isolation and the problems 
of decoupled single loop control are introduced through a two degree of freedom example problem. It 
is shown that center of mass uncertainty may result in instability when decoupled single loop control is 
used. This results from the ill conditioned nature of the feedback control design. The use of the Linear 
Quadratic Regulator synthesis procedure for vibration isolation controller design is discussed. 


1. INTRODUCTION 

Active vibration isolation for microgravity space 
experiments has generated much interest lately. A 
variety of disturbances on-board manned space 
orbiters contaminates the desired microgravity en- 
vironment. These accelerations cover a frequency 
band from d.c. to 100 Hz. Low frequency ( < 10“ 3 Hz) 
sources include drag, solar pressure oscillations, tidal 
effects, and gravity gradient forces. At the higher 
frequencies, manned activity, thruster firing, and 
orbiter systems contribute most significantly. A com- 
prehensive treatment of the orbiter accleration en- 
vironment is presented in [ 1 ] from which Fig. 1, a 
characterization of the environment, is taken. 

The need for the active isolation of materials 
processing and fluid science experiments in the fre- 
quency range 0.01-10 Hz has been demonstrated by 
Jones et ai [1-3]. Above this range passive isolation 
systems could be used. Below 0.0 1 Hz the rattlespace 
available for the experiment is not large enough to 
accommodate the relative motion. Therefore, these 
accelerations must be passed by the isolation system 
to the experiment. 

Active isolation systems for microgravity and 
pointing applications have been designed and con- 
structed by many investigators [3-5]. These systems 
generally use conventional PID control of a non- 
contacting actuator, either Lorentz or electromag- 
netic, to achieve 1ow f frequency disturbance attentua- 
tion. While an actual microgravity experiment may 
require umbilicals for cooling and power (at this 
point, it is not clear whether these functions can be 
performed otherwise as described in [4]) the isolation 
systems designed and tested so far preclude an umbil- 


ical from consideration. These systems achieve their 
performance by the very low stiffness made possible 
by low gain feedback of the relative position of the 
experiment to the mounting surface. Without an 
umbilical this stiffness may be set by the designer at 
will. However, when an umbilical is present, the 
umbilical stiffness presents a lower bound on achiev- 
able stiffness unless the feedback loop is used to 
introduce a negative stiffness. In this paper, the issues 
of control system design for the generic (i.e. with 
umbilical) microgravity experiment will be con- 
sidered. 

Previous research in the area of active microgravity 
vibration isolation has established the importance of 
the umbilical in control system design. Jones et ai [6] 
present a good preliminary examination of the single- 
degree-of-freedom control issues for intrusive and 
non-intrusive isolation systems. Grodsinsky [7] 
examined the use of acceleration and velocity feed- 
back. Many of the issues these researchers have 
discussed are revisited here from a control theory 
perspective. Analysis of the six-degree-of-freedom 
problem in the literature has been restricted to one- 
loop-at-a-time design. Generally the effects of cross 
coupling between the various degrees of freedom 
have been ignored. Owens and Jones [2] have investi- 
gated the effect of cross coupling due to center of 
mass displacement for a single loop based controller. 
Their work examines this important problem for the 
non-intrusive experiment platform case where rela- 
tive position feedback is sufficient. The authors con- 
cluded that satisfactory performance can be achieved 
if the control loops are designed for the decoupled 
degrees of freedom and not autonomously for each 
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Fig. 1. The microgravity environment (from [1]). 


local position. It should be noted that high gains are 
not required to achieve isolation for the umbilical- 
free case. An example is presented in this paper which 
shows that decoupled single loop design may not be 
sufficient for the generic isolation problem. 

Any microgravity isolation system design should 
meet the following specifications for translational 
axes: 

(1) Unity transmissibility from d.c. to 0.001 Hz so 
as to prevent the experiment from impacting its 
enclosure's walls. 

(2) At least 40 dB attenuation above 0.1 Hz [3]. 

(3) Both stability and performance robustness 
with respect to changes in umbilical/experiment prop- 
erties, non-collocation or misalignment of sensors 
and actuators, center of mass uncertainties, and 
unmodeled cross coupling between the degrees of 
freedom. 

Robustness refers to the ability of the control 
system to perform satisfactorily when the true plant 
varies from the nominal plant. Performance require- 
ments of the type (2) for rotational degrees of free- 
dom have not yet been specified to the knowledge of 
the authors. 

In this paper we shall examine the control system 
issues associated with active microgravity vibration 
isolation. The purpose here is not to develop new 
control theory but to apply existing concepts to the 
problem. We hope that this paper will serve a tutorial 
function for vibration engineers involved with the 
microgravity problem. The thesis of this paper is that 
control system design, not passive isolator design 
familiar to vibration engineers, is the proper tool for 
analysis and synthesis. First, the control theory re- 
quired for the examination is reviewed in Section 2. 
Section 3 reviews passive isolation and applies it as an 
analogy to control system design. In Section 4 classi- 


t The Nomenclature is given in the Appendix at the end of 
this paper. 


cal loop shaping is applied to the isolation problem 
and a controller is designed. A discussion of the result 
and a passive system analogy follow. An example 
multi-degree-of-frcedom system is explored in Section 
5 and system robustness is examined. Section 6 
concludes with an examination of the Linear 
Quadratic Regulator for the isolation problem. 


2. CONTROL THEORY PRELIMINARIES! 

We examine here the prerequisite control theory 
for the examination to follow. While the actual 
isolation problem is multi-dimensional, a single- 
degree-of-freedom example will be examined first. 

The one-degree-of-freedom microgravity vibration 
isolation problem, depicted in Fig. 2, consists of an 
experiment of mass m connected by an umbilical and 
an actuator to a w r all of the experiment enclosure. The 
umbilical is modeled here as a linear element with 
stiffness k and damping c. The wall’s motion (dis- 
placement y) is transferred through the umbilical to 
the experiment resulting in its motion (displacement 
x). Direct disturbances may also act on the exper- 
iment due to the experiment’s processes (e.g. motors, 
valves, shutters). While it may seem that there is no 
need to distinguish between umbilical and direct 
disturbances, they are indeed different. The distinc- 
tion lies in the fact that the actuator influences 
through the experiment's motion the force transmit- 
ted through the umbilical; direct disturbance forces, 
however, are independent of actuator force. This 
distinction carries through to both passive isolator 
performance and control system design. 

The equation of motion for the experiment is 


nix + cx + kx = cy + ky + d -hf (I) 

where d is the direct disturbance force and / is the 
actuator force. We assume here that the spacecraft 
wall is of sufficient impedance so as to not be affected 


EXPERIMENT ENCLOSURE 



Fig. 2. The onc-degree-of-frcedom microgravuy vibration 
isolation problem. 
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by the actuator force. Under Laplace transformation 
eqn (I) yields 

m-\ — r, +k r l^CO 

[_w 2 + « + /tj 

+ [^T77Tl} Z)(j) + ^)J < 2a > 

or 

*U-\ -Z +k r lr(i) 

+ [^ +C5 + fc ) g(j) + ^«- < 2t » 

This is illustrated in the block diagrams for the 
isolation system [Fig. 3(a) and (b)J. Here, H(s) is the 
feedback transfer function, T(s) is the feedforward 
transfer function, and i»i(r), c 2 (j) are measurement 
noises. The actuator force is therefore a linear func- 
tion of the wall and the experiment motion. The 
subscripts p and a throughout this paper refer to 
whether the model used is in position or acceleration 
form. 

If the umbilical properties are known explicitly and 
measurement noise is sufficiently small, then trans- 
mitted disturbances can be rejected with only feedfor- 
ward control. Note, however, that direct disturbances 
can only be attenuated through feedback. As always, 
the primary purpose of feedback here is to account 
for uncertainties, either in the disturbance or in the 
plant model. 

The price paid for this property of feedback is the 
requirement that the feedback be stabilizing over the 
range of uncertainties in the nominal plant , the plant 
model assumed for design. The nominal stability of 
the closed loop system may be checked by a variety 
of methods, the most popular for single-input— single- 
output (S1SO) systems being the Nyquist and Bode 



plots. Implicit in these methods are measures of 
system robustness. The Nyquist stability criterion can 
be generalized to multi-input-multi-output (MIMO) 
systems, however the robustness measures do not 
carry over as straightforwardly. 

Both Fig. 3(a) and (b) can be generically expressed 
in the form of Fig. 4 where G(s ) is the plant, P(s ) is 
umbilical’s pre-compensation of the wall disturbance 
and Z5(s) is the equivalent disturbance to the system. 
Figure 4 has been presented in unity feedback form 
so as to introduce the concept of loop shaping and the 
trade-offs inherent in control system design. Denote 
the transfer functions between D(s ) and A'(r), the 
sensitivity function, as 


5(s) = 




D(s) 1 + GH 


(3) 


and between v 2 (s) and x (,v ), the complementary sensi- 
tivity function, as 


C(i) = 


yu) 

l' 2 (s) 


GH 

1 Vgh' 


(4) 


(a) 


D 



(b) 


D 



Fig 3. Displacement (a) and acceleration (b) isolation system block diagrams. 
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Note that 5 (j) + C(j)= 1. Therefore, a feedback 
controller designed to attenuate external disturbances 
at a particular frequency 

IS(M)I <? 1.0 \GH(ju) 0 )\i> 1.0 

cannot attenuate the measurement noise signal at that 
frequency 

\C(ja> 0 )\ % 1.0 , 

Likewise, a controller designed to reject a certain 
frequency measurement noise, |C(ycu 0 )| <? 1.0, must 
pass the external disturbance at this frequency, 
l-HM)! ft 1.0. Classical design of control systems 
usually involves separating (if possible) the frequency 
spectrum into regions where input disturbances 
(measurement noise here) and output disturbances 
(external disturbance here) predominate. The meth- 
odology, known as loop shaping, consists of choosing 
ff(s) so that GH is large and therefore S(s ) is small 
at frequencies where output disturbances are domi- 
nant, and choosing H(s) so that GH is small and 
therefore C(s) is small at frequencies where input 
isturbances are dominant. This would be a relatively 
simple task if the designer only needed to be con- 
cerned with the magnitude of GH. However, stability 
of the feedback system requires that the argument of 
GH at crossover, where \GH(ja> 0 )\ = |,0, be 
> - 180°. That is, the system must have some phase 
margin. Since the phase of a transfer function is tied 
to the derivative of its magnitude (in dB) with respect 
to frequency, as was shown by Bode [8], the loop 
shaping’s results are fundamentally limited by the 
difference in frequency between the input and output 
disturbances. The designer may only change through 
shaping H(s) the magnitude in dB of GH at so fast 
a rate. Thus, the frequency bands where the magni- 
tude of the sensitivity function and complimentary 
sensitivity function may be small must be separated 
in frequency by a crossover region of a certain width. 
This width is dependent on G(s ) as well as on how 
small |C(s)| and |S(j)| must be. 

The trade-off between rejection of input and output 
disturbances through feedback is also inherent in 
passive isolation systems. Suppose we are capable of 
choosing the umbilical stiffness and damping of Fig. 2 
so as to design a passive isolator. Note that the 
transfer function relations 

-V(s)_ G(s) X(s) 1 

P(j) ]+G(j) D(s)/m ~ T + G(s) ^ 

apply where 



From this, it is easy to see that direct disturbances act 
as output disturbances while wall accelerations act as 
input disturbances. The difference between designing 
an isolation mount for base disturbances and for 
direct disturbances is well known and understood by 
vibration engineers. A soft mount is appropriate for 


isolating against base disturbances while a stiff mount 
is appropriate for direct disturbances excitation The 
loop shaping capability of springs and dampers is 
however, very restricted. Indeed, one cannot shape 
the loop to yield an unstable system. An active 
control system may have its loop shaped to an 
arbitrary specification provided it is possible to meet 
the specification without sacrificing system stability. 
Here lies the chief advantage of designing an isolation 
system from a control paradigm: the interaction of 
the conflicting specifications, stability and robustness, 
is clear throughout the loop shaping procedure. It 
should be noted here that sensitivity and complimen- 
tary sensitivity functions are extendable to MIMO 
systems through the use of singular values. 

Robustness in single-input-single-output con- 
troller design is measured by gain and phase margins 
The gam margin is the range of gain that can be 
introduced into the loop while maintaining stability. 
Similarly, the phase margin is the amount of phase 
that can be introduced into the loop while maintain- 
ing stability. The practical importance of the margins 
is that the gain and phase of the nominal plant is not 
the same as that of actual plant. These margins may 
be easily determined from Nyquist or Bode plots. 
Loop shaping also implies that a compensator H(s) 
should not be so large as to extend the crossover 
frequency of the compensated system into the higher 
frequency range where nominal models are very 
inaccurate. J 

Robustness for MIMO systems can also be 
specified in terms of the simultaneous gain and phase 
variations that may be introduced into the loops 
while preserving stability. However, this description 
does not account for unmodeled coupling in the 
dynamics. Uncertainty may be represented in terms 
of an additive (in parallel) or multiplicative (in series) 
transfer function matrix appended to the plant. 
(While these are the most common there are other 
representations.) Using either uncertainty represen- 
tation it can be easily shown by the small gain 
theorem that stability can be guaranteed if uncertain- 
ties in the plant are required to be bounded by a norm 
of the compensated plant. This is best represented in 
terms of the frequency dependent singular values of 
the plant and uncertainty transfer function matrices. 
This measure, however, is conservative since it allows 
cross coupling dynamics between channels that in 
actuality could never occur. The structured singular 
value methodology attempts to alleviate this conser- 
vatism through structuring the uncertainty model. 
Readers interested in a general trea'ment of MIMO 
stability and robustness should consult Ref. [9], 

3. PASSIVE ISOLATION: AN ANALOGY 

We now examine the design of an active vibration 
isolation system for microgravity space experiments 
from an analogy to passive isolators. Indeed, the 
primary reason for pursuing an active rather than a 
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Frequency (Hz) 

Fig. 5. Specifications (1) and (2) and uncompensated trans- 
missibility A*(.y)/T(.y). 

passive system is not the increased flexibility in loop 
shaping but the limitations of passive systems in 
attaining a stiffness low enough to meet the isolation 
requirements. This is true even when no umbilical is 
present. 

For the generic system model of eqn (1) with the 
nominal values m = 220 kg, fc=20N/m and 
c — 6.63 N • s/m (5% of critical damping). The trans- 
missibility curve between base and experiment accel- 
eration, shown in Fig. 5, is given by 

X(s) = 2 Z<o u s+<ol 

T(5) s 2 + 2£co n s + oj 2 n J) 


with 


co n = yjk jm — 0.3 rad/s = 0.048 Hz 


C = y/c 2 /4mk ^ 0.05 . 

Also depicted in Fig. 5 are the transmissibility 
specifications (1) and (2) discussed in Section 1. 
While the system satisfies the unity transmissiblity 
criterion, note that the natural frequency is not low 
enough to meet the 40 dB attenuation requirement. 
The system is also deficient in the magnification of 
disturbances at and near the resonance. Clearly any 
modification to the umbilical’s dynamics through 
feedback should include increased damping through 
a positive gain on experiment velocity. Feedback of 
inertial experiment velocity permits the damping co- 
efficient ( to be increased in the denominator of 
eqn (6) without changing it in the numerator. 
Thus, the resonance can be removed without affecting 
the roll-off rate [since the zero of eqn (6) is not 
changed]. 

If the umbilical were softer, say with k — 0.20 N/m, 
both specifications (1) and (2) could be met by the 
passive system. Unfortunately, a passive system 
cannot lower the stiffness with its inherently positive 
gains on position feedback. An active system, though, 
permits insertion of a negative stiffness spring in 
parallel with the umbilical. For example, for the 


nominal plant with the controller transfer functions 
of Fig. 3(a) equal to 


tfp(*)= -(6.0s + 19.8) 

7p(s)= “(6.0s + 19.8) (7) 

the natural frequency of the system is moved an 
order of magnitude lower. (Here, a negative damper 
has also been introduced so as to maintain the 
system’s 5% critical damping for the purpose of 
comparison. If less negative damping is introduced in 
order to remove the resonance, even more negative 
stiffness must be introduced to meet the 40 dB specifi- 
cation.) Note that this vibration engineering ap- 
proach, i.e. lowering the stiffness, requires the near 
cancellation of the umbilical’s stiffness with that 
introduced via feedback. If the negative stiffness 
exceeds that of the umbilical, the equivalent stiffness 
of the system will be negative and the system will be 
unstable. It is not surprising then that the introduc- 
tion of negative stiffness via the controller has no 
robustness whatsoever. The design using eqn (7) has 
<0.1° phase margin. The root locus for the system, 
shown in Fig. 6, clearly indicates this potential for 
instability. A focus on equivalent stiffness in isolation 
system design thus leads to control systems which 
sacrifice robustness for performance. In addition, a 
design which achieves isolation through lowering the 
system stiffness cannot attenuate direct disturbances 
over the same frequency band, as discussed in 
Section 2. 

From a vibration engineering viewpoint, an 
alternative means of achieving rejection of disturb- 
ances is to fasten the experiment rigidly to an inertial 
structure. While there is no such structure in space, 
it is possible to achieve this effect by a high positive 
gain feedback on inertial experiment position. (The 
inertial position must be obtained by integrating an 
accelerometer reading twice. This does pose a prob- 
lem since this procedure is marginally stable. How- 
ever, this problem may be ameliorated through 
replacing the integrators with a second order low pass 
filter. The authors are aware of this method being 
employed successfully on a six-degree-of-freedom 
magnetic suspension isolation rig at NASA Lewis 
Research Center.) This inertial position feedback acts 
as a very stiff spring tying the experiment to inertial 
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space. A controller and resulting transfer functions in 
this case are 

// p (s) = 2000 T p (s) = 0 

X(s ) _ 6.63s + 20 

?(*) 2205 2 -f 6.635 + 2020 

2205 2 


( 8 ) 


D(s)/m 220 5 2 + 6.635 + 2020 ‘ 

While this controller meets the 40 dB specification, it 
does not have unity transmissibility below 0.001 Hz. 
An experiment controlled in this fashion will collide 
into the wall. The feedforward transfer function may 
be adjusted to provide unit gain via 


T p (s) = 


-2000 
159 5 + 1 


This feedforward with the feedback term of eqn (8) 
effectively acts to base disturbances as a high relative 
stiffness up to 0.001 Hz changing to a large inertial 
stiffness at higher frequencies. The resulting trans- 
missibility X(S)fV(s) is presented in Fig. 7. Note that 
since the feedback loop introduces no damping, the 
original resonance is still present although ' less 
damped and at a higher frequency. This may be 
corrected by adding inertial damping into the feed- 
back loop. While this design method may be used to 
meet the specifications with robustness it has three 
faults: (1) it requires inertial experiment position, 
inertial wall position, and inertial experiment velocity 
measurements which are problematic to obtain, (2) it 
requires very high gains in both feedforward and 
feedback loops to obtain attenuation, and (3) an 
extension of the method to multi-degree-of-freedom 
systems would be difficult. It is also possible that 
when a flexible wall is considered, rather than the 
infinite impedance structure assumed, the system will 
be unstable. 

As another method of fastening the experiment to 
inertial space, one may employ inertial damping via 
feedback. By feeding back the inertial experiment 
velocity with a high gain, it is almost possible to 
achieve both the 40 dB and unity transmissibility 
specifications without resorting to feedforward. For 
example, with 

// p (5)= 10005 r p (5) = 0 




Fig. 8. Transmissibility X(s)/V(s) for inertial damping 
design. 

the resultant transmissibility is shown in Fig. 8. 
Unfortunately, the roll-off rate here is approx. 
20 dB/decade and therefore it is impossible to achieve 
both specifications simultaneously. This method has 
the advantage over the inertial spring of being a great 
deal simpler and requiring only one inertial measure- 
ment (experiment velocity which requires only one 
integration of accelerometer measurements). 

Another passive analogy is the lowering of the 
natural frequency of the umbilical by increasing the 
experiment mass. An increased experiment mass 
would attenuate direct disturbances as well as those 
transmitted through the umbilical. In addition, at 
frequencies below the natural frequency of the 
umbilical, the isolation system would have unity 
transmissibility. Of course, for space applications any 
additional mass is very costly. To lower the natural 
frequency by an order of magnitude would require 
increasing the experiment mass by a factor of one 
hundred. Clearly, it is not practical to accomplish 
increased isolation through the addition of real mass. 
However, it is possible to increase the effective mass 
of the system through feedback. This will be exam- 
ined in the next section, as this idea most properly 
evolves out of loop shaping. 

To summarize, the passive isolation analogy to 
active control system design yields some insight but 
falls short as a design tool on three counts: (1) it does 
not have the flexibility to shape the response with its 
simple analogical elements, stiffness, damping, and 
mass, so as to achieve the performance requirements, 
(2) it cannot be easily or effectively generalized to 
multi-degree-of-freedom problems, and (3) it com- 
pletely ignores the robustness problem inherent to 
active control systems. We advocate, therefore, that 
vibration engineers consider active isolation a con- 
trols problem and address it from an automatic 
controls perspective. 


Fig. 


7. Transmissibility 3t(sV?{s) for inertial stiffness with 
feedforward design. 


4. THE CONTROL SYSTEM APPROACH 

A simple controller is now designed for the system 
described by eqn (1) and the nominal values. The 
authors refer the reader again to Fig. 5 where the 
transmissibility curve between experiment and wall 
accelerations (or positions) is presented along with 
the design specifications (1) and (2). Pne goal is to 
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design a feedback control //„(*) that results in the 
closed loop transfer function 


<? d (j) = = - grWfj 

1 + G,(J)W P W 


(9) 


Real s 

Fig. II. Root locus of loop shaped design with respect to 
umbilical stiffness error. 

These two design specifications on H p (s) are shown 
in Fig. 9 along with a simple function satisfying these 
conditions, 


satisfying both constraints; i.e. 


H p (s)= 5000 s 2 . 


( 10 ) 


l G ciO"o)l % 1.0 ^ < 0.001 Hz 
2 n 


l G «t(M>)l <0.01 ^ > 0.1 Hz 

In 

Here, G p (s) and P(s) are as indicated in the block 
diagram of the system [Fig. 3(a)], Note that the 
uncontrolled system G p (s)P(s) already satisfies the 
first of these constraints; therefore, H p (s) should be 
very small in the low frequency band so that the 
closed loop system will continue to satisfy the unit 
transmissibility specification. Consequently, this 
specification yields a condition like 

\GH p (jio 0 )\ < 0.01 1 cu 0 

|tfp(M)l<0.2 /2^ <0-001 Hz. 

At and above 0.1 Hz, the attenuation of the uncon- 
trolled system is not sufficient. It is desirable to 
increase the attenuation by approximately two orders 
of magnitude. This may be accomplished by requiring 
H p (s) to be very large in this frequency range 
approximately * 

l^pO'Wo)! > I00{ w 0 

|tfp(M,)l>2000/2^ >01Hz - 


This controller design results in the closed loop 
transmissibility between experiment and wall acceler- 
ations which is plotted in Fig. 10. Note that both 
specifications (1) and (2) are met. Inertial damping 
should be added to this design to eliminate the 
resonance. It is easily seen from a root locus plot that 
tins design (S robust with respect to changes in 
umbilical/experiment properties (Fig. 11) an d 
actuator finite bandwidth (Fig. 12). In practice, this 
design would be improved by rolling off the controller 
gain. This limits the controller bandwidth so as to not 
affect possible unmodeled lightly-damped high fre- 
quency modes of the system (e g. wall flexure) A 
controller design would probably also include a weak 
position integral feedback to provide a slow centering 
force so that accelerometer bias and noise do not 
result in wall collision. 

The reader might object to the controller of eqn 
(10) since it is improper (i.e. has more zeros than 
poles). However, this controller is realizable. Note 
t . a ‘ / W mu l*iplies the position measurement to 
yield the control force. Since the factor r 3 in the time 
domain is equivalent to two differentations with 
respect to time, eqn (10) prescribes constant gain 
acceleration feedback. This, as discussed earlier, in- 
creases the effective mass of the system. [Of course, if 
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we modify eqn (10) to limit the controller band- 
width, then the mass analogy only holds within the 
band.] 

While both transmitted and direct disturbances are 
attenuated, the experiment acceleration level will be 
approximately the same as the accelerometer 
measurement noise level. This results from the trans- 
missibility between experiment acceleration and 
measurement noise being nearly one due to the high 
gain feedback. This is a fundamental issue as dis- 
cussed in Section 2; one must trade-off the rejection 
of disturbances to the system and the rejection of 
measurement noise. Since the disturbances may be up 
to 1000 times larger than the measurement noise 
(accelerometer resolution typically 1 ^g) the con- 
troller is designed to reject disturbances. The per- 
formance of the control system is thus directly a 
function of the quality of the accelerometer. 

Recently, an alternative approach to design of 
active vibration isolation control systems for micro- 
gravity experiments was presented in Ref. [10]. A 
desired transmissibility ratio G cX (s) is specified along 
with the plant model G p (s) and P(s). Equation (9) is 
then solved via algebraic manipulation for the feed- 
back controller H p (s) that yields the desired trans- 
missibility (feedback of relative position is also 
allowed and may be used; if used, a second condition 
must then be specified for solution). While this ap- 
proach resembles loop shaping in that it attempts to 
achieve a certain transmissibility, it is fundamentally 
different in that it does not properly consider the 
plant. The algebraic procedure in essence first elimin- 
ates the plant and then replaces it with one which will 
yield the desired transmissibility. As a control design 
procedure, this methodology has serious flaws: (l) the 
stability of the resulting system may be entirely 
dependent on perfect knowledge of the plant, (2) the 
procedure incorporates none of the known relation- 
ships and fundamental trade-offs between stability 
and attenuation; it implies that any specified trans- 
missibility is achievable, and (3) for systems with 
right half plane poles/zeros, the methodology may 
attempt cancellation with right half plane zeros/poles. 
For a simple controls problem, the algebraic manipu- 
lation method may result in a good controller. 
However, for more difficult problems, the method 
is questionable. An extension of this methodology 
to MIMO control would be plagued by many 
problems. 

To summarize, controller design for single-degree- 
of-freedom vibration isolation problems is best per- 
formed through the classical control framework of 
loop shaping where the natural interplay between 
performance, stability and robustness are evident. 
For multiple degree of freedom isolation problems, 
recent advances in controller design, such as the 
extension of loop shaping principles via frequency 
weighting and singular values [1 1] seems to be most 
promising. In order to emphasize the question of 
coordination in control of MIMO systems, we next 


examine a multiple degree of freedom isolation 
problem. 

5. A MULTIPLE-DEGREE-OF-FREEDOM system 

A common misunderstanding among many engin- 
eers unfamiliar with control system design is the 
nature of the differences between SISO and MIMO 
control problems. The relative ease with which the 
uninitiated comprehend the elimination of one error 
signal through negative error feedback yields the false 
impression that the MIMO control problem is little 
more than the feeding back of multiple error signals. 
This impression, however, is not totally groundless. 
Indeed, many MIMO controllers in use today were 
designed by a single-loop-at-a-time procedure. 
Design with this method can be quite difficult, time 
consuming, and non-intuitive. Robustness is difficult 
to check except by analyzing all the possible permu- 
tations to the nominal plant. The fundamental 
problem in MIMO design is the coordination of the 
control in coupled channels when the plant is not well 
know'n (poorly modeled or time varying). 

Easily decoupled active vibration isolation control 
problems may be deceptively simple. Unmodeled 
cross-coupling due to inaccuracies in center of mass, 
sensor, and/or umbilical locations can result in poor 
performance and even instability. An example iso- 
lation problem illustrates. Figure 13 shows a two- 
degree-of-freedom isolation system composed of an 
isolated platform (width 0.5 m and height 0.2 m, 
depth unspecified), two accelerometers, two actua- 
tors, an umbilical, and a translating base. The plat- 
form may translate vertically or rotate about its 
center of mass. The actuators and accelerometers are 
positioned a distance of q = 0.2 m symmetrically 
about the assumed center of mass location. An 
umbilical of stiffness k (no damping) runs between 
this location and the base. The platform has mass nt 
and inertia /. The equations of motion for the 
platform's translation jt(r) and rotation 6(t) are 

mx + k A6 +kx=f x +f 2 + 
fff + k& 2 0 +kAx = (q + AV 2 -( 9 - + d 2 (II) 



a Acc«!«romet0f 

Fig. 13. Two-degree-of-freedom active isolation system. 
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where d x and d 2 are the disturbances, and A is the 
error in the assumed center of mass. The accelerom- 
eter readings are 

Mi = X - (g - A)ff 

Mi = x +(q + A)ff. (12) 

The nominal system (A = 0) can be decoupled in 
terms of the degrees of freedom by the change in 
variables 

*■-/. +/i 

i) 

z i = (yi + vi )/2 

z 2 = ?(y2-y,)/2 (13) 

which are nominally the translational force, the 
moment, the translational acceleration, and the angu- 
lar acceleration for the platform, respectively. The 
nominal transfer functions for the system are then 

z,w “[s2Tl}' w + D '<'» 

Z2(J)=|^](A/(j) + Dj(5)). 

For translational mot ion, the natural frequency of 
the platform is yjkjm. The rotational motion of the 
platform is free since the umbilical is attached to the 


translational motion from 0.056 to 0.006 Hz with 
40% of critical damping. The controlled rotational 
motion has a natural frequency of 0.006 Hz with 26% 
of critical damping. This controller design would 
yield very effective isolation on the nominal system. 

The actual close loop transfer functions, however, 
will be different from the nominal due to the error in 
the center of mass, A. The transmissibility can be 
derived from eqns (1 1 HI 5) as follows 

[ms* + k\X(s) + [*A]0 (j) = F(s) + £,(,) 

[/s J + fcA 2 ]©^) + [k AJA'(s) = M(s) + A F(s) + D 2 (s) 

Z,(s) = [s 2 ]X(s) + lAs 1 ]e(s) 

= [■s 2 ]0(j) 


F(s)= ~[a+c/s)Z,(s) 
M(s)= — [n/j + b ls 2 ]Z l {s) 

yielding 


[ (m+fl)5 2 + c\y + A:l 

~p \Z\ fa) 4- [m &]Z 2 (s) = Z), ( 5 ) 

[ (as^s+k^ l 

l s 2 + [ J2 JZj(i) = Z)j(s). 


The poles of this system are given by the roots of the 
characteristic equation 


center of mass. To compensate the nominal system, 
feedback can be designed for each mode of the system 
separately, since the system is decoupled. Transla- 
tional acceleration and velocity feedback are first 
used to add effective mass and damping. 


F(j)= ^Z,(s). (14) 

This lowers the natural frequency of translational 
motion yielding the closed loop transfer function 


Z,(s) 


(m + +k 



Next, angular deflection feedback is used to constrain 
low frequency rotational motion and some damping 
is provided. 


yielding / 


05 ) 


[(m + a)s 2 + cs + k][/s J + ns + b] 

-[mA][A(as 2 -(-cj +k)] = 0 (16) 

For the nominal plant, A = 0, the roots of eqn (16) 
result in the prescribed natural frequencies and criti- 
cal dampings. However, as the center of mass error 
increases, the poles migrate and the system becomes 
unstable. For an error as small as 6 mm, instability 
occurs. A plot of the pole movement vs error in 
center of mass is shown in Fig. 14. This sensitivity 
results from the ill-conditioned character of the re- 
quired controller. Ill-conditioned here means that the 
controller s gain to an output signal varies strongly 
with the signal’s direction. This results in a control 
system which is not robust to this model’s uncertainty 
(center of mass) [12). A proper MIMO controller 
design might remedy this problem. In any case, an 
analysis of the problem from a MIMO control 


z 2( s)=r - > 2 ,) 

[_« + ns + 6J 

The following values are used to illustrate this 
example 

Platform Control system 

m = 400 kg a = 31600 kg 

k = SO N/m c = 1000 N ■ s/m 

/= 10 kg • m 2 b = 0.015 N • m 

n — 0.2 N • m • s 

where the control system values are in effective units. 
This control design lowers the natural frequency of 



Fig. 14. Root locus of two-degrec-of*frcedom with respect 
to center of mass error A. 




696 


C. R. Knospe el al. 


perspective would indicate the potential instability 
and the nature of the trade-off between performance 
and robustness. (The authors note that increasing the 
damping and stiffness for the rotational mode 
improves the system robustness significantly, while 
changing the damping or effective mass for the 
translational mode has little effect.) 


6. LINEAR QUADRATIC REGULATOR FOR ISOLATION 


MIMO control design, since it requires a high 
degree of coordination, must proceed by a synthesis 
procedure. One such method is Linear Quadratic 
Regulator (LQR) synthesis [13]. This produces a state 
feedback controller which is optimal with respect to 
the quadratic (two norm) performance function 


J = 


*00 

x J (joj)Qx(jaj) + u J (j(o)Ru{joj) dco 

J - X 


( 17 ) 


where Q and R are respectively the symmetric 
(usually diagonal) state and control weighting 
matrices, and x(jco) and u(ja>) are the Fourier 
transforms of the state and control vectors. The state 
(positions and velocities for vibration isolation) sa- 
tisfies the differential equation 


k — Ax + Bu . 


The quadratic performance function of LQR, eqn 
(17), is well suited to this problem since vibration 
isolation quality is usually measured in terms of 
root-mean-square. However, some modification of 
the performance function is necessary to apply this 
synthesis procedure to microgravity isolation con- 
troller design. The reader will note that state feedback 
for the isolation problem is feedback of experiment 
positions, velocities, angles, and angular velocities. 
Thus, LQR can only result in inertial stiffness and 
inertial damping feedback. As was shown in Section 
3, these isolation techniques cannot yield acceptable 
isolation performance. Thus, an LQR performance 
function of the form of eqn (17) will not yield a 
satisfactory controller. Note that the differential 
equation does not include a disturbance term. Conse- 
quently, the controller is optimal with respect to 
white noise. Since the power spectrum of the micro- 
gravity environment is not of this shape, the LQR 
controller will not be optimal with respect to rejection 
of the disturbance. Through the incorporation of a 
disturbance model (essentially a shaping filter) the 
LQR problem may be modified to yield an optimal 
disturbance accommodating (i.e. rejection) con- 
troller. This incorporates the addition of pseudo- 
states to the state variable model [14]. 

Closely related to disturbance accommodation is 
the concept of frequency weighted LQR performance 
functions [15]. Here, the Q and R matrices are chosen 
to be even rational functions of frequency. This 
results in the addition of pseudo-states to the state 
variable model. Through choice of the weighting 
functions, the designer can in essence shape the 


control loops [1 1]. This also permits the weighting of 
experiment acceleration. It should be noted that for 
successful application of LQR theory to the micro- 
gravity isolation problem frequency-shaped cost 
functions must be used. Without this, the control 
resulting from the synthesis procedure would attenu- 
ate the vibration at frequencies below 0.001 Hz (non- 
unity transmissibility). The reader should note that 
the well known robustness characteristics of LQR 
controllers do not apply to most frequency shaped 
designs or to plants with unmodeled cross coupling. 

7. CONCLUSIONS 

Successful active isolation for microgravity exper- 
iments can be achieved, but only if the problem is 
analyzed from a controls perspective. A passive iso- 
lation analogy, while useful for an understanding of 
the control problem, is not an effective design tool. 
Design of active vibration control systems can best be 
carried out through loop shaping. For intrusive iso- 
lation platforms, this results in a high gain accelera- 
tion feedback design. A two-degree-of-freedom 
example was used to illustrate the instability that can 
result under unmodeled cross coupling when the 
control system is designed via decoupling/single loop 
design procedures. The source of this sensitivity was 
ill-conditioning of the controller. The LQR was 
examined for the isolation problem. For synthesis of 
an effective controller, the procedure must be 
modified to include loop shaping. 
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APPENDIX 

Nomenclature 

A - system dynamic matrix 
a = acceleration feedback coefficient 
B = system input matrix 
b = rotational stiffness feedback coefficient 
C(j) = complimentary sensitivity function 
c = damping 


d,D(s) - direct disturbance force 
5{s) = equivalent disturbance 
fF,F(s) = actuator force 

G(s) = plant transfer function 
H(s) = feedback transfer function 
/ = moment of inertia 
k = stiffness 

A/,A/(0 = control moment 

n = rotational damping feedback 
P(s) = umbilical precompensation transfer function 
Q = state weighting matrix 
q = actuator placement 
R — control weighting matrix 
SCO = sensitivity function 
7X0 = feedforward transfer function 
u = control vector 
v(s) = measurement noise 
*,Y(0 = experiment position 
x = state vector 
y t YCO^wall position 

y - accelerometer measurements 
z,Z(s) = decoupled measurements 
A = center of mass error 
0,@(O = angular position 
cu n = natural frequency 
C = percent critical damping 


Subscripts 

a - acceleration 
cl — closed loop 
p = position. 
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Abstract 

Manned orbiters will require active 
vibration isolation for acceleration-sensitive 
microgravity science experiments. This 
paper investigates control aspects of the 
isolation problem and proposes a viable, 
robust control. Since umbilicals are highly 
desirable for many experiments, and since 
their presence greatly affects the vibration 
isolation problem, they must be considered 
in control synthesis. Experiment isolation 
involves reducing undesirable plant outputs, 
such as excessive accelerations and 
unacceptable relative displacements. The 
former are undesirable due to experiment 
dynamic environmental demands; the latter, 
due to rattlespace constraints. Ideally, 
controller design should minimize these 
outputs, while considering input— and 
output directionality and frequency content. 
This paper investigates applying modern 
control theory to the isolation problem, 
incorporating frequency— weighting and 
disturbance accommodation techniques. The 
resulting controller achieves excellent system 
performance, for plants within a reasonable 
range of variations from the nominal. 

Robust stability and performance guarantees 
were measured by singular value and 
structured singular value checks, yielding 
guarantees on allowable real plant parameter 
uncertainties, and on acceptable controller 
and sensor phase and gain variations. The 
problem of unmodeled high frequency modes 
was eliminated by using frequency weighting 
to reduce controller bandwidth. 


f-' 

This design method has been developed 
with the six-aegree-of— freedom isolation 
problem in mind, to which it is fully 
applicable. In this paper the method is 
applied successfully to the single— degree— 
of— freedom isolation problem. 

Introduction 

The vibration environment onboard 
current and planned manned orbiters will 
require isolation for many microgravity 
science experiments^ The disturbance 
frequencies are sufficiently low, and the 
attenuation requirements sufficiently great, 
so as to preclude a solely passive isolation 
systenfr 3 ! 

Since the disturbances to be attenuated 
are three-dimensional (3— D)0} the isolation 
actuator(s) must be capable of acting over 
six degrees of freedom. Although 
microgravity isolation systems have been 
developed and tested^} most controllers 
offered to date fail to take into account the 
effect of umbilicals in a 3— D isolation system 
design^. Since umbilicals are highly 
desirable for many experiments®(e.g., for 
evacuation, power transmission, coming, 
material transport), they must be considered 
in a generally applicable control scheme. 
Simple application of nonintrusive control 
methodologies is insufficient^ 

The requisite multi ple-degree- 
of— freedom (MDOF) controller is much 
more difficult to design than a 
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single— degree-of— freedom (SDOF) 
controller, because the isolation system has 
many inputs (actuator forces! and outputs 
(measured displacements ana accelerations). 
Multiple-^nput— multiple— output (MIMO) 
designs can be very susceptible to unmodeled 
cross— coupling between channels of input or 
output)®} a problem not encountered in 
SDOF design. The control forces used must 
therefore be properly coordinated if the 
controller’s performance is to be sufficiently 
insensitive to unmodeled dynamics (i.e., 
robust). The design of a robust MIMO 
control system is a nontrivial problem 
requiring the iterative use of synthesis and 
analysis tools, the former for controller 
design and the latter for system stability — 
and performance evaluation®! 

H 2 synthesis, the controls approach 
chosen for this work, has been used 
effectively in MIMO control problem^ 7 ! 
Additionally, the approach seems 
appropriate for the present problem, since 
the quadratic performance index relates well 
to root— mean— square (rms) statistics and 
power spectral densities, in which form 
present orbiter acceleration data is currently 
made available#} 

The microgravity isolation problem 
aboard manned orbiters involves reduction 
of different kinds of undesirable plant 
outputs, such as excessive accelerations and 
unacceptable relative displacements. The 
former outputs are undesirable due to the 
demands of the experiments themselves®); 
the latter, due to rattlespace constraint^ 3 * 8 ! 
Since some disturbances may be directional, 
since some undesirable plant outputs may be 
more important than others, and since the 
degree of undesirability of these outputs may 
also vary with direction or frequency, the 
design of an optimal controller ideally should 
incorporate these factors. 

Plant outputs cannot be minimized 
apart from consideration of the associated 
control costs, because active control both 
consumes power and releases heat. Since 
both of these costs are of concern in a space 
environment, the control effort used should 
not be excessive. And at higher frequencies, 
the control effort should be reduced in order 


to limit controller bandwidth for the sake of 
robustness concerns)®! 

H 2 synthesis, more commonly known as 
LQG (Linear Quadratic Gaussian) or LQR 
(Linear Quadratic Regulator) synthesis, 
allows the designer to develop an optimal 
regulator that consists of full state feedback 
using constant feedback gains. Inaccessible 
states are reconstructed by an asymptotic 
observer that uses constant observer gains. 

In the LQR synthesis approach the observer 
gains are chosen to produce a stable observer 
with poles placed as desired by the 
designer# 0 ! In the LQG synthesis approach 
the observer gains are selected to minimize 
the rms error of the observed states from the 
actual states, based on the assumption that 
the process— and measurement noise vectors 
are zero-mean white Gaussian# 0 ! 

H 2 synthesis can be extended to allow 
the controller (i.e., the regulator and 
observer) gains to be optimized for colored 
process noise# 0 * ll ! Another extension 
permits the states and controls 
independently to be "frequency-weighted", 
so that certain frequencies of each are more 
heavily penalized than others, in an rms 
sense# 2 * 11 ! A third extension allows for the 
inclusion of deterministic disturbance 
information into the optimized control 
expression C 13 » 14 > 15 1 

The application of basic H 2 synthesis to 
the tethered microgravity isolation problem 
was first proposed in 1990#°); but the most 
complete treatment to date, of the full 
extended H 2 synthesis approach, appeared a 
year later# 7 ! An algebraic introduction is 
published in reference (18), and the full 
deterministic solution is developed in (19). 
The present paper presents the results of a 
systematic application of the extended H 2 
synthesis method to a realistic 1-D isolation 
problem. The resulting controlled system is 
evaluated using singular values and 
structured singular values ("/i analysis") to 
determine guarantees of system performance 
with uncertainties in the umbilical, payload, 
sensor, and actuator models. An excellent 
introduction to much of the pertinent 
analysis methodology is contained in 
reference (20). 


2 


Problem Description 


The general 3-D vibration isolation 
problem for the tethered payload has three 
translational and three rotational degrees of 
freedom. In reference (20) Allan and Knospe 
presented a brief survey of several published 
3-D syspension designs. The extended H a 
synthesis - n analysis approach can be 
readily applied to a 3-D problem once the 
system model has been reduced to state- 
space form (i.e., consisting of a set of 1st 
order linear ordinary differential equations) 
However, the simpler 1-D problem offers the 
benefit of providing a simple model that is 
highly conducive to developing a physical 
intuition. Further, the specific 
mathematical model is much less 

d e pcndent. The mathematical 
theory summarized below is generally 
applicable to either the 1-D or 3-D 
problem; but for the reasons noted above, it 
will be applied only to the 1-D problem in 
the application section of this paper. 


Let the payload (e.g., experiment) be 
modeled as a lumped mass with inertial 
position x(t). Assume that the orbiter (i.e., 
experiment rack) has inertial position d(t) 
and that massless umbilicals characterized 
by a stiffness and a damping connect the 
payload and orbiter. Suppose further that a 
Lorentz actuator applies a control force 
proportional to the applied current u(t) with 
proportionality constant a. Such a model is 
shown in Figure 1. Typical parameter 



?„ W , e wy? ho ! en: mass = 75 lbm > stiffness 
// 'nlLT 1, damping = 0.01138 lbf-sec/ft 
(C = 0.3%), and a = 2 lbf/Amp. Clearly 
this model is inaccurate, since the plant will 
n0 V* a perfectly rigid mass and since the 
umbilical is neither massless nor accurately 
characterizable by a single stiffness-and- 
damping model. The differences between 
,® actual system and the nominal system 
will be addressed later. It will be seen that 
the system (consisting of payload, umbilical, 
sensor(s), controller, and actuator) is robust 
to anticipated uncertainties (i.e., remains 
stable and achieves good performance). 


The goal is to find a feedback controller 
such that the controlled system satisfies the 
following specifications: 


1. Above 0.1 Hz the payload acceleration 
x(t) should be 40 dB below the orbiter 
acceleration d(t). 

2. Below 0.001 Hz the payload vibration 

/ \ track the orbiter vibration 

d(t) to within 10 percent, in order to 
prevent collision of the payload with the 
walls of the experiment rack®) 

3. The payload should track perfectly the 
DC motion of the orbiter, where no 
relative motion can be tolerated. 

4. The loop gain of the system (plant and 
controller) should be less than 0.1 above 
200 Hz, to avoid controller excitation of 
orbiter- or payload flexible modes. 

5. The system should remain stable and 
exhibit good performance for anticipated 
inaccuracies in the system model. 

These specifications are derived from the 
available information concerning existing 
and required vibration environments, (Ownile 
recognizing that rattlespace constraints will 
require orbiter tracking at very low 
frequencies 


Figure 2 shows the open loop 
transmissibilitv of the plant (i.e., its 
amplification factor, plotted as a function of 
harmonic-input frequency), along with 
specifications (1) and (2). 
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Basic Solution 

The equation of motion can be expressed 
in the state-space form 

x = Ax + B u + E s w*. 

X = Cx + Du + Mri* 
where w., and n* are white noise vectors 

reflecting the power spectral densities of d 
and of the sensor noise, and where the vector 
x includes states derived from the equation 
of motion (e.g., inertial or relative positions 
or velocities, or accelerations), along with 
disturbance-accommodation pseudostates 
and frequency weighting pseudostates if 
desired^ 7 ! 

As outlined in reference (17), once the 
problem is expressed in the above state 
space form with a suitably chosen quadratic 
performance index, an optimal controller can 
be found readily by the solution of an 
algebraic Riccati equation. The solution and 
solution methods are well-known. The 
resulting controller uses a control current 
that employs full state feedback, where 
inaccessible states (or pseudostates) may be 
reconstructed using an asymptotic observer. 
If the type of asymptotic observer used is a 
Kalman-Bucy filter, the observer design also 
involves the straightforward solution of 
another algebraic Riccati equation*! 7 !. Thus, 
a controller can be synthesized in a 
straightforward manner. Considerable 
insight is necessary, however, in determining 


how to use the powerful tools of 
disturbance-accommodation and frequency 
weighting to one’s advantage. Once a 
controller has been synthesized, it must then 
be evaluated to see how it fares against the 
design specifications. This check is 
straightforward for specifications #1, 2, 3, 
and 4. However, for specification #5 the 
check is much more difficult. 

H analysis is a powerful tool that can be 
used to determine how much uncertainty can 
be tolerated at various locations in the 
closed loop system. The measures are 
conservative, due to aspects of the 
mathematics involved; but the results 
provide guaranteed minima of allowable 
uncertainty magnitudes that can be of 
immense value to the designer. In brief, the 
designer places one or more complex 
uncertainty blocks A(s) at appropriate 
locations in the system transfer function 
block diagram. By /x analysis methods he 
can then determine how large (2— norm) the 
uncertainty (or uncertainties) can be 
without driving the system unstable or 
exceeding specified performance limits. 
Typically uncertainty blocks fall into the 
categories of multiplicative input, 
multiplicative output, additive, and 
performance uncertainties. These can be 
used to provide conservative measures of 
allowable actuator gain and phase 
variations, sensor gain and phase variations, 
and unmodeled higher frequency plant 
dynamics, along with simultaneous 
performance guarantees. Feedback 
uncertainty blocks can be used as well to 
provide measures of allowable plant 
parameter variations from nominal values. 

The rest of this paper is primarily 
devoted to helping provide insight into how 
to use extended synthesis and /i analysis 
effectively, for microgravity isolation, by 
reviewing a logical application of the 
approach to a realistic 1-D problem. 

Controller Design 

To begin the design process, the 
one— dimensional problem was first expressed 
in state-space form, with payload relative 
position, relative velocity, and acceleration 
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selected as states. Although many other 
state choices could have been made, these 
three were chosen to minimize the number of 
states necessary and to maximize the 
physical intuition possible. The selection 
would result in a state feedback control that 
respectively modifies the effective umbilical 
stiffness and damping, and the effective 
payload mass all being familiar, accessible, 
and intuitive system parameters. Relative 
rather than inertial, position feedback would 
help to avoid exceeding rattlespace limits; 
and relative velocity feedback would provide 
a means of damping out system resonances. 

1 he selection of acceleration as a state was 
considered desirable due to insight gained 
from an analogy to passive isolation systems. 

A controller which increases effective 
payload mass (by negative acceleration 
feedback) would potentially be able to 
accomplish disturbance rejection without 
unnecessarily sacrificing stability- or 
performance robustness. 

A second important feature of the 
problem formulation was the decision to 
incorporate disturbances of two different 
kinds, the direct (i.e., onboard the 
experiment) and the indirect (i.e., acting via 
the umbilical). It had been observed by 
examining the pertinent transfer functions, 
^V^mg the effective umbilical stiffness 
could aid in indirect disturbance rejection 
oidy, but that increasing payload effective 
mass could help reject disturbances of both 
kinds. Although the primary type of 
disturbance was considered likely to be the 

\n£ eC l , vx? m , e f nS 1 w * s needed t0 force ^e 
LQR-KBF (also known as LQG) design 

machinery" to increase effective mass so as 

to result m a robust controller. Including a 

direct disturbance provided this mechanism. 

After completing the problem 
formulation, the next step was to develop a 
computer code for use in design and analysis. 

iSrA-TT A^ ed desi S n code w as written in 
MATLAB to allow for frequency accom- 
modation of both direct and indirect 
disturbances. A large selection of frequency 
weightings and disturbance accommodation 
filters was made available to the designer, 
ihe code computes both feedback and 
observer gains. A number of analysis 


routines were also written to allow the 
designer to evaluate the resultant designs for 
purposes of comparison. The number of 
system states, system performance, stability 
robustness, parameter sensitivity, and 7 
observer quality were items whose 
comparisons were facilitated by these 
routines. 

. W j‘ h the desi gn and analysis tools in 
place, the next step was to develop the 
desired controller. In order to make the 
controUer as simple as possible, it was 
decided to begin with the basic LQG 
approach (no frequency weighting, no 
disturbance accommodation, no direct 
disturbance) and to add complexity as 
needed. At each stage of additional 
complexity an iterative cycle of design and 
analyas was employed in an attempt to get 

ofcompl l eid a ty Vab ‘ e C ° mr0ller “ lhat level 

The basic LQG approach yielded a 
satisfactory controller in terms of 

mWt manC f’ b u t U had almost no stability 
robustness to changes in umbilical stiffness 

from the nominal (as measured by feedback 
uncertainty). This lack of robustness wm 
due to the fact that LQG found adding 

fn e | atl ?i t0 ** a "Reaper" means of 
indirect disturbance rejection than adding 

effective mass. No frequency weighting was 
found which could rectify this problem 6 

A direct white disturbance was added in 
an attempt to force the LQG design 
"machinery" to add effective mass 
Although there were some gains in stability 
robustness this was due entirely to changes 
m observer gain matrix L. The feedback 
gain matrix K remained unaffected (note 
that this is fundamental in LQG theory and 
is not a numerical problem), and the 
feedback stability robustness was still 
unsatisfactory. 

Disturbance accommodation, with a 
lowpass filter applied to a large direct 
(white) disturbance, resulted in a controller 
with excellent feedback- and multiplicative 
input stability robustnesses, as measured bv 
singular value checks. The multiplicative 
output stability robustness was unacceptably 
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low if cross-coupling was considered possible 
between states, but structured singular value 
checks indicated that without cross-coupling 
the allowable multiplicative output 
uncertainty was quite satisfactory. Since 
effective stiffness, effective damping, and 
effective mass of the controlled system are 
uncoupled for the true one-dimensional 
problem, the stability robustness measures of 
the system were considered acceptable. 
Further, the performance was excellent, 
easily exceeding the specifications. However 
the controller gains were still large at higher’ 
frequencies where unmodeled system modes 
were of concern (see specification # 4). It 
was therefore necessary to use state— and 
control frequency weighting to force the 
controller to "turn off' by approximately 
100 Hz (i.e., to reduce loop gain below a 
magnitude of one) so as to avoid exciting 
unmodeled flexible modes. 

To reduce the loop gain at the higher 
frequencies it was necessary in that range (1) 
to place a high weight on control, (2) to 
apply low weights to all three states, and (3) 
to reduce the direct disturbance. At low 
frequencies the control weighting was left 
constant (i.e., "flat"), in an attempt to 
minimize the number of added pseudostates. 
However, the resulting closed loop system 
now had very poor low frequency stability 
robustness to parametric uncertainties, even 
though it both retained its excellent 
performance and now provided the desired 
low controller bandwidth. 

A classical design approach to the 
problem provided a simple solution to the 
robustness issue. It was noted that for a s 
controller with acceptable nominal f 

performance the low frequency asymptote for ° 
controller gain could have slope -1 or 0 or 
greater (Bode-a, log-log scale). That is, 
control gain at DC could be zero (slope >1) 
finite (slope = 01, or infinite (slope — —1). ’ 
Zero DC controller gain would, of course, 
result in a closed loop system that would’ I 
achieve the unit transmissibility of the open j 
loop system at low frequencies, as desired. “• - 
But by using a control weighting filter with 
zero DC gain (slope ^ 1) the extended H 2 
synthesis "machinery" could be freed to 
consider finite or infinite DC controller-gain 


options. Consequently the control weighting 
was made to be zero at DC (at the expense 

» pseudostate). The result was a 
controller that satisfied the design 
specifications and exhibited good stability 
robustness to parametric and to 
multiplicative input- and output 
uncertainties. Considering (for the moment) 
only single-parameter uncertainties 
stability was guaranteed for umbilical 
stiffness to within ±99.7% of nominal, and 
umbilical damping could be essentially 
unknown. Payload mass needed to be 
known only to within ±65.2% of nominal. 
Having these initial favorable indicators of 
system robustness the next step was to 
reduce the controller size. Further 
robustness analysis would then be conducted 
on the reduced— order controller. 

The controller described above was a 
mnth^rder coutrdlef (i.e., had nine states), 
with payload acceleration as its only 
required input. Other states and 
pseudostates were reconstructed in the 
observer. To reduce the controller to a 
^^r order, a routine was written in 
MATLAB in order to permit removing high 
frequency modes (modal truncation) and 
weakly controllable and -observable system 
dynamic^) The result of applying this to 
the ninth— order controller was a third-order 
controller that has all the essential features 
of the ninth-order one. The loop gain 
controller, and transmissibility plots for this 
reduced controller are shown in Figures 
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Figure 3a.— open loop bode plot, froe control 
to -control (negative loop gain) 
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systee transnissibilities 

that the transmissibility is unity up to 10* 3 
Hz and that it is below 10 ‘ 4 at 0 1 Hz 

lmn ,he °P en >°°P and closed 
loop Bode plots merge at about 100 Hz. 

This is due to the fact that the controller has 
“u'fe 3b U r d o£f ' by that &e,U “ Cy 

There are four basic checks that must be 
made of any controlled system: nominal 
stability, nominal performance, robust 
stability, and robust performance. These 
lour checks are considered below 
consecutively. ’ 

, tended H 2 synthesis method used 
lor this portion provides an inherent 
guarantee of stability for a nominal plant 
with full state feedback. Further, the 

P rind P Ie " guarantees that for a 
perfectly known plant a stable asymptotic 


observer will not destabilize the system. 

J.n ®’ ?° r rm . naI stability is assured with the 

fSi?- der *°i Serv ® r - P rov ided the observer 
itself is stable. Reducing the controller 

2S? re 1 mov ® 8 this guarantee, but simple 
eigenvalue checks verify that both the 
j!iJ d third-order controller as designed 
and the associated controlled system <ure 

5 th!? in r the . no ? u ° al piant. A simple check 
of the loop gam Bode plot (Figure 3a) 

confirms the conclusion that the closed loop 
system is stable, since it is known that the 
loop gain is minimum phase. 

The second necessary check is of 
nominal performance. As indicated by the 
closed-loop transmissibility plot (Figure 3c) 
the nominal performance is quite y 

satisfactory Note that the "less than 10-3" 
spec at 0.1 Hz is surpassed by more than an 
order of magnitude. This overdesign was 
intentional, and necessary, since plant 
modeling errors (open loop system sensors 
and actuators) will certainly degS’de ’ 
performance margins. 6 

tn bU - t 8ta t m *y mea sures are necessary 
to determine whether the closed-loop system 
will remain stable given the anticipated^ 
sensor, actuator, and plant parameter 
imcertainties. Three different types of 
robust stability measures were used for 
guaranteeing system stability for ’ 

^V P ^k atl u e input > multiplicative output 
and feedback uncertainties below certain ’ 

ilpf 1S ‘M Tlie m ^ tlpIlcati ve input uncertainty 
allowable was found to be equivalent to a 
guaranteed phase margin (interval) of Ms* 

+48 J, and to a guaranteed gain margik ’ 
(interval) of (0.304, 5.434]. The actual 
margins are even larger (phase margins: 

1-55 , +55 J, gain margins: [0,+oofi. Since 
only one plant output is sensed (viz 
payload acceleration), the multiplicative 
input and output robust stability guarantees 
are identical. A feedback uncertainty 
measure was used to determine guaranteed 
mimmum stability bounds on uncertainties 
in umbilical stiffness and damping, and on 
payload mass. It was found, as noted 

P, re ™ us| y (P- 6), that closed-loop system 
stability was guaranteed for 

single-parameter uncertainties much larger 
than anticipated. It was found, by B 
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considering the feedback uncertainty 
structure, that for simultaneous mass, 
damping, and stiffness uncertainties of 
±20%, ±.100%, and ±69%, respectively, 
system stability could be assured. Higher 
frequency modes of the system were 
considered not to be a significant concern 
since the controller bandwidth was limited 
during design. 

Finally, measures were needed of 
performance robustness. Structured singular 
value plots were made to find conservative 
bounds on multiplicative input (and output) 
uncertainties that would not lead to plants 
with unacceptable performance. Below 10 ' 3 
Hz it was found that for combined sensor 
and actuator uncertainties of up to ±11* in 
phase or of ±19% in gain the performance 
can be guaranteed to remain acceptable. At 
higher frequencies the guarantees are much 
better, so that by 220 Hz uncertainties of up 
to ±180o in phase or of ±200% in gain are 
permissible. 

Structured singular value plots were also 
used in an attempt to find performance 
robustness guarantees in the face of known 
parametric uncertainties, but the effort was 
only partly successful. The checks led to the 
conclusion that for single-parameter 
uncertainties in stiffness of ±40% both 
stability and acceptable performance could 
be assured. However, single— parameter 
uncertainty bounds found by this method on 
damping and mass were too conservative to 
be useful. Consequently, real parametric 
studies were conducted on plant— uncertainty 
effects on closed— loop performance. It was 
found that closed loop performance appeared 
acceptable for the various combinations of 
parametric uncertainties examined, with 
mass and stiffness varied in the intervals 
[ — 50%, +100%] and [-20%, +100%], 
respectively, and witn damping varied by 
more than ten times its nominal value. 

Conclusion 

The above extended H 2 synthesis - /x 
analysis approach was found to produce a 
controller that easily satisfies the competing 
demands of the posed 1— D microgravity 
vibration isolation problem. Further, unlike 


the classical approach, it is readily 
extendable for use on a 3— D problem. 
Frequency weighting and disturbance- 
accommodation were both found to be 
necessary if H 2 synthesis is to be used in 
involving the posed isolation problem. Their 
inclusion, along with a judicious choice of 
states, provides the designer with a powerful 
and intuitive set of weapons for his design 
arsenal. Disturbance accommodation of a 
direct disturbance model was found to be 
effective in forcing the Hj synthesis 
machinery to avoid negative-stiffness 
solutions. The result was an actively 
controlled system that uses a "smart" form 
of acceleration feedback to overcome the 
robustness problems that commonly plague 
the basic LQG synthesis approach. 
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1 Introduction 

It IS generally accepted that microgravity space experiments will need to be 
isolated from the vibrations inherent on spacecraft in earth orbit[3]. The funda- 
mental constraint on any isolation system’s capability is the available working 
envelope[4]. Figure 1 shows the relationship between the envelope (peak-to-peak 
displacement) and frequency for several sustainable RMS acceleration levels. 

The graph is for a one degree-of-freedom case and assumes sinusoidal vibra- 
tions, but the relationships are acceptable for order of magnitude estimates 
even if these assumptions are relaxed. 

No definitive specification of the required isolation levels or frequency range 
exists. The proposed US Space Station usable specification^] is also shown in 
Figure 1. It is claimed that vibrations below this curve will not adversely affect 
microgravity experiments. We have pursued the design of an active isolation 
system with a ‘reasonable’ envelope of 4 inches of travel, and a sustained 1 fig 
RMS acceleration. It can be seen from the figure that this will offer isolation 
down to 0.002 Hz. The amplitude to which vibrations can be attenuated is con- 
strained only by controller design and available instrumentation. Operation at 
lower frequencies, however requires a larger envelope, which becomes prohibitive 
in terms of available spacecraft space. We have also required that the system 
be active in all six degrees-of-freedom, with a rotational range of 40 degrees. 

Redundant coarse-fine schemes with magnetic levitation for vibration iso- 
lation are discussed in the robotics literature[2]. This approach is particularly 
attractive in the microgravity application since it allows the use of magnetic 
levitation while overcoming range of motion limitations. We have chosen the 
Stew r art platform for our coarse stage and a novel magnetic bearing for the 
fine stage, ^he approximate regions of activity in the frequency-displacement 
plane of these two devices are shown in the figure. Both stages act to attenu- 
ate spacecraft vibrations, effectively reducing vibration amplitudes below their 
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Figure 1: Peak- to- Peak Displacement vs. Frequency for Various RMS Acceler- 
ations, U. S. Space Station Usable Specification, and Activity Regions of the 
Two Actuator Stages 


active regions on the plane. As an example, it can be seen in the figure that 
a vibration of the spacecraft with 10 inches of displacement at a frequency of 

1 Hz falls outside the active region and could only be partially attenuated. It 
should be noted that such a large vibration is unlikely. If the displacement was 
only 1 inch, however, the coarse stage would absorb all of it except about 0.005 
inches, and the remainder would be reduced down to the micro-g level by the 
fine stage. 

The combination of the Stewart platform and a magnetic bearing allows 
continuous isolation at frequencies above 0.002 Hz, and a compact, reliable 
package suitable for the application. These choices and some preliminary design 
concepts will be discussed in detail. 

2 Stewart Platform 

The Stewari platform is a six degree-of-freedom parallel manipulator first pro- 
posed by Stewart[5]. It has been extensively used in aircraft cockpit simulator 
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applications, and substantial design information is available in the literaturefl]. 
Figure 2 shows the mechanism in our proposed configuration. Six linear. ac- 
tuators (legs) connect a base (bottom) to a platform (top). The base will be 
mounted in the spacecraft and move with it, while the platform tracks an inertial 
reference frame. VVe propose the use of stepper motor driven ball lead-screws 
as actuators. 



This mechanism was chosen over other candidates such as a carriage/gimbal 
approach, or a serial linkage mechanism because it has the following features. 

• Inherent rigidity. The parallel connection of the actuators gives the mech- 
anism rigidity on the order of the extensional rigidity of the actuators. 
For the proposed actuators, this will allow controller design to ignore the 
dvnamics of the mechanism. The effects of ‘umbilical connection to the 
platform will also be negligible. 

• Determinate inverse kinematics. The actuator lengths required to achieve 
a prescribed orientation are found directly from a coordinate transforma- 
tion from the base to the platform frame. This is seldom the case for a 
serial linkage. This will also simplify control. 

• Compactness. The configuration proposed here places the fine stage on 
top of the platform for convenience in testing. A fully developed imple- 
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mentation could locate the fine system and microgravity experiment in 
the space between the base and platform, resulting in a compact package. 

The Stewart platform has some disadvantages that must be considered. It is 
nonlinear in its response to actuator lengths, its general direct kinematics have 
not been discovered in closed form, and it has singularities in its operational 
space. The fist two problems can be overcome with digital controls. The singu- 
larities, which are points or loci where the mechanism gains a degree of freedom 
and the actuators can lose control of it, must be addressed in design. 

A simulation code has been written to allow exploration the design alterna- 
tives. Figure 2 is an example of its output. Preliminary results indicate that our 
specification (4 inches translation, 40 degrees rotation) will be achievable with 
actuators 10.5 inches long in the retracted position, and with 9 inches of stroke. 
The simulation can confirm that singularities are safely outside the working en- 
velope. Commercial actuators with the required range, load capacity, speed and 
acceleration have been identified. 


3 Magnetic Bearing 

Two axes of a six axis magnetic bearing are shown in Figure 3, mounted atop 
the Stewart platform. A ferromagnetic cube is at the center of the bearing. 
Two pole pieces protrude from each of its faces (four shown) and each pole 
piece is surrounded by a coil. This part of the structure is called the core and is 
mounted to the platform with four posts. Three ferromagnetic bands surround 
the core (one shown) forming three independent magnetic flux paths. The core 
is capable of exerting three orthogonal forces, and three orthogonal torques on 
the bands. For the axes shown, equal currents in each pair of adjacent coils will 
cause magnetic flux to flow in a local circuit, causing an attractive force to the 
band. By controlling these currents a prescribed force can be exerted on the 
band along the axis that crosses the page form left to right. If the currents in 
adjacent coils are not equal, some flux will flow around the outside of the band 
and through the center of the cube. This will create a controllable torque on 
the band around the vertical axis. 

Similar pole pieces and coils will protrude from the other faces of the cube, 
and corresponding bands will surround them. These have been omitted so that 
all parts can be seen. Also, the size of the bearing and the gaps have been 
exaggerated for clarity. Flux sensors will be mounted in the pole pieces and 
this will allow the position of the bands relative to the core to be calculated 
for control. The microgravity experiment will occupy the space surrounding the 

bearing, and be attached to the bands. 

This configuration was chosen over other levitation approaches such as Lor- 
entz actuators or magnetic actuators located on the periphery of the experiment 
package because it has the following advantages: 


4 


Figure 3: Magnetic Bearing 


. Compactness. The high force capability of the magnetic bearing ■ relat1 ^ 
to a Lorentz actuator of similar size and power consumption suits the 
application. Testing in earth gravity will be facilitated, and levi a ion 
during launch to protect sensitive instrumentation be ^ 

the rigid structure required to mount actuators around the periphery 

avoided. 

. Force/torque balance and rotational range. Actuators capable of the re- 
quired forces mounted on the periphery of the exp 
torques far greater than is required, and they 

of the experiment. The proposed design approach brings the rela 
force/torque magnitudes closer to the requirement, and allows subst 
tial rotational range. 

. Integral sensor capability. Compact semiconductor magnetic flux sensors 
(hall effect or magneto- resistive) can be utilized to both s ... 

[em and infer relative position. No elegant integrated approach ,s known 

for Lorentz actuators. 

Magnetic bearings have typically been avoided in 'large gap’ applications be- 
cause of their nonlinearity (force is proportional to the squar of flux . We 
feel that emerging Digital Signal Processor technology and contro 
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Figure 4 : Coarse-Fine Actuator 
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1 Introduction 

The authors have previously presented a conceptual design for a coarse-fine actuator pair and discussed 
its efficacy in the microgravity vibration isolation appiication[l]. The coarse stage comprises a Stewart 
platform [2] which is mounted in a spacecraft and isolates low frequency, high amplitude vibrations. The 
fine stage is a novel magnetic bearing mounted on the Stewart platform (between the legs for compactness) 
and levitates the experiment to isolate all frequencies at low amplitudes. The combination is illustrated in 
Figure 1. 




Figure 1: Coarse-Fine Actuator Pair 

This paper will present a survey of published 6 DOF levitation designs and discuss a novel magnetic 
bearing in terms of design, predicted performance, and control issues. 

Supported in part by the NASA Lewis Research Center and the Commonwealth of Virginia’s Center for Innovative 
Technology. 
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2 Survey of Published Designs 

Several designs for 6 DOF levitation are discussed in the literature. A comparison of the specifications for 
these designs is given in Table 1. 


Group 

Trans . 

Rot. 

Force 

Envelope 

Weight 

Actuator 

Sensor 

Honeywell 
N. Wales 
NASA 
SatCon 
IBM 
Toshiba 

±5 mm 
±5 mm 
±4 mm 
±10 mm 
±5 mm 
±2 mm 

±1.6° 
< ±.2° 
±3° c 
±8° c 
±4° 
±1.5° 

43 N 
.04 N “ 
445 N 
4 N 
32 N 
20 N c 

27x34x50 cm 
100x100x100 cm b 
30x30x15 cm c 
40x40x12 cm c 
25x25x15 cm c 
25x25x20 cm 

36 kg 
? 

? 

4.9 kg 
? 

8 kg 

Mag. Brng. 

Lorentz 
Mag. Brng. 
Lorentz 
Lorentz 
Mag. Brng. 

Eddy & Flux 
Capacitive 
Eddy 
Eddy 
Optical 
Eddy 


Requirement, not limitation 
includes experiment package 
c Estimated by authors 


Table 1: Comparison of Published Designs 

Four designs specifically for microgravity isolation have been published. Honeywell [3] has a well devel- 
oped system called FEAMIS with which they have demonstrated impressive isolation performance. The 
system is designed for the Space Shuttle experiment configuration. The University College of North Wales 
[4] also has a well developed system designed for the European Space Agency experiment configuration. 
NASA [5] and SatCon [6] both have laboratory levitation systems. 

Two levitation designs were developed for different applications, but they are mentioned here because 
they are similar and could be easily adapted to the isolation application. IBM [7] has a laboratory levitated 
robot “wrist” which enhances robot accuracy and performance. Toshiba [8] has a satellite antenna pointing 
system which is fully developed. Both devices have demonstrated positional accuracies on the order of 1 pm. 

Isolation of vibrations with large amplitudes — typically occurring at low frequencies — requires a 
large translational range. SatCon’s system has the largest range, but there is a significant tradeoff with the 
device’s force capability. A coarse-fine approach would allow both a large range, provided by the coarse 
stage, and a high force capability, since the levitation gaps are small. There is no available data on the 
rotational range requirements of the application. Isolation with an umbilical disturbance requires a high 
force capability as is offered by the systems from Honeywell, NASA, IBM, and Toshiba. Space and weight 
should be minimized in any spacecraft. SatCon, IBM, and Toshiba’s systems offer advantages in envelope 
space and weight. 

The choice of the actuator technology between Lorentz force and magnetic bearings has no definitive 
advantage. Lorentz actuators offer linearity, simplicity, and compactness. Magnetic bearings offer higher 
force capability and lower power consumption, particularly if gaps are minimized. 

Four position sensor technologies offer promising performance. Eddy current position probes are simple 
and robust, but bulky and heavy for large gaps. Capacitive sensors are simple and light weight, but can 
be noisy in unconstrained environments. Optical lateral effect photo-diodes are compact and quiet, but 
they require substantial supporting electronics. Hall effect flux sensors can be used with magnetic bearing 
designs both to linearize the control problem, and to measure position. 

3 Design 

The magnetic bearing proposed has two parts: a stator which is attached to the spacecraft, and a sur- 
rounding “flotor” to which the experiment is attached. 



The stator is illustrated in Figure 2. It has twelve pole pieces and coils arranged around the surface 
of a cube. The cube and pole pieces are ferromagnetic. Each pair of pole pieces and the region of the 
cube to which they are attached comprise a typical “horseshoe” electromagnet causing an attractive force 
toward the nearby flotor. Magnetic flux through the center of the cube will cause an imbalance in the flux 
levels of a pair of pole pieces, resulting in a net torque on the flotor. Differential Hall effect sensors will be 
located in the cube side of each pole piece to measure the local flux. All electrical connections wiU be to 

the stator. 



The flotor is illustrated in Figure 3. Three ferromagnetic bands are rigidly attached to each other, but 
form independent flux paths. The bands are thicker in the region near the pole pieces to avoid saturation. 
Flux which passes through the center of the cube is returned through the remaining portion of the bands. 

Four mounting posts will attach to corners of the cube, and pass through clearances in the flotor. These 
posts could carry cooling fluid to be circulated through the stator if it is required. 

Design equations relating force and moments to the coil currents will be derived below referring to 
Figure 4 The figure shows a schematic slice through the stator and flotor with appropriate nomenclature 
and sign convention information. It should be noted that a complete model comprises three such systems, 
but they are identical and orthogonal, so only one will be analyzed. 

The relationship between coil currents (i, .,) and the force and moment generated m one slice 

of the stator (F„M.) can be derived from Maxwell's Equations. The first Maxwell equation , ( 1) .which 
relates magnetic field intensity (H) around a closed path to the electric current density (J) through that 
path, is discretized and applied to closed loops drawn through the slice. A is the number of turns in each 
coil, and Gi are the air gap lengths which are dependent on the stator s position relative to the flotor. 
The iron flux paths are ignored because their reluctance is low relative lo that of the air gaps. Many such 
equations can be written (2), but only three are independent. 


IJ (U = da 


( 1 ) 
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Figure 3 : Flotor 
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( 2 ) 


The second Maxwell equation ( 3 ), which ensures conservation of magnetic induction (B), is used to 
obtain a fourth independent equation ( 4 ). 


V • B = 0 ( 3 ) 

B x + B 2 — B 2 — B4 = 0 ( 4 ) 

We can assume linear magnetization in the air gaps ( 5 ), where /to is the permeability of free space, to 
obtain a relation between magnetic induction in the gaps and coil currents ( 6 ). 
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Figure 4: Schematic Cross Section of Magnetic Bearing 
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( 6 ) 


The four gaps are geometrically related to the offset of the stator with respect to the flotor by the 
relations (7) which assume small angles. 


G \ = Go — y — sB 

G 2 — Go — y + s 9 ( 7 j 

G*3 = Go + y — $8 
G 4 = Go + y + $0 


Go is the air gap length with the stator centered in the flotor. 
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The magnetic energy stored in the magnetic bearing (u m ) is found from (8) where A g is the area of 
the pole faces. 




■ \L*-* dv 

— 2 ^ + B^G? + B3G3 + BlG^j 


(8) 


The force and moment on the stator are found from the relations (9) and (10). 

F - dUm 

Fy ~ TjT 

if & 

M '"sr 


(9) 


( 10 ) 


After considerable algebraic manipulation, and introduction (without loss of generality) of the linear 
current transformations (11) we obtain the force and moment relations sought (12) and (13). 


h = t't - *2 + »3 - *4 (11) 

J3 = i\ — »2 — *3 + U 


yvVo (Gq] 2 + Osji + yjj) (Glj 3 - 8 2 s 2 j 3 + Gpyjg + Qsyji) 
4G o (G 2 -0V-y 2 ) 2 

AgNhs ^ o (G 0 ; 2 + 0*j\ + yhHGpji + G o 0sj 2 + 0syj 3 - y 2 j x ) 

4G 0 (G 2 -^ 2 -y 2 ) 2 


The current j 2 is analogous to the bias current in a conventional bidirectional thrust bearing and 
could be fixed at a constant value — nominally half of the maximum current. The force generated is 
predominantly driven by j‘3 and the moment by j\. The system is unstable (negative stiffness) in both 
translation and rotation. The currents i 1, . . 14 can be found by a pseudo-inverse technique from Ji, J 2 ij 3 * 
Closed form analytic inverses to ( 12 ) and (13) have been found for a known position. 


4 Predicted Performance 

The equations of the previous section were used to predict the performance of a specific design. The design 
has a center cube of 2 in. on a side, pole faces of 1 x .5 in., and pole length of 2 in. Maximum current 
is determined by allowing a coil current density of 5000 amp/in 2 which is known to be conservative from 
previous designs. The gap in the centered position was chosen to be .125 in. plus an allowance of .030 in. 
for inclusion of flux sensors and a protective layer on the inside of the bands. The resulting specifications 
for the design are presented in Table 2. The 53 N force is a continuous worst case, with the stator moved 
away from the flotor in the direction of the force. The continuous force capability in the centered position 
is 175 N. Intermittent force capability is limited only by the current capability of the amplifiers, and the 
saturation Limit of the magnetic material used. Using Vanadium Permadur with this design, saturation 
would occur at about 1000 N. Of the 1.5 kg weight, the flotar comprises only 1.2 kg. 


G 



Trans. 

Rot. 

Force 

Envelope 

Weight 

±3.2 mm 

±7° 

53 N 

15x15x15 cm 

4.5 kg 


Table 2: Specification of UVA Design 


When compared with the designs presented in Table 1, the UVA design has several advantages. The 
envelope is substantially smaller than any of the previous designs, while the performance is similar. In 
addition to saving space, this compactness allows the flotor to be naturally rigid, and thus avoids control 
problems with structural dynamics. The design is quite dense in comparison with the others, but it is 
lighter than the lightest for which data were available. 


5 Control 

A regulator has been designed to reject the disturbances caused by the umbilical connection to the exper- 
iment. A schematic is shown in Figure 5. Nonlinearities in the magnetic bearing are eliminated by using 
flux feedback in a minor loop [3]. Six accelerometers mounted on the flotor produce a generalized accelera- 
tion signal which is fed back through a linear controller. More details on the controller are available in [9]. 
The desired control force is processed through an inverse magnetic circuit model to obtain a desired flux 
signal. This model could be either a digital algebraic model, or an appropriately trained neural network. 



Figure 5: Control Schematic 

A relative position sensor has not been chosen but the optical scheme used in IBM’s design is a strong 
candidate. Alternatively, the current and flux signals could be processed to infer position [10]. The purpose 
of the relative position signal is only to prevent collision with the walls, so accuracy demands are relatively 
low. 
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6 Conclusion 


A design for a novel magnetic bearing, proposed as the fine stage of a coarse-fine actuator for microgravity 
vibration isolation, has been presented. The bearing is novel in that it uses a geometry that has just 
three independent flux path systems. This contrasts the twelve flux path systems (six bidirectional thrust 
bearings) used in conventional designs. The novel design results in compactness, light weight and high per- 
formance, when compared with the published designs. A control system is proposed to reject disturbances 
caused by an umbilical connection to the experiment. 

Future work will focus on building a laboratory version of the bearing and control system. 
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Nomenclature 

= matrix in state dynamics 
= vector in state dynamics 
= constants in optimal control 
= microscopic component of y(t) 

= cumulative distribution function 
= probability density function 
= final condition cost 
= transfer function of shaping filter 
= Hamiltonian 
= cost function 

= maximum stroke of experiment 
= fourth moment of a sample of T~ x 
= bandlimited white noise 
= power spectral density. 

= complex frequency variable 
- first-passage time 
= time 

~ tj 1116 at significant maximum 
= time at zero crossing 
= experiment acceleration 
= state vector, = (.*■, x 2 ) r 
= experiment position 
= experiment velocity 
= wall process 
= upper wall true position 
~ lower wall true position 
= upper wall constraint 
= lower wall constraint 
— maximum of wall process 
= macroscopic component at y(t) 

= macroscopic component velocity 
= significant maximum position 
= dirac delta function 
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X, 

- Lagrange multiplier for dynamic equations 

<t> 

= Lagrange multiplier for final condition 

V 

= expected frequency of uperossings (cycles) 

V 

- expected frequency of maxima 

O) 

= frequency, rad/s 

Ud 

= damped natural frequency, rad/s 

£ 

= damping coefficient 

Subscripts 

d 

- microscopic component 

f 

= final time 

max 

= at maximum 

N 

= normal 

n 

= bandlimited white noise 

T 

= first-passage time 

y 

= wall process 

y 

= macroscopic component 

z 

= significant maxima position 


Introduction 

T* FIE microgravity environment of space may permit 

X advances in material science experiments. Such experi- 
ments could aid in the understanding of basic physical phe- 
nomena, quantify the limitations and effects imposed by 
gravity, and spur application to Earth and space based 
processes and products. A microgravity environment could 
potentially eliminate buoyancy-driven convection, sedimen- 
tation and hydrostatic pressure as well as yield other advan- 
tages.' 

At this time, the acceleration environment requirements for 
various experiments are not well known. 2 An assessment of 
existing theoretical and experimental data available up to 1985 
indicated acceleration levels below 10 6 So would be required 
lor frequencies below 0.1 Hz for many of the processes. The 
requirements at higher frequencies are somewhat lower 3 
Work to determine the levels necessary is in progress by 
various materials experiment researchers. An example is a 
twin crystal growth experiment to be carried out on space- 
craft. 4 v 

An essential part of the development of a microgravity 
experiment program is the characterization of the low 
acceleration environment aboard spacecraft. NASA has 
carried out a series of measurements reported at various 
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Table 1 

g'go 


10- 7 

10-8 

10- 7 


2x IO' 2 
2 x 10- 3 
2 x 10" 4 


IO ' 4 

IO' 4 


Mierogravity space experiment acceleration environment 7 


F, Hz 

Quasisteady or “dc” 
0 to 10~3 
Oto IO' 3 
0 to 10-3 


____ Source 

acceleration disturba nces ~~~ 

Aerodynamic drag 
Light pressure 
Gravity gradient 


Period acceleration disturba nces 


Thruster fire (orbital) 
Crew motion 
Ku band antenna 


Nonperiodic acceleration disturbances 


Thruster fire (attitude) 
Crew push off 


EXPERIMENT ENCLOSURE 


EXPERIMENT 



UMBILICAL ACTUATOR 


conferences. 5,6 A summary of this data was presented in Ref. 7 
and is repeated in Table I. Additional results have been re- 
ported in Refs. 8 and 9. 

The vibration levels reported in the above literature for 
spacecraft are significantly greater than allowable for material 
science experiments. In order to achieve accurate and repro- 
ducible results in such experiments, vibration isolation will be 
required. 3 Acceleration disturbances in the orbiter environ- 
ment cover a wide frequency bandwidth, from 0 to 100 Hz. 
Sources include spacecraft drag, light pressure oscillations, 
manned activity, and thruster fire. The frequency and 
amplitudes of these accelerations are summarized in Table 1. 

The frequency and amplitude of any particular vibration 
source determines the level of isolation that can be achieved. 
At relatively high frequencies, above approximately 10 Hz, 
passive vibration isolation is normally possible. Two examples 
of such isolation systems are reported in Refs. 10 and 1 1. For 
lower frequencies, active vibration isolation is necessary. One 
of the few such systems is examined in Ref. 12. 

. ^ fundamental restriction on active microgravity vibration 
isolation systems is the limited available volume aboard space- 
crafts for experiments. This kinematic constraint cannot be 
overcome through improvement of the vibration isolation 
control system, sensors, or actuator. The purpose of this 
paper is to explore the limitations on vibration isolation sys- 
tems arising from the stroke restriction. Thus an ideal vibra- 
tion actuator is assumed, and the problem is solved in part by 
optimal control theory. This work is an extension of research 
on isolation limits under sinusoidal excitation. 13 In this paper, 
the excitation is a wideband, zero-mean stochastic process. A 
lower bound on the root-mean-square acceleration is deter- 
mined in terms of the maximum stroke. 

Optimal Control Formulation 

While the isolation problem for microgravity space experi- 
ments is multidimensional, this analysis examines the one-di- 
mensional case. Consider the system illustrated in Fig. 1 with 
experiment position x(t) and wall positions y { (t) and y 2 (t). 
The experiment is connected to the spacecraft by umbilicals, 
such as power or fluid lines, and by a vibration isolation 
actuator. Although the spacecraft has a Finite mass, it may be 
considered to have infinite impedance for this analysis since 
the spacecraft-to-experiment weight ratio is very large. The 
spacecraft acts as an external base motion transmitting forces 
through the umbilical and the actuator. 

An ideal actuator is assumed. Therefore, the acceleration of 
the experiment is the minimum acceleration possible given that 
the experiment stay between the two walls. The effects of the 
power/data/cooling umbilicals represented in Fig. 1 are re- 
moved through the ideal actuator. Thus, the problem is re- 
duced to a kinematic representation. The vibration isolation 
problem becomes an optimal control problem: find the opti- 
mal trajectory (minimum acceleration) given the constraint 
conditions (moving walls). This problem formulation was 
used previously by the authors to find the limitations on 
isolation for harmonic disturbances. 13 


Fig. 1 One-degree-of-freedom isolation problem. 



Fig. 2 Optimal path through stochastic walls. 


Significant Extrema 

The kinematic formulation allows the stochastic problem to 
be easily conceived. Figure 2 illustrates. An optimal trajectory 
x(/) is sought between two walls whose motion is described by 
a single zero-mean stochastic process y(t), 

=>(/) + l n 
.Flow = y (0 L> /2 

where L is the maximum stroke of the experiment between the 
two walls. 

Let the experiment acceleration be denoted as u(u = x). 
Then the cost function J to be minimized is 

(i) 

with the constraint 

, „ L L 

y(t)--%x(t)<y(t) + -, 0 (2) 

for a g i ven wall centerline motion y(t). The optimal trajectory 
will be in general smooth with as few extrema as possible. The 
trajectory will also cover as little distance as possible As pre- 
viously reported in Ref. 13, with a sinusoidal disturbance the 
optimal trajectory tangents the maxima of the lower wall and 
the minima of the upper wall when the amplitude of y(i) is 
greater than L/2. Given these tangencies, it is an easy task to 
, compute the minimum acceleration path connecting them Al- 
though this path may intersect a wall’s path, as was reported 
in Ref. 13, this “cheating” results in less than a 1% reduction 
in root-mean-square from the true optimal path. By dropping 
the inequality constraints from the optimal control problem 
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Solution in Terms of the Significant Extrema 
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The velocity boundary condition at time t 2 is left free since 
this condition depends on the next significant extremum the 
tocation of which cannot be readily derived even in a prob- 
abihstic sense. Therefore, the velocity at time t 2 is unspecified 
so that the return path root-mean-square acceleration will be a 
lower bound on the minimum acceleration. The optimal con- 
trol problem is thus reduced to that over the time /, to /, with 
cost function 2 


frequency 17 


E(y 2 ) = a] = 


Sy{w) do? 


( 8 ) 


T S h7m y Ji e , area und f «neath the power spectral density curve) 
The mean square of the velocity and position processes are 


fh 


/ = 


u l dt 


and the boundary conditions of Eq. (3). 

This problem is solved by the calculus of 
(Appendix A) yielding the optimal trajectory. 14 
The mean-square acceleration on this trajectory is 


(4) 

£(y 2 ) = o 2 = 

LV- 

(9) 

variations 

E(y 2 ) = a) = 

f 

J -00 0) 

(10) 


MS ( x) = f 7~~— f'\ 2 df]=^ 

Lte-f.) J,, J r 4 


(5) 


cr n r ‘ he tlme from s, 8 nl ficant extrema to zero 

south? 8 * ° the fi^t-passage time. This is the lower bound 
sought. Assuming that the wall process y(t) is ergodic the 

wthT^ | SqU3re ex P eriment acceleration attainable 

acreler?t!ot a tf S r°h at0r mUS ‘ ** greater than the me an-square 
acceleration of the composite trajectory, 
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where E[ ■ J is the exact value operator 15 


Y* fa (*./) d z dr 


( 6 ) 


> ) is the joint probability density function on the 

efo^o? irrtTh 3 ^- ^ S ' g " ifica "‘ “a * 

rtl A r g ' • (The s'Snificant extrema problem can be 
%e or ln terms of significant maxima because of the 
symmetry of the problem.) Equation (6) can be rewritten 
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dz 
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* h ? r 5 Fl S X/ T*)\z I is the conditional fourth moment of T~ l 
maxima ** ^ probabllity densit y function of the significant 


Distribution of Significant Maxima 

nr £®" Sl ^ ef tbe stationary, zero-mean, Gaussian random 
process describing the wall’s acceleration j >(/). This process 
can be characterized by its power spectral density S/u).' 6 The 
mean square of the wall acceleration is the integral of this over 
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fs(x) = (2t)~ l/i e - (Jr 2 /2) 


1 = 


fa(x) = f(a) da 


(14) 


bmTon.) 35 a " imerCSting di scussion on this maxima distri- 

The ratio ^/ Vy can be regarded as the average number of 
maxima between uperossings for the y(t) process 18 Prob 
ability density functions for the distribution of maxima of^' 
stationary Gaussian random process are shown for different 
values of this ratio in Fig. 6. uuierent 

A random wall process typical of the microgravity vibration 
environment will have many maxima pj zero cro« n? 
(Hy/v y > 1 ). Most of these maxima are clustered near the 
peaks of the process and can be considered the high frequency 
component.' 8 Figure 7 illustrates. The random process yU) 
as was shown by Crandall,' 8 can be decomposed into the sums 

ponenTrf( r ,) SCOP1C COmP ° nem 3nd 3 micros «>pic com- 


y(t)=HO + dt 


(15) 


denSi,y funf ' ion for ,hf distribution of maxima of 
a stationary Gaussian process. 


The macroscopic component has a much smaller root mean 
square a, and greater frequency of zero crossing „ ZZTc 

ZT C T S ,: gna ' ( °>' ^ The macroscopic component will 
have nearly the same amplitude as the original signal The 

bmh ,hT ° f maXima 6iVCn in Eq ' (13) is the distribution of 
both the micromaxima and macromaxima. Ciearlv the 

micromaxima do not contribute (substantially) to the cost 



function; therefore, they are, by definition, not significant 
maxima. 

The distribution of significant maxima is therefore derived 
from the distribution of macromaxima, the maxima of the .HO 
process, 


/w a >/*,) = <i -«*/„])«/ 


y/py 
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To obtain the density function for the distribution of 
significant maxima from the distribution of macromaxima it is 
only necessary to recall Eq. (3) with y m3LX replacing y miX 


r p 


z = 


t -4/2 


J > max>4/2 
^ max <^4/2 


(17) 


find the cumulative distribution on z in terms of the cumu- 
lative distribution on y mAX 


of the geometric configuration, Fig. 5, in representing typical 
wall position histories. The decomposition, or smoothing, of 
the y(t) process yields a macroscopic random process HO 
upon which to base the lower bound. If the random process 
HO has close to one maxima per upcrossing on average (it 
must have at least one since it is continuous), the typical 
geometric configuration will resemble Fig. 5. Thus, for a 
narrow-band random process y(t ), the bound will be very tight. 
(In the limit this approaches a sinusoid with slowly changing 
random amplitude and phase 22 ; the resulting bound would be 
only a few percent low.) For a wideband macroscopic process 
the bound will be considerably looser but non-trivial. It is 
therefore desirable that the decomposition produce a macro- 
scopic process with an average number of maxima per up- 
crossing (jiy/vy) close to one. The smoothed process HO 
should also have nearly the same variance as the original signal 
HO- A method for the decomposition of HO and d(t) will be 
discussed later with the application. 

First-Passage Problem 

An expression for E 1(1 /T A )\z ) can be found explicitly if 
the conditional probability density function f T (T\z) is 
known 15 


F z (z) = 



z >0 
z <0 


and to differentiate with respect to z 15 


(18) 


/x(0» 


\r, tm (L/2Mz)+/ fm Jt*+L/2) 

0 


Z >0 

(19) 

z <0 


where 5( - ) is the dirac delta function. Note that the portion of 
the Him distribution between - oo and 4/2 maps into a dirac 
delta at z = 0. The integration of this part of f z {z) in Eq. (7) 
will be zero 


3 z 2 Fy mix (£/ 2) 8(z) E 




dz = 0 


( 20 ) 


The nonsignificant maxima, therefore, do not increase the 
expected mean square. (Note that E[(\/T 4 )\z = 0) is finite 
since T is defined as the time for the HO process to reach y - 0 
starting from rest at y mix . When z = 0, y m3LX = 4/2; therefore, 
T is greater than zero.) 

The tightness of the lower bound on experiment mean 
square acceleration given in Eq. (7) depends on the accuracy 



Fig. 7 Macroscopic maxims and microscopic maxima. 


E\(\/T*)\z\ = fr(T\z) dT (21) 

Jo 1 

The distribution of the first passage time of a stochastic 
process T is a classical problem probability that has been 
tackled with limited success by many authors. 23-25 No exact 
theoretical solution exists. Many approximate methods have 
been used including series solution, 23 Poison approximation, 24 
numerical probability diffusion, and Monte Carlo tech- 
niques. 25 The particular variant of the first-passage problem 
of interest here is nonstationary due to its initial conditions 


y(t\) — y’mix 


y('i) = 0 (22) 

and can be expected to have a crossing in the near future since 
HO is zero mean (This, in contrast to the first passage of a 
level y - a > a y , which will occur extremely infrequently.) The 
nonstationary, quick-crossing nature of this problem suggests 
Monte Carlo simulation as the preferred method of solution. 
It is not necessary to generate the conditional probability 
density function of Tin this manner; only the E ( (1/T 4 ) \z ) is 
needed. 

To perform the simulation, a forced differential equation 
with the proper probabilistic characteristics must be found. 
More specifically, the wall acceleration in the simulation 
should have a power spectral density that matches that of the 
acceleration disturbance found aboard the spacecraft in the 
region where the experiment is to be located. Given a mea- 
sured power spectral density, the dynamic model, or shaping 
filter, that transforms Gaussian white noise to the desired 
random process can be found by spectral factor- 
ization. 17 - 26 This dynamic model can then be used in Monte 
Carlo simulation (Appendix B). 


Application to Microgravity Vibration Isolation 

The method described above is now applied to the 
microgravity problem. As described in Ref. 7, vibration 
aboard spacecraft is essentially of two frequency regions: 
0-10 3 Hz and above 1 Hz. Consequently, the power spectral 
density for analysis here is modeled independently in these two 
regions. The low frequency portion was not obtained from an 
experimental power spectral density since such information is 
unavailable. This part of the spectral density, illustrated in 
Fig. 8, has therefore been devised for the purpose of this 
example to fit the known environmental disturbances in this 
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Fig. 8 Acceleration power spectral density in the low frequency 
range. 

o 



Fig. 9 Acceleration power spectral density in the high frequency 
range. 


frequency range. The spectral density of the analytical model 


and by Eqs. (! 1) and (12) 

v y - 0.18836582 x " 3 upcrossings/s 

H y = 0.13153971 x 10' maxima/s 

Hy/vy = 6983.2 maxima/upcrossing (25) 

Clearly there are many micromaxima per peak of the macro- 
scopic process. The two distinct regions of the power spectral 
density permit an easy decomposition of the wall processes 
y(t), XO, and 'y(t) into macroscopic (low frequency) and 
microscopic (high frequency) components. For a problem with 
a continuous power spectral density, a smoothing filter as 
examined by Crandall 18 may be used for the decomposition. 

The decomposition of the Gaussian signal yields two inde- 
pendent random processes 21 

E\mdU) 1=0 a 2 y =al + o 2 d 

E\7(t)d(t)\ =0 o*-o* + ol 

E[7(t)d(t)) = 0 <7? = a 2 + a 2 d (26) 

with 


o y = 0.16503623 x 10' m 
<ty = 0.19030049 x 10 “ 2 m/s 
<*y = 0.25429066 x 10" s m/s 2 

The resulting upcrossing and 
v y = 0.18351888 x lO ’ 3 
Hy = 0.21267216 x 10 “ 3 
Hy/vy = 0.1 1588571 x 10' 


o d = 0.63558716 x lO ’ 4 m 

o d — 0.44024768 x 10 ~ 3 m/s 

o d = 0. 1614351 1 x 10 " 1 m/s 2 
(27) 

maxima frequencies are 
v d = 0.11024073 x 10 1 
= 0.58360775 x lO ’ 1 
Hd/va = 5.2939393 (28) 


Sy{<J) 


75.75 

(u/0.001 16) 4 - 1.99 (y/0.001 16) 2 + 1 

0 


0< <0.006283 


0.006283 <03 <3. 14159 


(23) 


has most of its power concentrated near orbital frequency 
(1.851 x 10 ” 4 Hz, 90 min orbit) with a power corresponding 
to a root-mean-square acceleration of approximately 0.2 pg. 
It should be pointed out that the accuracy of the lower bounds 
generated by this method depends heavily upon the model 
used in this region. 

The high frequency region of the power spectral density, 
shown in Fig. 9, is from experiments aboard Spacelab . 33 An 
analytical model could easily be fitted to the curve shown; 
however, this is not necessary for this application. Note that 
the method of analysis described previously does not require 
that disturbance spectral density occupy two separate regions 
as in this case. 

The mean square wall process acceleration, velocity, and 
position can be found by evaluating Eqs. (8-10) respectively. 
This may be done through numerical quadrature yielding, 

a, = 0.016143511 m/s 2 

Oy — 0.0019532654 m/s 

o y ~ 1.6503623 m (24) 


Therefore, the smoothed wall process has an average of 
1.158857 maxima per upcrossing (or cycle). Since this is dose 
to one, the wall process geometry closely resembles Fig. 5 and 
the lower bound will be tight. Further note that 

Oy/Oy = 1.0 v y/ ¥ y — 0.97426848 (29) 

Thus, the smoothed process retains nearly all of the signal’s 
amplitude and has on average fewer cycles per unit time. The 
apparent drop in “frequency” is due to wall process 
smoothing removing spurious microcycles. Therefore, the 
decomposition does not alter the original signal’s significant 
maxima amplitude-time characteristics. Finally, unless the 
maximum stroke is of the same magnitude as the microscopic 
component, the smoothed wall acceleration J(t) establishes an 
upper bound oy on the minimum root-mean-squarc experi- 
ment acceleration. 

With the significant maxima distribution found, attention is 
now turned to finding the needed first-passage time moment. 
The low frequency part of the power spectral density, 
Eq. (23), is converted to an equivalent dynamical system via 
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spectral factorization yielding the shaping filter 

8.64581 


Cy(5) = 


(s/0.001 16) 2 + 0. 1(5/0.001 160) + 1 
with a bandlimited white noise excitation n(t) described by 
1 0 < cj < 0.006283 


(30) 


S n ( w) = 


0 0.006283 <o> 


(31) 


Converting Eq. (30) into two first order differential equations 
produces 


=J >2 (32a) 

7 2 = - 1.3456 x 10 - 6 j>! - 1.16 x 10- 4 ^ 2 

+ 1.16338 x lO" 5 n (32b) 

with initial conditions from Eqs. (15) and (22) 

7i(0) =7i 0 = ^ + L/2 

J> 2 (0) = 0 (33) 

The bandlimited white noise process n(t ) is approximated 

by the sum of a large number of sinusoids at nearly evenly 
spaced frequencies between zero and the cutoff frequency via 
the method advocated by Shinozuka and Jan. 27 The phases are 
random while the amplitudes are determined by the power 
spectral density desired. As the number of sinusoids is 
increased, the approximated signal approaches n(t). It was 
determined through testing that 40 was a sufficient number of 
sinusoids for this example. 

The differential equations of Eq. (32) are numerically 
integrated from the initial conditions until the first crossing. If 
the process rises above the initial position the simulation is 
stopped and restarted from the initial condition again. This 
insures that the crossing time obtained is the first-passage time 
from the maxima. The first-passage time of each simulation is 
recorded, and the fourth moment of T~ 1 is calculated from 
the collection of k first-passage times 

"»-„(.p, 0 ) = l ET,-* (34) 

K l m l 


o 



Fig. 10 Numerically computed fourth moments of T~ l . 
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Fig. 11 Lower bound on root-mean-square experiment acceleration 
under prolonged exposure to specified vibration environment. 


For an ergodic process, 16 as the sample size k grows, this 
statistic approaches £(r~ 4 I.Vi 0 |. The statistic is 

found for a range of y ]Q values. The results of the simulations 
for the example are shown in Fig. 10. Note that the data are 
asymptotic to 


= 0.26105 x 10 “ 12 s - 4 


as 7i 0 approaches infinity where is the damped natural 
frequency of Eq. (32). This is because the noise is essentially 
negligible when the system’s energy is very high. Once the 
maximum stroke L is specified, the £ ( 7" ~ 4 17i 0 ) yields the 
E\T ~*\z\. The data points are cubic spline interpolated to 
provide an approximation of the function E(T~*\z) suitable 
for numerical integration. 

The experiment acceleration mean square lower bound 
sought can now be calculated from the significant maximum 
probability density and the conditional moment on the 
first-passage time by numerical quadrature of Eq. (7) using 
Eqs. (16), (19), and (28) along with the interpolation of the 
inverse moment data. This is done for varying values of 



maximum stroke. The lower bound on minimum root-mean- 
square experiment acceleration that can be obtained under 
prolonged exposure to the example vibration environment 
given the stroke constraint is thus computed and is plotted in 
Fig. 11 vs maximum stroke. For example, with a maximum 
stroke of 1 m, according to the figure the root-mean-square 
acceleration of an experiment must be greater than 0.17 n g t 
whereas if a 2-m stroke is permitted, this can be reduced to 
0.13 /ig. Note that the lower bound is non-zero for finite 
maximum stroke. This is in contrast to the previously reported 
results for harmonic disturbances where the maximum accel- 
eration was zero when the maximum stroke was greater than 
twice the harmonic amplitude. 


Conclusions 

In this paper, the microgravity vibration isolation problem 
was formulated as a one-dimensional kinematic problem. The 
geometry of the stochastic wall trajectories was defined in 
terms of their significant extrema. An optimal control solution 
for the minimum acceleration return path determined a lower 
bound on experiment mean-square acceleration. This bound 
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was expressed in terms of the probability density function on 
the significant maxima and the conditional fourth moment of 
the first passage time inverse. For an example given, the first 
of these was found analytically w hile the second was found via 
Monte Carlo simulation. The experiment root-mean-square 
acceleration lower bound as a function of available space was 
then determined through numerical quadrature. 

The method of analysis is quite general and intuitive. The 
authors feel that this could be applied to other problems with 
stochastic constraints. While such an analysis does not yield a 
controller, it does aid in the selection of system parameters 
(for example, the maximum stroke). 

Lower bound plots of the type developed here may assist 
microgravity experiment designers as well as vibration 
isolation engineers. The lower bound depicted in plots of this 
type could not be achieved by real systems for several reasons. 
The levels derived are beneath the theoretical minimums. In 
addition, a real system is causal and cannot base its current 
control on unknow n future disturbances. Also, any real active 
control system will have some non-ideal characteristics. The 
sensors employed to provide feedback will have some error as 
well. In spite of these comments, the authors of this paper are 
optimistic about attaining microgravity isolation levels close 
to the levels depicted here with an active vibration isolation 
system. 


and 
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where T - (t 2 - / ,), the time from significant extrema to zero 
crossing, or the first-passage time. 

Appendix B: Shaping Filter for Simulation 

The power spectral density is first approximated by a 
rational polynomial in w 2 . (This can be done to any required 
accuracy for a given power spectral density by using higher 
order polynomials.) This representation of the true power 
spectral density is then factored into two terms 


Sy(uj 2 ) = C(/o>)C( (Bl) 


Appendix A: Solution for Optimal Trajectory 

The derivation of the optimal trajectory proceeds as 
follows. Define the system state variables as 

Xi=x X\ = x 2 x 2 = u (Al) 

and the equations of motion become 


The first term G{jw) has all its poles and zeroes in the left half 
plane, whereas the poles and zeroes of the second term, 
G(-ju), are the mirror images in the right half plane 
Replacing yw in G(Juj) with the complex frequency variable s 
yields the transfer function that produces a Gaus- 
sian random process with power spectral density S(w) from 
Gaussian white noise input: 


x = Ax + Bu 

white noise 

0(s) 

J(S) 





x r = [x, x 2 ] 



1 

0 


B T = [ 0 1 ] 


(A2) 


The cost function is adjoined by the equations of motion and 
the final condition, utilizing three Lagrange multipliers. 14 The 
result is 

J~4>x lf =\ (u 2 + X,(jf 2 -3r 1 ) + X 2 (u-x- 2 )] dt (A3) 
J t i 


It is a straightforward process to convert the transfer function, 
or shaping filter, with n poles and m zeroes to n + m first 
order linear differentia] equations with white noise input via 
the inverse Laplace transform and algebraic manipulation. 17 26 
(Note that the derivative of white noise is nonexistent.) For 
most applications only the model of the low frequency 
component j>(r) as a set of differential equations is needed to 
generate the £ | (l/r J )lz j. 


which is the general functional for this problem. Define the . , 

Hamiltonian and final condition cost 14 as Acknowledgment 
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Employing the calculus of variations, the minimization equa- 
tions are 

\ m — - n \ dH ^ „ dH 

x ' a*r° o-^-2» + ma5) 

with natural boundary conditions 

X V = G *\j = *• X V = C X2/ = 0 (A6) 

Solving this gives 

= C 0 , X 2 = - C 0 t + C|, u = Vi C 0 t- Vi C, (A7) 

Imposing the prescribed and natural boundary conditions 
yields 

n 6z ^ 6 ^2 

c ° = yi ■ Ci = yr ( A «) 
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Nomenclature 

= matrix in state dynamics 
= vector in state dynamics 
= amplitude of base motion 
= coefficients of optimal control u Q t 
= force on experiment platform ° P ‘ 

= Hamiltonian 

— cost function 

= experiment mass 

= stroke (maximum translation of experiment) 
= half-period 
= time 

= boundary constraint exit time 

- control acceleration 

= state vector, = (jc,,jc 2 ) r 
~ position of experiment platform 
= velocity of experiment platform 
= position of base 
= frequency of base motion, Hz 
= coefficient of optimal control u 
= Lagrange multiplier ° p,J 

= frequency of base motion 
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Table 1 Microgravity space experiment 
acceleration environment (from Ref. 7) 


g/g 0 

F, Hz 

Source 

Quasisteady or dc acceleration disturbances 

to - 7 

0- 10~ 3 

Aerodynamic drag 

I0" 8 

0- 10" 3 

Light pressure 

10- 7 

0 - to - 3 

Gravity gradient 

Periodic acceleration disturbances 

2x I0- 2 

9 

Thruster fire (orbital) 

2 x 10~ 3 

5-20 

Crew motion 

2x 10- 4 

17 

Ku-band antenna 

Nonperiodic acceleration disturbances 

10~ 4 

1 

Thruster fire (attitude) 

io - 4 

l 

Crew pushoff 
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u = F t /m 

Fig. 1 Kinematic representation of the experiment spacecraft system. 


The purpose of this paper is to explore the limitations on 
vibration isolation for space experiments, rather than to de- 
velop the actual control algorithm. Thus, the ideal vibration 
actuator is assumed, and an optimal control is formulated. 
The optimal control problem is solved for a sinusoidal excita- 
tion to obtain the minimum acceleration trajectory. A subop- 
timal solution that gives results close to optimal is also 
explored. 


II. Kinematic Formulation 

For this analysis, a one-dimensional theory is developed. 
Clearly, the actual system required will be multidimensional 
so this work is preliminary in nature. Consider a one-degree- 
of-freedom system, the experiment, as illustrated in Fig. 1 
with position x(f)- It is connected to the spacecraft by umbil- 
icals, such as power or fluid lines, and by a vibration isolation 
actuator. A similar geometry is discussed by Genkin et al., 13 
with stiffness and damping as well as an active vibration 
isolation actuator. However, that system has one side fixed to 
the ground and the forced mass in motion. Although the 
spacecraft actually has a finite mass, it may be considered to 
ha v e infinite impedance for this analysis since the spacecraft- 
to-experiment weight ratio is very large. Thus, the spacecraft 
acts as an external base motion y(t) transmitting forces 
through the umbilicals and the actuator. 

This representation reduces the problem to a kinematic 
one. Onboard the spacecraft, available interior space for the 
experiment is limited. The walls around the experiment, which 
should not be contacted, constrain the maximum translation 
of the experiment, or stroke, to a fixed distance L. The base 
motion y(t) imposed on the walls, spaced to permit a stroke of 
L, forms the problem constraints. The problem of vibration 


isolation/attenuation becomes one of finding the optimal tra- 
jectory (minimum acceleration) given the constraint condi- 
tions (moving walls). 


III. Optimal Control Formulation 

The objective is to formulate and solve the optimal control 
problem for minimum experiment acceleration trajectory 
in time. Let the experiment acceleration be denoted as u 
( u = F v lm). Then the cost function / to be minimized is 

Z = j* « 2 d/ (1) 

with the constraint 

AO - L <> x(t) £ ></), 0^/ (2) 

for a given base motion >{f). 

This problem is examined for harmonic base motion at a 
single frequency. Let ></) have the form 

AO =A[ 1 cos(cu/)J (3) 

with the half-period T = n/ cu. The cost function J now sim- 
plifies to 

y = j^ u 2 dt (4) 

due to the periodicity of the problem. Also, the constraint 

becomes 

A[ \ — cos(7tr/r)] -L £*(/) <: A[\ -cos( 7 rr/r)] 0 £ t <L T 

(5) 


over the half-period. 

This problem may be viewed as finding the optimal path 
through sinusoidally oscillating walls, as illustrated in Fig. 2. 
If the base travel 2 A is smaller than the space L, the minimum 
acceleration is zero and the problem trivial. However, if the 
base travel 2A is larger than the space available for vibration 
isolation L, then the optimization problem has active inequal- 
ity constraints on the state variables. The solution to this 
problem may be attempted using the calculus-of-variations 
approach by adjoining the Hamiltonian with a second-order 
state- variable inequality constraint. This method requires the 
satisfaction of two interior boundary conditions (position and 
velocity continuity) at the junction points of constrained and 
unconstrained path arcs. 14 Because the wall motion is sinu- 
soidal, these tangency constraints require the solution of 
several transcendental equations. Therefore, no closed-form 
solution to the general problem is available. As will be shown, 
under a certain condition the problem can be solved to yield 
an analytic solution. When this condition does not hold, a 
suboptimal solution may be employed. Thus, easy-to-use 
equations and plots for determining vibration isolation limits 
are made available to microgravity experiment designers. 

IV. Analytic Solution 

To obtain an analytic solution to the problem, the con- 
straints are simplified to the boundary conditions 

*(0)=0, x(T) =2A ~L £0 

i(0) = 0, x(T) = 0 (6) 

which an optimal solution clearly must satisfy. 

Define the system state variables 13 as 

x { = x, jc, = x 2 , x 2 = u 


(7) 
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iK7*T 


and the equations of motion become 
i = Ax + Bu 

xr -'- ^ "-[« 0} 


\ 

y(t) — L 


36,00 ^ i 

time ISEC) 

Fig. 2 Optimal path through harmonic walls. 


* r =[0 I] (8) 


,, -™ C . COS ‘ function 's adjoined by the constraint, Eq (7) 
using two Lagrange multipliers. 16 The result is 4 1 ' 

J = j[ [ w2 + ^i(^2-jc,) + x 2 (u — x 2 )] dr (9) 

which is the general functional for this problem. 

Define the Hamiltonian as 

// = w 2 + x 2 + X 2 u ( 10 ) 

cquat, onTar? ^ Ca ' CU,US ° f Variations ’ the 

f dH 

*1 = -t— = 0 

dx, 

f dH 

<3*2 

n <3# 

0 = ^ =2 " +/l * (II) 

Solving these gives 

^1 = c 0 

^2 = C 0 I + Cj 

u = I<V-jCi (12) 

Imposing the boundary conditions of Eq. (6) yields 

f„= -24(2 A -L)/T 2 

C| = — 12(2/1 — L)/T 2 (1 3) 


(OT TOUrse if 2A L< 0. the optimal trajectory is clear! 
onhis pajcctoty b r «>o'-™“"-.qua re (rats) acceleratio, 
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V. Conditions on the Analytic Solution 

The condition under which Eq. (14) satisfies the inequalif 
constraints can be obtained by expanding Eq. (5) as a Taylo 


1 

= ~(-\ _±( n ‘\* 

At) 2a\t) 


A t ) 2 A t ) + ( ,6 

If we combine Eqs. (5), (14). and (16), x{,)Zy«) becomes 

(17) 

For small t/T, 

3(2/1 — L) 5 An 1 /2 
which yields the condition 

T>[2-(n 2 j(>)]A ( , 8) 

The symmetry of the optimal trajectory and inequality con- 

STjrsrrr “it* ,his is a . lso the sufficiem ^ 

time * satisfy the inequality constraint near the final 


"Wop, - — 2(2 ^ — , + 6(2 A -L) 
r 3 t 2 


VI. Suboptimal Solution 



tkfi« S l l K b a. Op,,mal . SOlUti0n ‘° Eq (4) that automatically sa- 
tisfies the inequality constraints of Eq. (5) is y 


(14) 

*. U b(0 =(A -L/2)[l -cos(a//DJ, 2A - L > 0 ( 19 ) 

...L.L.L 1 ... 


2A -L 2.0 


M .ub = (-4 - L/2)(n 2 /T 2 ) cos(jr;/T) 


and rms acceleration 


rms(*, ub ) = ( N /2/4)n 2 ((2/f - L)/T 2 ] (21) 

When Eqs. (15) and (21) are compared, it is clear that the 
suboptimal solution is only slightly inferior to the optimal 

[rms(x opt )/rms(i, ub )] = (4^/6 /ji 2 ) ss 0.9927 (22) 


A comparison of the optimal and suboptimal accelerations 
and trajectories when condition (18) holds is shown in Fig. 3. 
The optimal solution when condition (18) is invalid is a 
combination of a linear control law and wall following trajec- 
tories. 


“opi,(0 = i 


(An 2 /T 2 ) cos(nt/T) y 
a -2a t/T t * 

( An 2 /T 2 ) cos (nt/T), 


0 ^ t <; t* 
t* <t <T-t* 
< T 


(23) 


o 
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Fig. 3 Comparison of optimal and suboptimal acceleration histories. 


where /* and a are determined by the solution of several 
transcendental equations expressing position and velocity 
continuity at /*. It is clear that the rms of this optimal 
solution must be bounded between that of the suboptimal and 
the optimal with the inequality constraints of (5) dropped 

<yn> ^ 

3 4641 ^ <rms ( u op. 1 ) < 3.4895 < —yr-^ <24) 


The suboptimal solution can obviously be employed as an 
excellent approximation to the optimal, Eq. (23), when condi- 
tion (18) is not satisfied. 

VII. Limitations on Isolation 

The primary purpose of this paper is to determine theoreti- 
cal limits to vibration isolation. Although a dimensionless 
plot could have been produced, it was felt that a few typical 
dimensional plots would be of more use to designers of 
microgravity materials science experiments. Figures 4-6 
present the curves of the minimum experiment acceleration vs 
base acceleration at constant frequencies for stroke limits of 
5, 10, and 20 cm, respectively. 

The horizontal axis gives the rms base acceleration calcu- 
lated from 

nns( y) = (J 2I2)co 2 A (25) 

The vertical axis is the minimum experiment acceleration 
from Eq. (15) when Eq. (18) holds and Eq. (21) when it 
does not. 

In Figs. 4-6, the minimum experiment rms acceleration at 
any given frequency is zero (isolation) until the base displace- 
ment amplitude equals one-half the maximum stroke possible. 
The minimum experiment rms then quickly rises with in- 
creases in base acceleration and asymptotically approaches 
the zero-vibration-reduction line, at 45 deg. Along this line, 
the base and experiment act as if they were rigidly coupled 
together and have the same acceleration. 

The primary limitation is the length of stroke allowed 
between the experiment and the base. As an example calcula- 
tion, consider an rms base acceleration of 1 x 10" 3 g 0 at a 
frequency of 0.06 Hz. The base travel is given by 

2A =2 v /2Inns(> ; )/c t > 2 ] (26) 


O 
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Fig. 4 Minimum experiment acceleration given disturbance frequency, 
stroke => 5 cm. 


Fig. 5 Minimum experiment acceleration given disturbance frequency, 
stroke = 10 cm. 





o 



Fig. 6 Minimum experiment acceleration given disturbance frequency, 
stroke = 20 cm. 1 J * 


which has the numerical value of 19.5 cm for this example 
For a stroke L = 5cm, Fig. 4 shows that the experiment 
acceleration will be at least 7 x 10- go- When the stroke is 
increased to L = 10 cm, Fig. 5 gives a value of 4.5 x 10— g a 

The last case is a stroke of L = 20 cm. Figure 6 indicates 
that the minimum acceleration is zero. In this case, L > 2A so 
the stroke is large enough to accommodate the full sinusoidal 
motion without wall contact. 

As an alternative to the plots, Eqs. (15), (18), (21), and (25) 

may be used directly. These can be simplified for this 
purpose to 


rms(x) 

rms(y) 


M) 


0,A < 

2 A 


H) 
KH 


< 1 (27) 


0 , 




with 


2A =[2 v /2rms(y)/co 3 ] 


VIII, Conclusions 

This paper has developed a kinematic formulation for the 
microgravity space experiment problem in one dimension, 
t-urther, two solutions, one optimal and the other suboptimal 
but very close to optimal, have been obtained. These permit 
plots of vibration attenuation for given levels of available 
space. For the sinusoidal oscillation assumed here, the experi- 
ment could be completely isolated if sufficient space were 
ava'lable. Unfortunately, the low-frequency motions (0- 
,.° Hz ) would require motions with a length much larger 
than possible aboard spacecraft. 

Plots of the type developed here are intended to assist 
microgravity experiment designers as well as vibration isola- 
tion engineers. These plots represent the ideal vibration isola- 
tor. Real systems will not be able to attain the ideal for 
several reasons. The actual motion will have several frequency 


components as well as a random component. The random 
component alone will ensure that the full space L cannot be 
employed. Some safety space will have to be allocated to 
prevent occasional wall contact. Any real active control sys- 
tem will have some nonideal characteristics. The sensors 
employed m the active control feedback loop will have some 
errors as well. The authors of this paper are aware that the 
development of a very low-frequency accelerometer is 
difficult. In spite of this, we are optimistic about the levels of 
vibration isolation discussed here. It seems reasonable to 
believe that approximately 75% of ideal isolation is possible 
with an actual control system and actuator. 
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mega— Gauss— Oersted). This design was not only large and heavy, but could not be built from 
a single piece since magnet manufacturers do not make sizes larger than 2 inches. The cost 
and difficulty of assembly ruled out an actuator using multiple magnet segments. It therefore 
became necessary to design the Lorentz actuator using a smaller core gap than is conventionally 
used. Usually this gap is large to reduce magnetic flux leakage across it so as to yield an 
actuator that will produce a force independent of coil position. It was hypothesized that this 



leakage could be substantially reduced by saturating the actuator’s core. This could only be 
verified, short of building a prototype, via finite element analysis. A commercially available 
finite element analysis package, MAGGIE, 
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was used to "test" a large number of designs. The final actuator, shown in Figure 1, has an 
outer diameter of 1.95 inches and a 4 inch length. This actuator has the following features [1] : 

* Long Stroke: The actuator has two inches of stroke. 

* Position Independence: Over the entire stroke, the actuator’s gain is almost 

independent of position. For a constant coil current, this means that the actuator 
force is the same irrespective of the axial position of the coil. This is achieved by the 
design since the maximum flux density across the core gap is only 7% of the maximum 
flux density across the pole face gap. 

Current Linearity: The average flux density in the effective air gap remains constant 
with variations in the coil current between the upper and lower limits. This is achieved 
through the large reluctance of the permanent magnet in the electromagnetic flux circuit 
and the saturation of the core. 

Force: A maximum force of 1.50 lbs is produced by this actuator with a coil current of 
2.5 A. 

* Materials: The permanent magnet is neodymium iron boron, which has a very high 
maximum energy density product of 35 MGOe. The circuit material is a high 
permeability nickel-iron alloy that saturates at 1.50 Tesla. These materials permit a 
compact design. 

The experimental results have confirmed the soundness of the design approach [1] . 
Figure 2a shows the actuator force plotted versus position for a number of values of coil 
current. Note that the actuator’s force is fairly independent of the coil position over the 
actuator’s operating range (0.5 to 2.5 inches). Note also that the actual forces are larger than 
the predicted forces, but still within 20%. Figure 2b shows the same data in terms of actuator 
force plotted versus coil current for different positions. As shown in the figure, the actuator has 
a high degree of linearity with respect to current. Note that the actuator gain (slope of the line 
in Figure 2b) is fairly independent of coil position and is approximately 0.6 lbf/amp. 
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Force Produced by the Lorentz Actuator 

Legend Indicates Coil Position 


Force (Ibf) 



Current in Coil (A) 


Figure 2: Actuator force vs. current. 
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3. MULTIPLE DEGREE OF FREEDOM ACTUATOR DESIGN 

3.1 Introduction 

The University of Virginia also examined the design of multiple-^egree-^f-freedom 
actuators for microgravity vibration isolation. The fundamental constraint on isolation 
performance to be considered during actuator design is the available working envelope [2,3]. 
The implications of this constraint on active isolation were examined by the University in two 

journal publications [2,3]. 

Figure 3 shows the relationship between the envelope (peak-to-peak displacement and 
frequency for several sustainable RMS acceleration levels [4], The graph is for a 



Figure 3: Peak— to— peak displacement vs. frequency 


one-degree-of-freedom case and assumes sinusoidal vibrations, but the relationships are 
acceptable for order of magnitude estimates even if these assumptions are relaxed. No 
definitive specification of the required isolation levels or frequency range exists. The proposed 
US Space Station usable specification is also shown in Figure 3. It is claimed that vibrations 
below this curve will not adversely affect microgravity experiments. The design examined in 
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this section is an active isolation system with a "reasonable" envelope of 4 inches of travel 
and a sustained 1 ug RMS residual acceleration. It can be seen from the figure that this will 
offer isolation down to 0.002 Hz. The amplitude to which vibrations can be attenuated is 
constrained only by controller design and available instrumentation. Operation at lower 
frequencies, however, requires a larger envelope, which becomes prohibitive in terms of 
available spacecraft space. Another specification for the six— degree— of— freedom system 
considered is a rotational range of 40 degrees. 

A redundant coarse— fine scheme with magnetic suspension was chosen. This design is 
particularly attractive for microgravity applications since it allows the use of magnetic 
suspension while overcoming range— of— motion limitations. The design uses a Stewart 
platform for the coarse stage and a novel magnetic bearing for the fine stage [4,5]. The 
approximate regions of activity in the frequency— displacement plane of these two devices are 
also shown in Figure 3. Both stages act to attenuate spacecraft vibrations, effectively 
reducing vibration amplitudes below their active regions on the displacement vs. frequency 
plane. As an example, it can be seen in the figure that a vibration of the spacecraft with 10 
inches of displacement at a frequency of 1 Hz falls outside the active region and could only be 
partially attenuated. It should be noted that such a large vibration is unlikely. If the 
displacement was only 1 inch, however, the coarse stage would absorb all of it except about 
0.005 inches, and the remainder would be reduced down to the micro— g level by the fine stage. 

The combination of the Stewart platform and a magnetic bearing allows continuous 
isolation at frequencies above 0.002 Hz, and a compact, reliable package suitable for the 
application. These choices and some preliminary design concepts are discussed below in detail 
after a survey of other candidate designs. 

3.2 Survey of Published Designs 

Several designs for 6 DOF levitation are discussed in the literature. While these designs 
do not have the envelope of the proposed coarse— fine design, they might be suitable if a 
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coarse stage is not required. They also deserve examination as alternative designs for the fine 
stage. A comparison of the specifications for these designs is given in Table 1. 


Group 

Trans . 

Rot. 

Force 

Envelope 

Mass 

Actuator 

Sensor 

Honeywell 

5 mm 

±1.6° 

43 N 

27x34x50cm 

36kg 

Mag. Brag. 

Eddy & Flux 

N. Vales 

+5 mm 

<±.2° 

.04 N a 

100xl00xl00cm b 

? 

Lorentz 

Capacitive 

NASA 

+4 mm 

1 + 
CO 

o 

o 

445 N 

30x30xl5cm c 

? 

Mag. Brag. 

Eddy 

SatCon 

+10 mm 

±8° c 

4 N 

40x40xl2cm c 

4.9kg 

Lorentz 

Eddy 

IBM 

+5 mm 

±4° 

32 N 

25x25xl5cm C 

? 

Lorentz 

Optical 

Toshiba 

+2 mm 

±1.5° 

20 N c 

25x25x20cm 

8kg 

Mag. Brag. 

Eddy 


El 

Requirement, not limitation 

b Includes experiment package 
c 

Estimated by authors 


Table 1: Comparison of Published Designs 

Four designs specifically for microgravity isolation have been published. Honeywell has 
a well— developed system called FEAMIS [6] with which they have demonstrated impressive 
isolation performance. The system is designed for the Space Shuttle experiment configuration. 
The University College of North Wales also has a well-developed system [7] designed for the 
European Space Agency experiment configuration. NASA [8] has a well-tested laboratory 
system and has done testing in a weightless environment aboard an aircraft in a parabolic 
trajectory. They also have demonstrated impressive isolation performance for a feedforward 
control system. SatCon [9] also has a laboratory magnetic suspension system. 

Two actuator designs were developed for different applications, but they are mentioned 
here because they are similar and could be easily adapted to the isolation application. IBM 
[10] has a laboratory levitated robot "wrist" which enhances robot accuracy and performance. 
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Toshiba [11] has a satellite antenna pointing system which is fully developed. Both devices 
have demonstrated positional accuracies on the order of 1 /xm. 

Isolation of vibrations with large amplitudes — typically occurring at lew frequencies — 
requires a large translational range. SatCon’s system has the largest range, but there is a 
significant tradeoff with the device’s force capability. A coarse-fine approach would allow both 
a large range, provided by the coarse stage, and a high force capability, since the levitation 
gaps are small. There is no available data on the rotational range requirements of the 
application. Isolation with an umbilical disturbance may require a high force capability , as is 
offered by the systems from Honeywell, NASA, IBM, and Toshiba. Volume and weight should 
be minimized in any spacecraft. SatCon’s, IBM’s, and Toshiba’s systems offer advantages in 
envelope volume and weight. 

The choice of the actuator technology between Lorentz force and magnetic bearings for 
MDOF isolation systems is not a clear one. Lorentz actuators offer linearity, simplicity, 
open loop neutral stability, and compactness. Magnetic bearings offer higher force capability 
and lower power consumption, particularly if gaps are minimized. 

Four position sensor technologies offer promising performance. Eddy current position 
probes are simple and robust, but bulky and heavy for large gaps. Capacitive sensors are 
simple and lightweight, but can be noisy in unconstrained environments. Optical lateral effect 
photo— diodes are compact and quiet, but they require substantial supporting electronics. Hall 
effect flux sensors can be used with magnetic bearing designs both to linearize the control 
problem and to measure position. 

3.3 Coarse Stag e 

The Stewart platform is a six degree— of— freedom parallel manipulator which has been 
used extensively in aircraft cockpit simulator applications. Figure 4 shows the mechanism in 
the proposed configuration [4]. Six linear actuators (legs) connect a base (bottom) to a 
platform (top). The base would be mounted in the spacecraft and move with it, while the 
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platform would track an inertial reference frame. Stepper motor driven ball lead-screws are 
proposed as actuators. 



Figure 4: Coarse stage isolation actuator 

This mechanism was chosen over other candidates such as a carriage/gimbal assembly, 
or a serial linkage mechanism, because it has the following features: 

* Inherent rigidity: The parallel connection of the actuators gives the mechanism rigidity 
on the order of the extensional rigidity of the actuators. For the proposed actuators, this 
will allow controller design to ignore the dynamics of the mechanism. The effects of 
"umbilical" connection to the platform will also be negligible. 

* Determinate inverse kinematics: The actuator lengths required to achieve a prescribed 

orientation are found directly from a coordinate transformation from the base to the 
platform frame. This is seldom the case for a serial linkage. This will also simplify 


control. 



14 


* Compactness: The configuration proposed here places the fine stage on top of the 

platform for convenience in testing. A fully developed implementation could locate the 
fine system and microgravity experiment in the space between the base and platform, 
resulting in a compact package. 

The Stewart platform has some disadvantages that must be considered. It is nonlinear in its 
response to actuator lengths, its general direct kinematics have not been discovered in closed 
form, and it has singularities in its operational space. The first two problems can be overcome 
with digital controls. The singularities, which are points or loci where the mechanism gains a 
degree of freedom and the actuators can lose control of the platform, must be addressed by the 
design. 

A simulation code has been written to allow exploration of the design alternatives. 
Results indicate that our specification (4 inches translation, 40 degrees rotation) will be 
achievable with actuators 10.5 inches long in the retracted position, and with 9 inches of stroke. 
The simulations have confirmed that singularities are safely outside the working envelope. 
Commercial actuators with the required range, load capacity, speed and acceleration have been 
identified. 

3.4 Fine Stage 

The magnetic bearing proposed has two parts: a stator which is attached to the Stewart 
platform, and a surrounding "flotor" to which the experiment is attached. The proposed stator 
[5] is illustrated in Figure 5. It has twelve pole pieces and coils arranged around the surface of 
a cube. The cube and pole pieces are ferromagnetic. Each pair of pole pieces and the region of 
the cube to which they are attached comprise a typical "horseshoe" electromagnet causing an 
attractive force toward the nearby flotor. Magnetic flux through the center of the cube causes 
an imbalance in the flux levels of a pair of pole pieces, resulting in a net torque on the flotor. 
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In the proposed design differential Hall effect sensors are located in the base of each pole 
piece to measure the local flux. All electrical connections will be to the stator. 

The flotor concept is illustrated in Figure 6. Three ferromagnetic bands are rigidly 
attached to each other, but form independent flux paths. The bands are thicker in the region 
near the pole pieces to avoid saturation. Flux which passes through the center of the cube is 
returned through the remaining portions of the bands. 

Four mounting posts are attached to corners of the cube, and pass through clearances in 

the flotor. These posts could carry cooling fluid to be circulated through the stator if it is 
required. 

This configuration was chosen over other suspension approaches such as Lorentz 
actuators or magnetic actuators located on the periphery of the experiment package because it 
has the following advantages: 

Compactness: The high force capability of the magnetic bearing relative to a Lorentz 
actuator of similar size and power consumption suits the application. Testing in earth 
gravity will be facilitated, and suspension during launch to protect sensitive 
instrumentation may be feasible. Also, the rigid structure required to mount actuators 
around the periphery is avoided. 

Force/torque balance and rotational range: Actuators capable of the required forces 

mounted on the periphery of the experiment are capable of torques far greater than is 
required, and they limit the rotational range of the experiment. The proposed design 
approach brings the relative force/torque magnitudes closer to the requirement, and 
allows substantial rotational range. 

Integral sensor capability: Compact semiconductor magnetic flux sensors (Hall effect or 
magneto-resistive) can be utilized both to stabilize the system and to infer relative 
position. No elegant integrated approach is known for Lorentz actuators. 
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3.5 Predicted Performance 

The specific design examined at the University has a center cube of 2 in. on a side, pole 
faces of 1 x .5 in., and pole length of 2 in. Maximum current is determined by allowing a coil 
current density of 5000 amp/in 2 which is known to be conservative from previous designs. The 
gap in the centered position was chosen to be .125 in. plus an allowance of .030 in. for inclusion 
of flux sensors and a protective layer on the inside of the bands. The resulting performance of 
the design is presented in Table 2. The 53 N force is a continuous worst case, with the stator 
moved away from the flotor in the direction of the force. The continuous force capability in the 
centered position is 175 N. Intermittent force capability is limited only by the current 
capability of the amplifiers, and the saturation limit of the magnetic material used. Using 
Vanadium Permandur with this design would enable 1000 N force before saturation. Of 
the 4.5 kg mass, the flotor comprises only 1.2 kg. 


Trans . 

Rot. 

Force 

Envelope 

Mass 

+3.2 mm 

+7° 

53 N 

15x15x15 cm 

4.5 kg 


Table 2: Specification of UVA Design 


When compared with the designs presented in Table 1, the UVA design has several 
advantages. The envelope is substantially smaller than any of the previous designs, while the 
performance is similar. In addition to saving space, this compactness allows the flotor to be 
naturally rigid, and thus avoids control problems with structural dynamics. The design is 
lighter than other designs for which data were available [5,6] . 
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4. CONTROL SYSTEM DESIGN ISSUES 

4.1 Introduction 

The control issues of active microgravity vibration isolation were another area of 
investigation at the University. The thrust of this research has been the design of 
feedback/feedforward controllers using modern control synthesis. As part of this investigation 
we also examined passive vibration isolation analogies. In addition, a control architecture for 
the six-degree-of-freedom actuator discussed in the last section was proposed. 

Active isolation systems for microgravity and pointing applications have been designed 
and constructed by many investigators. These systems generally use conventional PD control 
of a noncontacting actuator, either Lorentz or electromagnetic, to achieve low frequency 
disturbance attenuation. While an actual microgravity experiment may require umbilicals, the 
isolation systems designed and tested so far cannot provide isolation for such an 
experiment. These systems achieve their performance by the very low stiffness made possible 
by low gain feedback of the relative position of the experiment to the experiment rack. 
Without an umbilical, this stiffness may be set by the designer at will. However, when an 
umbilical is present, the umbilical stiffness presents a lower bound on achievable stiffness unless 
the feedback loop is used to introduce a negative stiffness. The University has concentrated its 
work on the design of control systems for the generic (i.e. with umbilical) microgravity isolation 
problem. The University has set the following specifications for an active microgravity 

isolation system [12] : 

(1) Unity transmissibility from D C. to 0.001 Hz so as to prevent the experiment from 

impacting its enclosure’s walls. 

(2) At least 40 dB attenuation above 0.1 Hz. 

(3) Both stability and performance robustness with respect to changes in umbilical 
experiment properties, non-collocation or misalignment of sensors and actuators, 
center-of-mass uncertainties, and unmodeled cross coupling between the degrees of 


freedom. 
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Robustness refers to the ability of the control system to perform satisfactorily when the 
true plant varies from the nominal plant. Performance requirements of the type (2) for 
rotational degrees of freedom have not yet been specified by NASA or microgravity users, to 
our knowledge. 

4.2 Passive Isolation: An Analogy 

The design of an active vibration isolation system for microgravity space experiments 
was examined from an analogy to passive isolators [12]. It should be noted that the primary 
reason for pursuing an active rather than a passive system is not the increased flexibility in 
loop shaping accompanying active control, but the limitations of passive isolation systems. The 
stiffness of the umbilical precludes achieving a soft enough support so as to meet the isolation 
requirements for indirect (transmitted through the umbilical) disturbances. Also, a passive 
isolation system cannot isolate the payload from both indirect and direct (onboard the 
experiment) disturbances. An active system allows these limitations to be overcome. For 
example, an active system permits the insertion of a negative stiffness spring in parallel with 
the umbilical. Note, however, that this approach, i.e. lowering the stiffness, requires the near 
cancellation of the umbilical’s stiffness with that introduced via feedback. If the negative 
stiffness exceeds that of the umbilical, the equivalent stiffness of the system will be negative 
and the system will be unstable. It is not surprising then that the introduction of negative 
stiffness via the controller has no robustness whatsoever. A focus on equivalent stiffness in 
isolation system design thus leads to control systems which sacrifice robustness for 
performance. In addition, a design which achieves isolation through lowering the system 
stiffness cannot attenuate direct disturbances over the same frequency band. 

From a vibration engineering viewpoint, an alternative means of achieving rejection of 
disturbances is to fasten the experiment rigidly to an inertial structure. While there is no such 
structure in space, it is possible to achieve this effect by high gain feedback on inertial 
experiment position. This inertial position feedback acts like a very stiff spring tying the 
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experiment to inertial space. While such a controller may meet the 0.1 Hz 40 dB specification, 
it wiU not satisfy the specification on unit transmissibility [12]. If an inertial position 
feedforward loop is added, this problem can be eliminated. However, this method would be 
difficult to use effectively for multiple-degree-of-freedom isolation. 

Another method of fastening the experiment to inertial space examined by the 
University is the use of inertial damping via feedback. By feeding back the inertial experiment 
velocity with a high gain it was shown for an example problem that it is almost possible to 
achieve both the 40 dB and the unity transmissibility specifications without resorting to 
feedforward. Unfortunately, the roll-off rate is approximately 20 dB/decade, so that both 

specifications can not be simultaneously achieved [12] . 

Another passive analogy examined was the lowering of the natural frequency of the 
umbilical by increasing the experiment mass. An increased experiment mass would attenuate 
direct disturbances as well as those transmitted through the umbilical. In addition, at 
frequencies below the natural frequency of the umbilical-mass system, the isolation system 
would have unity transmissibility. Of course, for space applications any additional mass is 
very costly. To lower the natural frequency by an order of magnitude would require increasing 
the experiment mass by a factor of one hundred. Clearly, it is not practical to accomplish 
increased isolation through the addition of real mass. However, it is possible to increase the 

effective mass of the system through feedback [12] . 

To summarize, the passive isolation analogies examined yield some insight but they fall 
short as design approaches on three counts: (1) they do not have flexibility to shape the 
response so as to achieve the performance requirements, (2) they cannot be easily generalized to 
multi— degree— of— freedom problems, and (3) they completely ignore the robustness problems 
inherent with active control systems. 

4.3 Classical Control Design 

A one-dimensional isolation problem, shown in Figure 7, was examined using a classical 
controls loop-shaping approach, to gain insight into controller design and limitations. System 
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Figure 7: One-dimensional isolation problem 

parameters were chosen to be representative values which yielded a low natural frequency 
(k/m = 0.1, w n a 0.05 Hz [0.316 rad/sec]), and damping was assumed light (C = 01)- In the 
following discussion the variables d, x, and u represent experiment rack position, payload 
position, and control force, respectively; and it is assumed that the only available measurement 
is payload acceleration. The problem is to design a feedback controller, satisfying the following 

specifications: 

1. Above 0.1 Hz the payload acceleration x(t) should be 40 dB below the spacecraft 
acceleration d(t). 

2. Below 0.001 Hz the payload vibration x(t) should track the spacecraft vibration d(t) to 
within 10 percent, in order to prevent collision of the payload with the walls of the 
experiment rack surrounding it. 

3. The payload should track perfectly the DC motion of the spacecraft, where no relative 
motion can be tolerated. 
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4. The loop gain of the system (plant and controller) should be less than 0.1 above 200 Hz, 
to avoid controller excitation of spacecraft- or payload flexible modes. 

5. The payload acceleration should be less than or equal to 1.1 times the spacecraft 
acceleration at all frequencies. 

6. Large phase margins should be attempted at all crossover frequencies. 

The system equation of motion is 

rax + cx + kx = cd + kd — u 

and the system transfer functions are 


s a(s) = 


cs+k 

ms 2 +cs+k I 


s 2 D(s) + 


-s 


ms +cs+k 


U(s) 


with a system block diagram as shown in Figure 8. R(s) represents the input disturbance (rack 

acceleration, Laplace domain), C(s) represents the payload acceleration, H^s) represents the 

2 

— s 

controller, and U(s) represents the control force. G(s) and c s + k G ( s ) are the tw0 P lant 
transfer functions, and H(s) is defined as indicated in Figure 8 for convenience. 


G(s) 



H(s] 


Figure 8: One-dimensional isolation system block diagram 
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The first five specifications can now be re-expressed, respectively, in the following form: 

1. C ( s ) 0.01 above 0.1 Hz (0.628 rad/sec) . 

m 

2. 0.9 < C ( S ) < 1.1 below 0.001 Hz (6.28- 10 -3 rad/sec) . 

R^J 

3. lim C ( s ) =1. 
s-+0 R[sJ 

4. | H (s) G(s) | < 0.1 above 200 Hz (1256 rad/sec) . 

5. C ( s ) <1. 1 for all frequencies. 

R[sJ 

In order to use the classical approach efficiently, the above specifications must be 

reduced to loop-gain form. This reduction yields, respectively, the following: 

1. G < 0.01 above 0.628 rad/sec, which in turn requires roughly that |H| >22 

1+HG 

(i.e., greater than 100/G1) at that point. 

2 The second specification, 0.9 < G < 1.1 below 6.28*10 rad/sec, is roughly 

1+HG 

equivalent (since |G| ~ 1 below to the requirement that |HG| < 0.1 below 
6.28 *10 -3 rad/sec. 

3. lim | HG | = 0. 
s-* 0 

4. | HG | < 0.1 above 1256 rad/sec (same form as before). 

5. | HG | > 7 in the vicinity of w n (where | G | « 6.5) to reduce the transmissibility to about 
unity in that region. 

Standard loopshaping methods can now be used in a straightforward manner. See 
Figure 9 for asymptotic Bode— a plots of the specifications and of G(s); and for "first— pass 
plots of the loop gain "L(s)" [i.e., H(s)G(s)] and of the controller Hj(s). 




Figure 9: Asymptotic Bode-a plot of plant G(s) and "Grst pass loop gain L(s) 

The controller developed first led to a transmissibility resonance at u> n (not shown) so a filter of 
2 

f orm (S+0-316L was added, resulting in the following controller: 

7 ^+ 0.1 

1.2-10 3 (s + 0-25) ( s + 12) (s + 0 . 316 ) 2 (s 2 + 0-063 s + M l 
H 1 (S) ( s 7 0.009) (s + 1 T (s + 4) (s 2 + 0.1) 

Figures 10a, b.c represent loop gain, controller, and transmissibility plots, respectively. The 
control meets all specifications except for the goal of no more than a transmissibility of 1.1 at 
all frequencies; and this speciDcation is almost met. The two phase margins associated with the 
above controller are 59° and 88°, respectively. 
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TRANSMISSIBILITY 




Figure 10: (c) Closed-loop transmissibility plot designed using classical control methods 
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From the above analysis the following conclusions can be drawn. 

1 . The requirement that jg| be less than some fraction 0 2 above some frequency at, [spec 
#1] means that the open loop gain L(s) [i.e., H(«)G(.)] must be greater in magnitude 
than l G(s) above u> 2 . This means that there is a tradeoff between 0 2 and PMj. The 

smaller 0 2 is, the smaller PM, can be. That is, the better the disturbance rejection 
above w 2 , the lower the achievable phase margin PM,. 

P i $ PM. |. (Lowering 0 2 will also tend to reduce PM 2 , but not as directly.) 

- 2 1 

2. Raising cj 2 will improve PM, but degrade PM 2 - 

w 2 t^pm 1 t, pm 2 i 

3. The requirement to keep |L(s)| below 0 3 above n> 3 [spec #4] (so as to avoid exciting 
higher modes) has a cost in terms of PMj : /JjWPMji. 

4. Raising o/ 3 raises PM 2 : t^ PM 2 T- 

5. The requirement to hold |£| within some fraction /?, of unit transmissibility below 
some frequency u, [spec # 2] means that |L(s)| must be less than ^ |G(s)| below 

There is, then, a tradeoff between 0, and PM,: /?,HPM,i. (Changing 0, does not 
significantly affect PM 2 -) 

_ 6. Lowering u, will improve PM,: w,t^PM,T. 

7. Lowering the natural frequency « n eases the difficulty in obtaining adequate PM, by 

lowering the constraint at u 2 (see Figure 9) at u> 2 : 1* PMJ- This means that 

reducing the physical umbilical stiffness or increasing the physical payload mass will 

make for an easier control problem. 

8. The problem can be simplified, and both PM, and PM 2 can be increased, if the 
umbilical is damped such that the resonance near is small. (Refer to spec #5, p. 

~ 23.) 

9. The controller need not have zero gain at DC to be acceptable, as long as 

lim s 2 H (s) = 0. The controller may have a low frequency asymptote with slope 
s-dO 1 

— -1, 0, or greater. 
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Although the classical approach is not readily extendable to the MIMO problem, it does 
provide some useful insights for informing the extended H 2 synthesis approach that we will 
examine in Section 5. Weighting x(t) more heavily above o> 2 is analogous to lowering 0 2 (see 
conclusion #1 above), so that better disturbance rejection is achieved at the expense of phase 
margin (esp. PM 1 ). Weighting the control u(t) more heavily at higher frequencies corresponds 
to trying to reduce /? 3 , so that a reduction in controller bandwidth is purchased at the expense 
of phase margin (PM 2 ) (see conclusion #3). At the lower end of the frequency spectrum, 
increased weighting of relative displacement (x-d), reduced weighting of absolute acceleration 
(x), or increased weighting of the control (u) each corresponds to attempting to lower /?p at the 
cost of reducing PMj (see conclusion #5). Since an acceptable controller can have large, even 
infinite, DC gain (see conclusion #9) it is not necessary to weight u(t) highly at low 
frequencies. In fact for phase margin considerations (PM^ it may be best to have "cheap 
control at low frequencies, as previously noted (see conclusion #5). Unity transmissibility, 
then, could be "requested" at low frequencies by a relatively high low-frequency weighting of 
relative displacement. 

4.4 Extending to the Multiple Degree— of— Freedom Probl em 

The University has extensively examined the design of multiple— input-multiple-output 
(MIMO) controllers for the multiple-degree-of-freedom active isolation problem. This work 
will be examined in detail in the next section. Here, we will introduce some of the problems of 
extending single-input-single-output (SISO) methods to MIMO problems by examining a 
simple multiple — degree — of— freedom benchmark problem [12], shown in Figure 1 1. 

This problem illustrates how controller design via decoupling an isolation problem into its open 
loop modes, designing controllers for each mode, and recoupling back into the actuators, will 
often result in poor robustness due to unmodeled cross-couplings. This method of design, 
converting a MIMO control problem to a series of SISO problems, is often practiced. The 
example system is composed of an isolated platform (width 0.5 m and height 0.2 m, depth 

OF Pc Of: QUALITY 
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• Assumed Center of Mass 
® Center of Mass 
□ Accelerometer 


Figure 11: Simple multiple-degree-of-freedom benchmark isolation problem 

unspecified), two accelerometers, two actuators, an umbilical, and a translating base. The 
platform may translate vertically or rotate about its center-of-mass. The actuators and 
accelerometers are positioned a distance of q = 0.2 m symmetrically about the assumed 
center— of— mass location. An umbilical of stiffness k (no damping) runs between this location 
and the base. The platform has mass m and inertia I. The equations of motion for the 

platform’s translation x(t) and rotation 0(t) are 

mx + kA0 + kx = fi + f 2 + di 

I'O + kA 2 0 + kAx = (q 4- A)f 2 — (q - A)f t 4- d 2 

where d, and d 2 are the disturbances, and A is the error in the assumed center of mass. The 
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accelerometer readings are 

yi = x — (q— A)# 
y 2 = x - (q - A)0 

The nominal system (A = 0) can be decoupled in terms of the degrees of freedom by the change 
in variables 


F = fi + f2 


M = q(f2 - fi) 


Zi = (y 1 + ya)/2 
Z2 = q(y2 — yi )/ 2 

w hich are nominally the translational force, the moment, the translational acceleration, and the 
angular acceleration for the platform, respectively. The nominal transfer functions for the 

system are then 


Z,(s) = 


m 


s 2 + k 


(F(s) + D,(s)) 


Z 2 (s) 


1 

I 


(M(s) + D 2 (s)) 


For translational motion, the natural frequency of the platform is |k/m . The rotational 
motion of the platform is free since the umbilical is attached to the center-of-mass. To 
compensate the nominal system, feedback can be designed for each mode of the system 
separately, since the system is decoupled. Translational acceleration and velocity feedback are 
first used to add effective mass and damping. 


F(s) = - 



Z,(s). 
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This lowers the 
function 


natural frequency of translational motion, yielding the closed loop transfer 

,D,(s). 


Zl( s ) = 


( m + a ) s^ + cl + kj 


Next, angular deflection feedback is used to constrain low frequency rotational motion and 
some damping is provided. 


M(s) = - 


n b 

7 + T 


Z 2 (s) 


yielding 


Z 2 (s) = 


I s + n s + bj 


D 2 (s) 


where the control system values are in effective units. A control system was designed to lower 
the natural frequency of translational motion from 0.056 to 0.006 Ha with 40% of cr.t.cal 
damping. The controlled rotational motion has a natural frequency of 0.006 Ha with 26% of 
critical damping. This controller design would yield very effective isolation on the nominal 

system. 

The actual close loop poles, however, will be different from the nominal due to the error 
in the center-of-mass A. The poles of the actual system are given by the roots of the 

characteristic equation 

[(m + a)s 2 + cs + k][Is2 + ns + b] - [mA][A{as 2 + cs + k)] = 0 

For the nominal plant (A = 0) , this results in the prescribed natural frequencies and critical 
dampings. However, as the center-of-mass error increases, the poles migrate and the system 
becomes unstable. For an error as small as 6 mm for this system, instability occurs [12], A 
plot of the pole movement vs. error in center-of-mass is shown in Figure 12. This sens.t.vtty 
results from the ill-conditioned character of the designed controller. A proper MIMO 
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controller design might remedy this problem. In any case, an analysis of the problem from a 
MIMO control perspective would indicate the potential instability and the nature of the 

trade-off between performance and robustness. 

In the next section, the BIBO design methods developed at the University of 

Virginia are examined in detail. Special attention is given to the issue of 
robustness . 
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Figure 12: Loci of closed— loop poles as a function of center— of mass error 
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5. MODERN CONTROL THEORY DESIGN 

5.1 Modern Control Methods: An Overview 

Researchers at the University have investigated the use of Linear Quadratic Regulator 
(LQR) and the Kdman-Bucy Filter (KBF) synthesis methods for the design of controllers for 
microgravity vibration isolation platforms [13,14,15]. The LQR method produces a state 
feedback controller which is optimal with respect to the quadratic (two norm) performance 

index 

rOO 

J = x T (jw)Qx(jw) + u T (jw)Ru(jw) d u 

— J 00 

where Q and R are respectively the symmetric (usually diagonal) state and control weighting 
matrices, and x(jw) and u(j u) are the Fourier transforms of the state and control vectors. The 
state (positions and velocities for vibration isolation) satisfies the differential equation 

x = Ax + Bu 

The quadratic performance index of LQR is well suited to this problem since vibration 
isolation quality is usually measured in terms of root-mean-square. However, it has been 
shown by researchers at the University that some modification of the performance function is 
necessary to apply this synthesis procedure to microgravity isolation controller design. State 
feedback for the isolation problem is feedback of experiment positions, velocities, angles, and 
angular velocities. Thus, LQR can only result in (inertial or relative) stiffness and damping 
feedback. As was discussed previously, these isolation techniques cannot yield acceptable 
isolation performance. Thus, an LQR performance index will not yield a satisfactory controller 
unless frequency weighted Q and R matrices are used, or the plant model is changed so as 
to have an acceleration pseudo-state [12]. Either of these methods results in the 
addition of pseudo-states to the state variable model. Frequency weighted Q and R matrices 
are also necessary to achieve robustness. Through choice of the weighting functions, the 
designer can, in essence, shape the control loops. 
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The differential equation above does not include a disturbance term. Consequently, the 
resulting controller is optimal with respect to white noise (a weakness of the LQR machinery). 
Since the power spectrum of the microgravity environment is not of this shape, the LQR 
controller will not be optimal with respect to rejection of the disturbance. Through the 
incorporation of a disturbance model (essentially a shaping filter), the LQR problem may be 
modified to yield an optimal disturbance accommodating (i.e. rejection) controller. This also 
incorporates the addition of pseudo-states to the state variable model. Disturbance 
accommodation may also aid in increasing the controller’s robustness through loop 
shaping. Through the incorporation of the pseudo— states for frequency weighting and 
disturbance accommodation, controllers have been designed by University researchers using 
the standard Algebraic Ricatti equations of LQR— KBF. These calculations have been done 
using batch files written in the MATLAB language [15]. These controllers are then tested for 
robustness with respect to structured and unstructured uncertainties using singular value and 
structured singular value analysis. These analysis tools are the MIMO equivalent of the 
familiar gain margin, phase margin, and root locus robustness tests. Results for a 
one— degree— of— freedom problem are discussed below. MIMO vibration isolation research is 
ongoing at the University. These modern control methods require a considerable degree of skill 
and insight to employ properly. 

5.2 Modern Control Results 

The one-dimensional problem was first expressed in state-space form, with payload 
relative position, relative velocity, and acceleration selected as states. Although many other 
state choices could have been made, these three were chosen to minimize the number of states 
necessary and to maximize the physical intuition possible. The selection would result in a state 
feedback control that respectively modifies the effective umbilical stiffness and 
damping, and the effective payload mass — all being familiar, accessible, and intuitive system 
parameters. Relative, rather than inertial, position feedback would help to avoid exceeding 
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rattlespace limits; and relative velocity feedback would provide a means of damping out system 
resonances. The selection of acceleration as a state was considered desirable due to insight 
gained from the passive control studies. A controller which increases effective payload mass 
(by negative acceleration feedback) would potentially be able to accomplish disturbance 
rejection without unnecessarily sacrificing stability— or performance robustness. 

A second important feature of the problem formulation was the decision to incorporate 
disturbances of two different kinds, the direct (i.e., onboard the experiment) and the indirect 
(i.e., acting via the umbilical). It had been observed that reducing the effective umbilical 
stiffness could aid in indirect disturbance rejection only, but that increasing payload effective 
mass could help reject disturbances of both kinds. Although the primary type of disturbance 
was considered likely to be the indirect, a means was needed to force the LQR-KBF (also 
known as LQG) "machinery" to increase effective mass so as to result in a robust controller. 

Including a direct disturbance provided this mechanism. 

After completing the problem formulation, the next step was to develop a computer 
code for use in design and analysis. A PC-based design code was written in MATLAB to allow 
for accommodation of both direct and indirect disturbances. A large selection of frequency 
weightings and disturbance accommodation filters was made available to the designer. The 
code computes both feedback and observer gains, and also determines the constant feedforward 
(preview) gains for which the theory was developed in [16]. Although the feedforward option 
remains available, subsequent analysis determined that for the present application the 
feedforward gains do not make a significant enough contribution to warrant the additional 
controller complexity required. A number of analysis routines were also written to allow the 
designer to evaluate the resultant designs for purposes of comparison. The number of system 
states, system performance, stability robustness, parameter sensitivity, and observer quality are 
items whose comparisons are facilitated by these routines. 

With the design and analysis tools in place, the next step was to develop the desired 
controller. In order to make the controller as simple as possible, it was decided to begin with 
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the basic LQG approach and to add complexity as needed. At each stage of additional 
complexity an iterative cycle of design and analysis was employed in an attempt to get the 
"best” achievable controller at that level of complexity. 

The basic LQG approach (no frequency weighting, no disturbance accommodation, no 
direct disturbance) yielded a satisfactory controller in terms of performance; but it had almost 
no stability robustness to changes in umbilical stiffness from the nominal (as measured by 
feedback uncertainty). This lack of robustness was due to the fact that LQG found adding 
negative stiffness to be a "cheaper" means of indirect disturbance rejection than adding 
effective mass. No frequency weighting was found which could rectify this problem. 

A direct white disturbance was added in an attempt to force the LQG design 
"machinery" to add effective mass. Although there were some gains in stability robustness this 
was due entirely to changes in observer gain matrix L. The feedback gain matrix K remained 
unaffected (note that this is fundamental in LQG theory and is not a numerical problem), and 
the feedback stability robustness was still unsatisfactory. 

Disturbance accommodation, with a lowpass filter applied to a large direct (white) 
disturbance, resulted in a controller with excellent feedback— and multiplicative input stability 
robustnesses, as measured by singular value checks. The multiplicative output stability 
robustness was unacceptably low if cross— coupling was considered possible between states, but 
structured singular value checks indicated that without cross— coupling the allowable 
multiplicative output uncertainty was quite satisfactory. Since effective stiffness, effective 
damping, and effective mass of the controlled system are uncoupled for the true 
one— dimensional problem, the stability robustness measures of the system were considered 
acceptable. Further, the performance was excellent, easily exceeding the specifications. 
However, the controller gains were still large at higher frequencies where unmodeled system 
modes were of concern (see specification #4). It was therefore necessary to use state— and 
control frequency weighting in an attempt to force the controller to turn off by approximately 
100 Hz (i.e., to reduce loop gain below a magnitude of one) so as to avoid exciting unmodcled 
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flexible modes. To reduce the loop gain at the higher frequencies it was necessary in that range 
(1) to place a high weight on control, (2) to apply low weights to all three states, and (3) to 

reduce the direct disturbance. 

At low frequencies the control weighting was left constant (i.e., "flat"), in an attempt to 
minimize the number of added pseudostates. However, the resulting closed loop system now 
had very poor low frequency stability robustness to parametric uncertainties, even though it 
both retained its excellent performance and now provided the desired low controller bandwidth. 

A classical design approach to the problem provided a simple solution to the robustness 
issue. It was noted that for a controller with acceptable nominal performance the low 
frequency asymptote for controller gain could have slope -1 or 0 or greater (Bode-a, log-log 
scale). Therefore, the control weighting at DC could be zero (filter slope > 0) and the extended 
H 2 synthesis "machinery" could be freed to consider finite or infinite DC controller-gain 
options. This results, however, in the addition of a pseudostate. This change yielded a 
controller that satisfied the design specifications and exhibited good stability robustness to 
parametric and to multiplicative input- and output uncertainties. Considering (for the 
moment )only single-parameter uncertainties, stability was guaranteed for umbilical stiffness to 
within ±99.7% of nominal, and umbilical damping could be essentially unknown. Payload 
mass needed to be known only to within ±65.2% of nominal. Having these initial favorable 
indicators of system robustness the next step was to reduce the controller size. Further 
robustness analysis would then be conducted on the reduced-order controller. 

The controller described above was a ninth-order controller (i.e., had nine states), with 
payload acceleration as its only required input. Other states and pseudostates were 
reconstructed in the observer. To reduce the controller to a smaller order, a routine was 
written in MATLAB in order to permit removing high frequency modes (modal truncation) and 
weakly controllable and -observable system dynamics [17]. The result of applying this to the 
ninth-order controller was a third^rder controller that has all the essential features of the 
ninth-order one. The loop gain, controller, and transmissibility plots for this reduced 
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„ , M ntf , f rom the transmissibility plot that the 
controller are shown in Figures 13a, b,c. Note lrom 

COni • • u i m-4 at 0 1 Hz Notice further that 

transmissibility is unity up to 10-= Hr and that rt rs below ^ ^ ^ ^ 

the open loop and closed loop Bode plots merge at about 100 Ha. - 

the controller has essentially "turned off" by that frequency (see F.gure 13 ). 

There are four basic checks that must be made of any controlled system, nomr al 

.tability, nominal stability, robust stability, and robust performance. These fonr c ec 
rnnsidered below, consecutively. 

The extended Hr synthesis method used for this portion provides an inherent guaran 
of stabrlity for a nomina, plant with full state feedbac. Further, the 

g — that for a perfectly known plant a ^ the 

the system. Thus, nominal stab.l.ty rs assure w. 

Observer itseif is stable. Reducing the controller order removes thrs guara t , 

alue checks verify that both the reduced third-order controller as desrgned and 

agenval checks y the i 00 p gain 

associated controlled system are stable for the normna. plant. A srmple 

i * j. t, Hosed loop system is stable, since 
Bode plot (Figure 13a) confirms the conclusron that the closed P 

known that the loop gain is minimum phase. abated by the closed-loop 

The second necessary check is of nominal performance. As rndrca 

transmissibility plot (Frgure 13c) the nominal performance is quite — o ***** 
,.,ess than 10- spec at 0.1 Ha is surpassed by more than an order of magn tude. Tins 
idlign was intentional, and necessary, srnce plant modeiing errors (open loop system, 

sensors and actuators) will certainly degrade performance margms. 

’Robust stability measures are necessary to determine whether the —op sy^ 
.1,1 remain stable given the anticipated sensor, actuator, and plant parameter — - 
Three different types of robust stability measures were used, for guarantee, ng system s y 
• • t multiplicative output, and feedback uncertainties below certain levels. 

for multiplicative input, multiplicative P pua ranteed 

The multiplicative input uncertainty allowable was found to be equrvale g ^ 

phase margin (interval) of H*‘ , +«M. t0 1 g uaranteed gal " marS1 " ” 
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5 4341 The actuU matins are even iarger (phase margins: [-55-, +»']. **» margins: 
[0,+oo]). Since only one plant output is sensed (viz., payload acceleration), the multiplicative 
input and output robust stability guarantees are identical. A feedback uncertamty measure 
was used to determine guaranteed minimum stability bounds on uncertainties » ***** 
.tiffness and damping, and on payload mass. It was found, as noted previously (p. 34) that 
closed-loop system stability was guaranteed for single-parameter uncertaint.es much larger 
than anticipated. By considering the feedback uncertainty structure, rt was shown t at or 
simultaneous mass, damping, and stiffness uncertainties of ±20%, ±100 ., an 
respectively, system stability could be assured. Higher frequency modes of the system were 
considered not to be a significant concern since the controller bandwidth was bruited urmg 

8 Finally, measures were needed of performance robustness. Structured singular value 
plots were made to find conservative bounds on multiplicative input (and output) unce 
that would not lead to plants with unacceptable performance. Below 10- Hz it was found t at 
for combined sensor and actuator uncertainties of up to ±lf in phase or of ±19% in gam t e 
performance can be guaranteed to remain acceptable. At higher frequencies the guarantees are 

. . r in iVhcisG or of +200% 

much better, so that by 220 Hz uncertainties of up to ±180 ,n phase 

permissible. 

Structured singular value plots were dso used in an attempt to find performance 
robustness guarantees in the face of known parametric uncertainties, but the effort was on y 
partly successful. The checks led to the conclusion that for single-parameter uncertainties in 
stiffness of ±40% both stability and acceptable performance could be assured. However, 
single-parameter uncertainty bounds found by this method on damping and mass were too 
conservative to be useful. Consequently, real parametric studies were conducted on 
plant-uncertainty effects on closed-loop performance. It was determined that closed oop 
performance appeared acceptable for the various combinations of parametric uncertainties 
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examined, with mass and stiffness varied in the intervals [—50%, +100%] and [-20%, +100%], 
respectively, and with damping varied by more than ten times its nominal value. 

The above extended H 2 synthesis — n analysis approach produced a controller that easily 
satisfies the competing demands of the posed 1— D microgravity vibration isolation problem. 
Further, unlike the classical approach, it is readily extendable for use on a 3— D problem. 
Frequency weighting and disturbance— accommodation were both found to be necessary if H 2 
synthesis is to be used in involving the posed isolation problem. Their inclusion, along with a 
judicious choice of states, provides the designer with a powerful and intuitive set of weapons for 
his design arsenal. Disturbance accommodation of a direct disturbance model is necessary to 
force the H 2 synthesis machinery to avoid negative-stiffness solutions. The result was an 
actively controlled system that uses a "smart" form of acceleration feedback to overcome the 
robustness problems that commonly plague the basic LQG synthesis approach. 
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6. EXPERIMENTAL RIG 

6.1 Introduction 

The University began construction of a one-<legree-of-freedom experimental rig to 
demonstrate active microgravity isolation in the fall of 1990. The rig, now completed, was 
designed so as to illustrate active isolation of a tethered mass down to very low frequencies 
(0.01 Hz). This required both a large-stroke actuator and acceleration feedback as discussed in 

Sections 2 and 4. To our knowledge, this is the first microgravity rig to address either tethered 
or large— stroke active isolation. 

6.2 Rig Description 

The experimental rig built at the University of Virginia is shown in Figure 14. The rig 
consists of a 75 lb. steel cylinder representing a microgravity experiment, two air dashpots 
representing umbilicals, an electrodynamic shaker representing the vibrating experiment rack, 
and the large-stroke Lorentz actuator. The steel cylinder is suspended with magnetic supports 
so that it may freely move horizontally along its axis [16]. Similar to radial magnetic bearings, 
each support consists of four horseshoe electromagnets. Eddy current probes sense the radial 
position of the cylinder and complete the magnetic suspension feedback loops supplying 
current to the electromagnets. The supports hold the cylinder firmly in place but produce no 

friction. When the electromagnetic support system is turned off, the cylinder rests on a pair of 
touchdown pedestals. 

The electrodynamic shaker (representing the experiment rack aboard the orhiter) has 
a long peak-to-peak stroke of 6.25 inches. This is the vibration source from which the steel 
cylinder (experiment) must be isolated. The shaker is mounted, via aluminum plates, on a 
concrete block resting on the laboratory floor. The shaker can generate sinusoidal, random or 

impulse waveforms at frequencies down to DC, thus simulating the disturbances typically 
produced on a manned orbiter. 
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Fig. 14. Microgravity isolation rig at the University of Virginia 

The umbilicals connecting a microgravity experiment to the orbiter are expected to be 
flexible hoses and wires. These are modeled by air dashpots with adjustable stiffness and 
damping coefficients. The vibration isolation test rig at the University has been designed so 
that different kinds of umbilicals may be employed, including actual hoses like those used for 
fluid transfer. The large-stroke Lorentz actuator connects the levitated steel cylinder to a 
plate connected to the concrete base. 

The axial acceleration of the cylinder is sensed off a sensory plate using a very low 
frequency accelerometer with a resolution of approximately 1 Mg- The accelerometer signal is 
fed through a low pass filter and a transconductance bipolar linear amplifier to produce the 
required current. This current is applied to the Lorentz actuator to isolate the cylinder from 
the disturbances generated by the shaker. 


45 


The background vibration levels on the concrete base on which the cylinder is mounted 
have been measured over several twenty-four-hour periods, in both the horizontal and the 
vertical directions. These vibrations are of the order of miUi-g's, the quietest period occurring 
from late in the night to earl, in the morning [18]. Operating at this time will yield the highest 

degree of reproducibility in our results. 


6.3 Experime ntal Results 

Preliminary results have been obtained for vibration isolation in the (1-3) Hz 
range, in air dashpot (umbilical) was the only direct connection between the shaker 
armature (space platform) and the cylinder (science experiment requiring isolation) . 

An HP Structural Dynamics Analyzer was used for data acquisition. Figure 15 is a 
typical example illustrating the isolation obtained using simple lowpass acceleration 
feedback. The shaker generated a sinusoidal armature motion at a frequency of 
For this case, the shaker's acceleration had an amplitude of 14,000 *. The cylinder 
had a peak acceleration amplitude of approximately 7,000 pg with the controller "off" 
and 465 pg with the controller "on". Therefore, a fifteen-fold reduction of vibration 

has been obtained through acceleration feedback . 

The control system is now being modified to improve the isolation capability of 


the controller. 
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7. CONCLUSION 

The University has made substantial progress in many areas of active microgravity 
isolation in the last three years. We have primarily addressed the design of actuators and 
control systems for the active isolation of tethered experiments. In actuator research, our work 
has examined electrodynamic and electromagnetic actuators for single and 
multiple-degree— of— freedom isolation, and the use of coarse-fine systems for the practical 
extension of electromagnetic isolation to large strokes. For control system design, we have 
addressed performance limitations, robustness issues, and the use of H 2 methods for synthesis. 
Finally, we have constructed a single-degree-of-freedom test rig and demonstrated active 
isolation of a tethered mass through acceleration feedback. Our research is ongoing and 
several important results are still to be achieved. The University looks forward to continuing 
its work in microgravity vibration isolation and to continued collaboration with NASA Lewis 
Research Center. 

To make a microgravity environment available for space experiments in the near 
future, we recommend the following: 

* The umbilicals to be used to service the experiments need to be identified and 
their properties need to be examined. As the research conducted at the 
University over the last three years demonstrates, the difficulty of achieving a 
microgravity environment is very directly related to the umbilical s 
properties. For multiple— degree— of— freedom isolation, the uncertain coupling 
of degrees— of— freedom through the umbilical may present a challenge to 
controller design. For this reason, it is also recommended that controlled 
umbilical:; be examined. 

* The issue of direct disturbances needs to be addressed. Acceleration feedback, 
like that developed in our work, will be effective against direct disturbances 
as long as the frequencies of these disturbances are below that of the first 
flexible mode of the experiment structure. Perhaps a specification for 
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experiment designers on the frequencies and amplitudes (or power spectrum) of 
allowable direct disturbances can be written. Such a specification would 
require that direct disturbances be acceptable without active vibration control 
for frequencies near and above the first flexible mode. This may help focus 
attention on the issue of direct disturbances and experiment design so that any 
required technology development may begin soon. For example, such a 
specification may result in the inclusion of passive vibration isolation mounts 
onboard the experiment package to isolate the sensitive process from high 
frequency direct disturbances produced by auxiliary equipment (e.g., pumps, 
fans, shutters, valves). 

* The isolation frequency and amplitude requirements of microgravity experiments 
and the microgravity vibration environment of the space shuttle and space 
station need to be better characterized. This is very important in the low 
frequency (0-1 hz) range. Only when these quantities are specified can the 
required stroke of the actuator be determined. If strokes larger than 1 cm are 
necessary, a coarse-fine actuation system should be used. In this case, a 
technology development program needs to be started. The authors believe that a 
significant degree of development may be required for such a coarse-fine 

actuation scheme. 

* A six-degree-of— freedom microgravity isolation system needs to be flown aboard 
the space shuttle in the near future. Only when we start developing actual 
hardware and software for an orbiting isolation system will we make significant 
progress toward practical isolation for space experiments. Vhile we have 
learned a great deal from the experiments conducted so far, many of the 
difficulties that remain cannot be fully simulated or anticipated using ground 


based hardware. 
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ABSTRACT 

This report summarizes the research performed at the University of Virginia under 
Grant No. NAG-3-909 from NASA Lewis Research Center. This research on microgravity 
vibration isolation was focused in three areas: (1) the development of new actuators for use in 
microgravity isolation, (2) the design of controllers for multiple— degree— of-freedom active 
isolation, and (3) the construction of a single— degree— of— freedom test rig with umbilicals. 
Described herein are the design and testing of a large stroke linear actuator; the conceptual 
design and analysis of a redundant coarse— fine six— degree— of— freedom actuator; an 
investigation of the control issues of active microgravity isolation; a methodology for the design 
of multiple-degree— of— freedom isolation control systems using modern control theory; and the 
design and testing of a single— degree— of— freedom test rig with umbilicals. 
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1. INTRODUCTION 

The University of Virginia began research in microgravity vibration isolation in 1988 
under a three-year grant from NASA Lewis Research Center. The goals of this project were 
(l) to develop new actuators for use in microgravity isolation, (2) to investigate the design of 
controllers for multiple — degree — of— freedom (MDOF) active isolation, and (3) to construct a 
single— degree— of— freedom (SDOF) test rig with umbilicals. The Principal Investigator for the 
first two years of this project was Dr. Paul E. Allaire. Dr. Carl R. Kno.^pe became the 
Principal Investigator for the final year of the grant. Other faculty working with the project 
included Dr. Robert H. Humphris and Dr. David W. Lewis. Several graduate students have 
aided in this effort: Mr. Bibhuti Banerjee in actuator and test rig design, Mr. R. David 
Hampton in the multivariable control theory of active isolation, and Mr. A. Peter Allan in 
actuator design and instrumentation. Several University personnel have worked at NASA 
Lewis during the summers, and contact between NASA and the University has been frequent. 
This final report on the contract reviews the research performed over the three years of the 
grant. Both experimental and theoretical work in microgravity isolation continues at the 
University. The University aspires to become a center of excellence in active vibration 
isolation systems for space applications. 

The next six sections discuss the research efforts of the University in meeting the 
isolation technology needs of the microgravity community. In Section 2, the design and testing 
of a single— degree— of— freedom, large-stroke actuator are reviewed. Section 3 examines 
multiple— degree-of— freedom actuator design. A survey of published designs is presented and a 
new coarse— fine actuator is proposed and analyzed. In Section 4, the design of active isolation 
control systems is examined. Results and discussion of the University’s design methodology, 
frequency-shaped Linear Quadratic Regulator and Kalman-Bucy Filter synthesis, are given in 
Section 5. In Section 6, the design and early experimental results of the microgravity vibration 
isolation test rig are examined. Experimental results demonstrate that active vibration 
isolation of experiments with umbilicals can be obtained using loop— shaped acceleration 
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feedback. Section 7 concludes this final report, with a summary of our results and 
recommendations for future research. An appendix contains all papers presented, published, or 
submitted during the grant. 



